Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,747	0	T	G	42.4%	115.2 / 11.4	66	T255P (ACC→CCC)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (7/31);  new base G (28/0);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 2.29e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646601‑646896                                                                                                                                                    |                                                                                                                                                        tatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCgg                                                                                                                                                     >  1:246416/1‑150 (MQ=255)  atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCg                                                                                                                                                      >  1:246657/1‑149 (MQ=255)  atTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGa‑a                                                                                                                                                   <  1:495743/151‑1 (MQ=255)    ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑At                                                                                                                                                  <  2:421737/150‑1 (MQ=255)    ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGa‑a                                                                                                                                                   >  1:295832/1‑149 (MQ=255)         aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGccccc                                                                                                                                            >  1:895739/1‑151 (MQ=255)            ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCtt                                                                                                                                          <  2:246421/150‑1 (MQ=255)                cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTg                                                                                                                                      <  1:445541/150‑1 (MQ=255)                 aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCa                                                                                                                                    <  2:884426/151‑1 (MQ=255)                  aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTCG‑‑ATGGCCCCttttt                                                                                                                                       >  2:210014/1‑147 (MQ=255)                      cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCTTTTTTCaataa                                                                                                                                >  2:697367/1‑150 (MQ=255)                       ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGG‑ATGCCCCCcttttttgataa                                                                                                                                >  2:62713/1‑142 (MQ=255)                        tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATg                                                                                                                              <  2:289167/150‑1 (MQ=255)                         aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGt                                                                                                                             <  1:34576/150‑1 (MQ=255)                         aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTCTATAAtttt                                                                                                                            >  2:220594/1‑151 (MQ=255)                         aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTCAATAAtttt                                                                                                                            >  2:692715/1‑151 (MQ=255)                           ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCCG‑ATGCCCCCCTTTTGCCATAATAttt                                                                                                                           >  2:178107/1‑150 (MQ=255)                                tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGcc                                                                                                                      <  1:343380/150‑1 (MQ=255)                                 aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGCGG‑ATGCCCCCCTTTTTCAAAAATGtttttgcccc                                                                                                                    >  1:169559/1‑144 (MQ=255)                                       gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTCAAAAATGTTTTCGccccccg                                                                                                                 >  1:903161/1‑143 (MQ=255)                                          cTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGG‑ATGCCCCCTTTTTTTAATAATATTTTTGCCTCGCGCTGt                                                                                                             >  1:279211/1‑149 (MQ=255)                                           tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTGG‑AAGCCCCCCTTTTTTAAAAAAGTTTTTGCCCCCCGCTGt                                                                                                             >  1:398883/1‑148 (MQ=255)                                              aaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGG‑ATGCCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCTGTA‑Gcc                                                                                                        >  2:557892/1‑149 (MQ=255)                                                  cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCa                                                                                                  <  1:756007/151‑1 (MQ=255)                                                  cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGG‑ATGGCCCCCTTTTTTAATAATGTTTTCGCCCCCCGCCGTA‑Gcc                                                                                                        >  2:821337/1‑145 (MQ=255)                                                    ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGG‑ATGCCCCCCTTTTGCAATAATATTTTCGGCTCCCGCCGTA‑GCCCGGGCAgg                                                                                                >  2:92660/1‑151 (MQ=255)                                                    ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTTAATAATGTTTTCGCCTCGCGCCGTA‑GCCCGCGCAg                                                                                                 >  1:466924/1‑150 (MQ=255)                                                        aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTTATTAATGTTTTTGCCCCGCCCTGTA‑GCCCGCGCAGGcccc                                                                                            >  2:867328/1‑151 (MQ=255)                                                           aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCGGGCCCCCCt                                                                                         <  2:28139/151‑1 (MQ=255)                                                           aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCt                                                                                         <  1:560427/151‑1 (MQ=255)                                                           aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTGGA‑ATGCCCCCCTTTTTCTATAATTTTTTTGCCCCCCCCTGTAGGCCACCCAAGGcccccc                                                                                          >  2:501179/1‑151 (MQ=255)                                                             cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGGCCCCTTTTTGCAATAATGTTTTTGCCTCGCGCTGTA‑GCCCGCGCAGGcccccccc                                                                                        >  2:888061/1‑148 (MQ=255)                                                                  ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCAcc                                                                                  <  1:820665/151‑1 (MQ=255)                                                                     cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                               <  2:481517/151‑1 (MQ=255)                                                                     cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACcgc                                                                                <  1:611998/150‑1 (MQ=255)                                                                         aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTTAATAATTTTTTTGCCCCGCGCCGTA‑GCCAGCGCAGGcccccccccccccccccc                                                                              >  2:363606/1‑136 (MQ=255)                                                                           tGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTGGACAAGCCCCCTTTTT‑CAATAATGTTTGCGCCCCTCGCTGTA‑GCCAGGCCCGGCCCCCCTCTCCACCGCCCGATg                                                                          >  2:813684/1‑150 (MQ=255)                                                                            gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATg                                                                          <  2:676324/149‑1 (MQ=255)                                                                            gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCa                                                                        <  1:271353/151‑1 (MQ=255)                                                                                   ggACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGtt                                                                   <  1:62711/149‑1 (MQ=255)                                                                                        gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                             <  2:518353/150‑1 (MQ=255)                                                                                          gaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGG‑CTGCCCCCCTTTTTTTATAATTTTTTTGCCCCCCGCGGTG‑GCCCGCGCGCGCCCCCCCCCCCCCCCCGCGATGCCAccgtggggt                                                              >  1:880868/1‑139 (MQ=255)                                                                                           aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTCAATAAAGTTTTCGCCTCCCGCTGTA‑GCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGTGCAACGTTTGTATCGcc                                                          >  2:838765/1‑150 (MQ=255)                                                                                           aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCCG‑ATGCCCCCCTTTTTCACTAATGTTTTTGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCCCCCCCCCCCGCGACGCAACGTTGGGATCGcc                                                          >  1:221868/1‑150 (MQ=255)                                                                                             ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGt                                                        <  1:206931/150‑1 (MQ=255)                                                                                                  tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                                   <  1:383894/150‑1 (MQ=255)                                                                                                       acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATc                                             <  2:857697/151‑1 (MQ=255)                                                                                                         acTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGG‑AAGCCCCCCTTTTTCAATAATTTTTTCGCCTCCCGCTGTA‑GCCAGCGCAGGCCCCCCCCCCCACCCCCCGCTGCAACGTTGGGATTGCCGGTTATCGcc                                                >  2:294332/1‑146 (MQ=255)                                                                                                            cccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTa                                         <  1:788704/150‑1 (MQ=255)                                                                                                             ccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGC‑ATGGCCCCCTTTTTTAATAATGTTTTTGCCCCCCGCTGGA‑GCCACCGCAGGGCCCCCCCCCCCCCCCGCCGTGCAAAGGTGGGGACGCCCGTGATCGCCATCCGCAgg                                       >  1:202272/1‑151 (MQ=255)                                                                                                                gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                     <  2:86992/150‑1 (MQ=255)                                                                                                                gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGC‑ATGCCCCCCTTTTTCAATAATGTTTTCGCCTCCCGCTGTA‑GCCAGCGCCGGCCCCCCCCCCCCCCCCCCGGTGCCACGGTGGGATCGCCCTTGATCGCCCTCAATAgg                                       >  2:906490/1‑148 (MQ=255)                                                                                                                  cGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                                   <  2:645070/150‑1 (MQ=255)                                                                                                                    aaCTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAgg                                <  1:491255/151‑1 (MQ=255)                                                                                                                      cTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGA‑ATGCCCCCCTTTTTCCATAATATTTTTGCCCCCCGCTGGA‑GCCAGCGCCGGGCCCCCCCCCCCCCCCCCGACGCAACGGTGGTATCGCCGGTGGTCGCCATCAGGAGGGGGAAggg                               >  1:685569/1‑149 (MQ=255)                                                                                                                         ccGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCg                           <  2:886277/151‑1 (MQ=255)                                                                                                                          cGGAGATAATCGAGATCGGC‑GGGGGGTTGA‑ATGCCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCTGTA‑GGCCGCGCCGGGCCCCCCCCCCCCCCCCCGATGCCACATTGGGATCGCCCGTGATCGCCCTCAGTAGGAGGAAgg                                >  2:371232/1‑145 (MQ=255)                                                                                                                              gaTAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCaa                      <  2:192470/151‑1 (MQ=255)                                                                                                                                         tCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGGATCGAGCAATGCCAGTTc             >  1:645390/1‑149 (MQ=255)                                                                                                                                           ggc‑ggGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCGAGTTCAGGa         <  1:946058/151‑1 (MQ=255)                                                                                                                                            gc‑ggGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt        <  2:895757/151‑1 (MQ=255)                                                                                                                                               gggggTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc       <  1:451558/150‑1 (MQ=255)                                                                                                                                               gggggTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      >  2:242309/1‑151 (MQ=255)                                                                                                                                               gggggTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      <  2:237484/151‑1 (MQ=255)                                                                                                                                                ggggTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      >  2:793597/1‑150 (MQ=255)                                                                                                                                                   gTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc  >  1:893763/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                        TATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCGA‑ATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTA‑GCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646601‑646896

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑