Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,747	0	T	G	52.2%	65.4 / 110.9	67	T255P (ACC→CCC)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (2/30);  new base G (35/0);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 4.97e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.01e-12

AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAC  >  NC_000913/646599‑646895                                                                                                                                                     |                                                                                                                                                     aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtt                                                                                                                                                     <  2:915370/150‑1 (MQ=255) aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGtt                                                                                                                                                     <  1:351788/150‑1 (MQ=255)     tttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTccg                                                                                                                                                  >  2:282670/1‑147 (MQ=255)      ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAt                                                                                                                                                <  2:235978/150‑1 (MQ=255)        ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGc                                                                                                                                              <  2:556247/150‑1 (MQ=255)           aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATAccccc                                                                                                                                          >  1:724836/1‑151 (MQ=255)           aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATGgcccc                                                                                                                                          >  2:218416/1‑151 (MQ=255)            aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGccccc                                                                                                                                          >  1:772417/1‑150 (MQ=255)            aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGccccc                                                                                                                                          >  1:968690/1‑150 (MQ=255)              ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGccccc                                                                                                                                          >  1:769207/1‑148 (MQ=255)               tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGACGAATGGACCCtttt                                                                                                                                      >  2:709378/1‑151 (MQ=255)                 tCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCAAACCCCCctttt                                                                                                                                     >  2:658828/1‑150 (MQ=255)                  cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTg                                                                                                                                    <  1:105956/150‑1 (MQ=255)                  cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGc                                                                                                                                   <  2:125613/151‑1 (MQ=255)                   aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCa                                                                                                                                  <  2:998936/151‑1 (MQ=255)                         ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCtataa                                                                                                                              >  1:349117/1‑149 (MQ=255)                             ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGATGGCCCCCTTTTTCAATAAtttt                                                                                                                          >  2:88364/1‑149 (MQ=255)                                 ttAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCTTTTTTCAATAATGTTTTCGcc                                                                                                                    >  1:396313/1‑151 (MQ=255)                                     aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTTGCCCcgc                                                                                                                >  1:303911/1‑151 (MQ=255)                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTg                                                                                                            <  1:695699/151‑1 (MQ=255)                                         gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGGAGTCCCCTTTTTGTTATAATATTTTCGCCTcg                                                                                                                 >  1:517700/1‑146 (MQ=255)                                          tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGCCCCCCTTTTGCCATAATGtttttgccccccgc                                                                                                              >  2:720054/1‑137 (MQ=255)                                           gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGt                                                                                                           <  2:301216/150‑1 (MQ=255)                                           gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATTTTTTTGCCTCGCGCTg                                                                                                            >  1:979082/1‑149 (MQ=255)                                              aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGcc                                                                                                       <  2:91540/151‑1 (MQ=255)                                                 aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgc                                                                                                    <  2:484594/151‑1 (MQ=255)                                                      ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCAATAAAGTTTTCGCCCCGCGCCGTAGCCAgcg                                                                                                   >  1:946920/1‑147 (MQ=255)                                                      ggTAAGGATCAGTAGGTCAGCACTGCCGACTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTTGCTGTCCCCGTTTTTTAAAAACGTTTTCTCGGCCCCCTGTAGCCAGCCCAg                                                                                                >  1:558427/1‑150 (MQ=255)                                                       gTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGc                                                                                              <  2:10037/151‑1 (MQ=255)                                                            gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCTTTTTTGAATAATTTTTTCGCCCCGCGCTGGAGCCAGCGCAGGgccccc                                                                                         >  1:139090/1‑151 (MQ=255)                                                             aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGCAAGCCCCCCTTTTTCCATAATGTTTTCGCCCCCCGCCGTAGCCAGCGCCGGgcccccc                                                                                        >  1:127110/1‑150 (MQ=255)                                                             aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATGCCCCCCTTTTTTAATAAAGTTTTTGCCCCCCCCGGTAGCCCCCGCAGGccccccc                                                                                        >  2:831836/1‑150 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTc                                                                                       <  2:927407/150‑1 (MQ=255)                                                               cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATTCCCCCCTTTTTCTATAAATTTTTCGCCCCCCGCCTTAGCCAGCGCAGGcccccccc                                                                                       >  1:934759/1‑148 (MQ=255)                                                                  tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAAATTTTTCGCCCCGCGCGGTTGCCAGCGCCGGcccccccccccc                                                                                   >  2:848115/1‑145 (MQ=255)                                                                     gTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcg                                                                                <  1:365629/151‑1 (MQ=255)                                                                     gTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcg                                                                                <  1:177926/151‑1 (MQ=255)                                                                      tCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGATGGCCCCTTTTTTCAATAAAATTTTTCCCTCGCCCTGTGGGCAACGCAGGcccccccccccccccc                                                                               >  1:607568/1‑141 (MQ=255)                                                                       cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGGCCCCCTTTTTTAATAAAATTTTTGCCCCCCCCTGTTGCCCCCGCACGcccccccccccccccc                                                                               >  2:831223/1‑140 (MQ=255)                                                                         gCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGCTGCCCCCCTTTTTGAATAATGTTTTTCCCCCGCGGAGTTGCAAGCGCCGGCCCCCCCCGCcccc                                                                                 >  2:263864/1‑146 (MQ=255)                                                                          cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcgcg                                                                            <  1:125719/150‑1 (MQ=255)                                                                            cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTCGATAAAGTTTTCGCCTCGCCCTGTAGCCCGCGCAGGCCCCCCCCGCCACCCCGCGGTg                                                                         >  1:104593/1‑151 (MQ=255)                                                                            cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGGCCCCCTTTTTCAATAAATTTTTCGCCCCGCGCCGTAGCCAGCGCAGGCCCCCCCCGCCCCCCCGCGGTg                                                                         >  1:613157/1‑151 (MQ=255)                                                                                cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAATTACCCCCTTTTTTAATAATTTTTTTGCCCCGCGCCGTAGGCAGCCCCGGCCCCCCTCcccccccccc                                                                             >  2:729287/1‑139 (MQ=255)                                                                                 gccTGGACTGAGATTTCGTTCGATACTCTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACg                                                                    <  2:418777/151‑1 (MQ=255)                                                                                    tGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCGAGTGCCCCCTTTTTTCAATGATGTTTGTCGGGCGCGCTGTAGGCGGCGCCCGCCCCCCCCCCCCCCCCGCGGGGCAACggt                                                                  >  2:429498/1‑148 (MQ=255)                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                            <  2:793565/150‑1 (MQ=255)                                                                                          gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCg                                                           <  1:21850/151‑1 (MQ=255)                                                                                               ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGATGCAAAGTTGGTATCGcc                                                         >  1:80131/1‑148 (MQ=255)                                                                                                 cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTg                                                     <  2:853609/150‑1 (MQ=255)                                                                                                   ttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTTAATAATTTTTTTGCCTCGCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGGGCAACGGTGGTATc                                                            >  2:577743/1‑141 (MQ=255)                                                                                                    tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGAATGCCCCCCTTTTTTAATAATATTTTTGCCCCGCGCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGATGGAACGGTTGGATTGCCGGTg                                                     >  1:491061/1‑147 (MQ=255)                                                                                                    tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAAGCCCCCCTTTTTCCATAATTTTTTCGCCCCGCGCGGTAGGCAGCGCAGGCCCCCCCCGCCACCCCGCGAAGCAAAGGTGGGATTGCCGGTggt                                                   >  1:784547/1‑147 (MQ=255)                                                                                                        ttcgcgcccGGTCGAACGGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCGCCGCGCGATGCAACGTTGGTGTCGCCGTTGATCg                                                 <  2:183612/141‑1 (MQ=255)                                                                                                        tACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAATGCCCCCCTTTTTCGATAATATTTTTGCCCCGCGCCGGAGCCCGCGCCGGCCCCCCCCCCCCCCCCCCGGTGCAACGCTGGGGTCGCCCTTGTTCGc                                                >  2:9329/1‑148 (MQ=255)                                                                                                            cTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAg                                          <  2:261861/150‑1 (MQ=255)                                                                                                             tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGCATGCCCCCCTTTTTCAATAATGTTTTCGCCTCGCGCTGTAGCCAGCGCAGGGCCCCCTCGCCACCCCGCGCGGCGACGTTGGGATCGCCGTTGAtag                                                 >  1:865459/1‑140 (MQ=255)                                                                                                               ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAgg                                      <  2:270840/151‑1 (MQ=255)                                                                                                                  gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                    <  2:740859/150‑1 (MQ=255)                                                                                                                       aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGt                              <  1:332179/151‑1 (MQ=255)                                                                                                                              gagaTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                       <  2:968605/151‑1 (MQ=255)                                                                                                                                 aTAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAAt                    <  2:127099/151‑1 (MQ=255)                                                                                                                                            cggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAgg         <  2:330954/151‑1 (MQ=255)                                                                                                                                             ggcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGa        <  2:214427/151‑1 (MQ=255)                                                                                                                                             ggcgaGGGTTCGGATGACCCCTTTTTCCAATAATGGTTGCGCCTCGCGCTGCAGCCAGCGGCGGCCCCCCACGCCACCGCGCGATGCAACGCTTGTATCCCCGATGATCGCCATCAGTAGGAGCTAGGTATCGAGCCATGCCAGt              >  2:790351/1‑145 (MQ=255)                                                                                                                                              gcggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt       <  2:772332/151‑1 (MQ=255)                                                                                                                                               cggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc      <  1:88378/151‑1 (MQ=255)                                                                                                                                               cggGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc      <  2:971583/151‑1 (MQ=255)                                                                                                                                                   ggTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAc  >  1:69254/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAC  >  NC_000913/646599‑646895

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