Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	650,454	0	T	G	41.9%	98.7 / 11.4	63	T13P (ACC→CCC)	citE	citrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base T (10/26);  new base G (26/0);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.30e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAC  >  NC_000913/650313‑650603                                                                                                                                              |                                                                                                                                                         aaaCCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGtttt                                                                                                                                                 <  1:1101161/151‑1 (MQ=255)   aCCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTAc                                                                                                                                               <  1:374470/151‑1 (MQ=255)      aTGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCgggcgttgcgtt                                                                                                                                            >  2:1154264/1‑140 (MQ=255)            cgggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCa                                                                                                                                      <  1:442894/151‑1 (MQ=255)             gggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGATTTTCggtgtt                                                                                                                                         >  1:654449/1‑147 (MQ=255)                    gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGc                                                                                                                               <  2:9656/150‑1 (MQ=255)                        ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGaa                                                                                                                          <  2:617487/150‑1 (MQ=255)                        ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTTGTTGTTGCAGCGAAGA‑ggga                                                                                                                          >  1:303594/1‑148 (MQ=255)                          ttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAAt                                                                                                                        <  1:1034865/150‑1 (MQ=255)                            cACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTTGTTGTTTCAGCGAAGC‑GGAAATCa                                                                                                                      >  2:779658/1‑150 (MQ=255)                            cACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGCAGCGAAGC‑GGGAATc                                                                                                                       >  1:526888/1‑149 (MQ=255)                              cgcgATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑tt                                                                                                                   <  2:482487/150‑1 (MQ=255)                                 aaTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCa                                                                                                               <  1:154191/151‑1 (MQ=255)                                  aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCaa                                                                                                              <  1:706011/151‑1 (MQ=255)                                       aCGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAAtcctc                                                                                                         <  2:77273/151‑1 (MQ=255)                                          gAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCa                                                                                                      <  1:602212/151‑1 (MQ=255)                                           aGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGAGTTTTTCGTTGTTGCAGCGAAAC‑GGAAAACA‑TTGGCAATCCTccc                                                                                                       >  1:522003/1‑149 (MQ=255)                                             tattCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGCGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGCGGGGGGTTTTCGGTGGTTCAGCGAAGC‑GGAAACAA‑TTGGGAATACTacc                                                                                                       >  2:1179134/3‑147 (MQ=255)                                             tCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCGTTGTTGCCGAGAAGA‑GGGAATCA‑TTGGGAATCCCCCCATgg                                                                                                   >  2:812223/1‑151 (MQ=255)                                              cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGGTTTTGGAGCCGAGC‑GGAAATCA‑TTGGCACTCCTccccttgc                                                                                                  >  2:1059274/1‑146 (MQ=255)                                                tCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCa                                                                                                 <  1:1119800/150‑1 (MQ=255)                                                 cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTAGCAATCCTCCCATGGCag                                                                                                <  1:738814/150‑1 (MQ=255)                                                  gAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGGTTTCGTTTTTGCAGCGAAAC‑GGAAAATA‑TTGGCAAACCTCCCaacg                                                                                                   >  2:1224145/1‑143 (MQ=255)                                                     gTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGTTTTTCGTTGTTGCAGCGAAGC‑GGGAATCA‑TTGGCAATCCTCCCCTGGGAGAGcc                                                                                            >  1:1247438/1‑150 (MQ=255)                                                        aaaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGGAGCGAAGC‑GGGAATCA‑TTGGGCAACCTCCCATGGCAGAGCCggg                                                                                         >  1:833399/1‑150 (MQ=255)                                                         aaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGAGGTTTTCGTTTTTGGAGCGAAGC‑GGAAAACC‑TTTGGAATCCTCCCCTCGCagaggccgg                                                                                         >  2:788171/1‑144 (MQ=255)                                                               aGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGTTTTAGTGTGTTGGAGCGAAGG‑GGGAATCC‑TTTGGAATTCCCCCCTGGCAGAGCCGGGa                                                                                        >  2:398778/1‑144 (MQ=255)                                                                  gCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGgc                                                                             <  1:60241/151‑1 (MQ=255)                                                                      cAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCGTTGTTGCAGCGAAAC‑GGAAAACA‑TTTGCAATCCTCCCCCTGCAGAGCCGGGGTACCG‑CTGgcg                                                                            >  1:360499/1‑148 (MQ=255)                                                                         ccGGGTATATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGCAGGTTTTTTCTGTTGCAGCGGAGC‑GGAAATTC‑TTTGCCATACCCCCATGGCAGAGACGGGGTAACG‑CTGGCACGTg                                                                        >  2:1295029/1‑149 (MQ=255)                                                                          cGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCca                                                                     <  1:862809/151‑1 (MQ=255)                                                                          cGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCca                                                                     <  2:832820/151‑1 (MQ=255)                                                                                 aTGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTACGTTGTTTCAGCGAAGC‑GGGAATCA‑TTGGCAAACCTCCCCTGGCAGAGCCGGGGTACCG‑CTGGCGCGGGCCAGCAggg                                                               >  1:1349948/1‑150 (MQ=255)                                                                                  tGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTTGTTGTTGGAGCGAAGC‑GGAAAACA‑TTTGCAATCCTCCCATGGCAGAGCCGGGATACCC‑CCGGCGCGGGCCCGCAggg                                                               >  1:1080007/1‑149 (MQ=255)                                                                                  tGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTACGTTGTTTCCGCGAAGG‑GGGAAAAA‑TTGGGAATCCTCCCCTGGCAGAGCCGGGATAACG‑CTGGCCGGTGCCAGCAggg                                                               >  1:1252254/1‑149 (MQ=255)                                                                                  tGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGAGGTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCCTGGCAGAGCCGGGATAACG‑CTGGCGCGGGCCAGCAggg                                                               >  1:601221/1‑149 (MQ=255)                                                                                    aaGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGGTGTTGCAGCGAAGC‑GGAAATAA‑TTGGCCATCCTCCCATGGCAGAGCCGGGGTACCC‑CCGGCGCGCGCCCGCCggg                                                               >  2:264164/1‑147 (MQ=255)                                                                                      tgAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGTTTTTTGTTGTTGGAGCGAAGC‑GGAAAACA‑ATGGCAACCCACCCCTCGCAGAGACGGGGGACCG‑CCGGCGCGGGCCCGCAGGGCTTccc                                                         >  2:234956/2‑150 (MQ=255)                                                                                       gAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAg                                                        <  1:1336018/151‑1 (MQ=255)                                                                                       gAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTTTTGCAGCGAAAA‑GGAAAAAATTTGCCAACCCCCCCTTGGC‑GAGCCGGGGTACCG‑CCGGGGCGTGCCCCCCGGGCCTccc                                                         >  2:371243/1‑149 (MQ=255)                                                                                            gCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTCGTTTTTTCAGCGAAGC‑GGGAAACA‑TTTGCAAACCTCCCCTGGCAGAGCCGGGGTACCG‑CTGGCGCGCGCCCGCCGGGCGTCCCGt                                                       >  1:643362/1‑147 (MQ=255)                                                                                             cTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGAGGTTTTGGTTGGTGCAGCGAAGCTGGAAACAC‑TTGGC‑ACCCTCCCACGGCCGGGCCCGGCGCCCGCCGGGCCCGGGCC‑GCAGGGCTTCCAGTCGTg                                                   >  2:470189/1‑150 (MQ=255)                                                                                               gACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGGTTTACGTTGTTGCAGCGAAAC‑GGAAATCA‑TTGGGAATCCTCCCATTGCAGAGCCGGGGTAACG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCAc                                                >  1:659932/1‑151 (MQ=255)                                                                                                 ccATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTCGGTGTTGCAGCGAAGC‑CGAAAACA‑TTGGCAAACCTCCCCTGGCAGAGCCGGGGTAACG‑CTGGCGCGGGCCAGCAGGGCCTCCCGGCGt                                                    >  2:843977/1‑145 (MQ=255)                                                                                                    tCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGCAGCGAAAC‑GGGAAACA‑TTGGGAATCCTCCCCTGGCAGAGGCGGGGGACCG‑CTGGGGCGGGCCCGCAGGGCTTCCCGGCGTGGACGTaa                                            >  1:914806/1‑150 (MQ=255)                                                                                                      gcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGGGGTTTTCGTTTTTGCAGCGAAAA‑GGAAAACC‑ATGGCAAACCTCCCCTGGCAGAGACCGGGGACCC‑CTGGGGCGGGGCCGCAGGGCTCCCCGCCGTGccccgaaaa                                          >  2:955750/1‑142 (MQ=255)                                                                                                      gcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTTGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCCATGCAGAGCCGGGGTACCG‑CTGGCGCGGGGCCGCAGGGCTTCCAGGCGTGCACGTaaaa                                          >  1:1086203/1‑150 (MQ=255)                                                                                                           aTTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGt                                    <  1:1137646/151‑1 (MQ=255)                                                                                                              ggCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCa                                 <  1:1018582/151‑1 (MQ=255)                                                                                                                cACCAGGCACAAACAACATGCTGCGGCGGGGGCGAGTTTTACGTTGTTGCCGCGAAAC‑GGAAAACA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCC‑CTGGCGCGTGCCAGCCGGGGTTCCCGTCGTGCACGTAAAAAGCAGTccc                                 >  1:184138/1‑148 (MQ=255)                                                                                                                 aCCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTg                              <  1:472675/151‑1 (MQ=255)                                                                                                                 aCCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTg                              <  1:1297888/151‑1 (MQ=255)                                                                                                                  ccAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTACGTTGTTGCAGCGAAGC‑GGAAATAA‑TTGGCAATCCTCCCATGGCAGAGCCGGGGTACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTg                              >  1:1340093/1‑150 (MQ=255)                                                                                                                   cAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTgcg                            <  1:367759/151‑1 (MQ=255)                                                                                                                    aGGCACAAACAACATGCTGCGGCGGGTGCGAGGTTTTCGTTTTTGCAGCGAAAA‑GGAAATCA‑TTGGGAAACCCCCCCTGGCAGAGCCGGGGTACCC‑CTGGGGCGGGCCAGCAGGGCATCCAGTCCTGGACGTAAAAAGCAGTCCAGGgcgc                           >  2:1202748/1‑151 (MQ=255)                                                                                                                       cacaAACAACATGCTGCGGCGGGGGCGGGTTTTTTGTTGTTTCCGCGAAGC‑GGAAATCA‑TTTGCAATCCTCCCCTGGCAGAGCCGGGATACCC‑CTGGCGCGTGCCCGCCGGGCGTCCCGTCGTGCCCGGAAAAAGCAGTCCAGTGCGcc                          >  1:410848/1‑149 (MQ=255)                                                                                                                       cacaAACAACATGCTGCGGCGGGGGCGAGGTTTACGTTGTTGCCGCGGAGC‑GGAAAACA‑TTGGGAAACCTCCCATGGGAGAGGCGGGATACCG‑CCGGCGCGTGGCAGCCGGGGTTCCCGTCGTGCACGTAAAAAGCAGGCCAGTGCGcc                          >  1:1212902/1‑149 (MQ=255)                                                                                                                        acaaacaaCATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCttt                       <  2:745622/151‑1 (MQ=255)                                                                                                                              aacaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCAt                  <  1:192936/150‑1 (MQ=255)                                                                                                                               acaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCg                <  2:107444/151‑1 (MQ=255)                                                                                                                                      gcggcgGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCa         <  2:526770/151‑1 (MQ=255)                                                                                                                                          cggGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTg     >  1:284438/1‑151 (MQ=255)                                                                                                                                             gTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAc  >  1:447314/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                                         AAACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAC  >  NC_000913/650313‑650603

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