Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	650,464	0	A	T	41.3%	81.6 / 24.9	63	R9R (CGT→CGA)	citE	citrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (11/26);  new base T (26/0);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 2.45e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGT  >  NC_000913/650315‑650614                                                                                                                                                      |                                                                                                                                                          aCCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTAc                                                                                                                                                          <  1:374470/151‑1 (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                                                                                <  2:832820/151‑1 (MQ=255)                                                                        cGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCca                                                                                <  1:862809/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:1297888/151‑1 (MQ=255)                                                                                                               aCCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTg                                         <  1:472675/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)                                                                                                                                              cGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgcc          >  1:438654/1‑151 (MQ=255)                                                                                                                                                     tACGTTGTGGGCGCGAAGC‑GGGGATCA‑TGTGCGATCCTCCCATGGCAGAGCCGGGATGCCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCGCGTAAAATGCAGTCCAGTGGGCCTTTGTCAGCGACATTCAGCTGTACGGCGCGCAc    <  2:550689/150‑1 (MQ=255)                                                                                                                                                      aCGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACg   <  2:579263/150‑1 (MQ=255)                                                                                                                                                      aCGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGt  <  2:1004212/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                          ACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGC‑GGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCG‑CTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGCACGT  >  NC_000913/650315‑650614

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