Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,628,625	0	A	C	23.5%	155.3 / 26.1	68	V105G (GTC→GGC)	pdeF	cyclic di‑GMP phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base A (47/5);  new base C (0/16);  total (47/21)
Fisher's exact test for biased strand distribution p-value = 1.38e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628487‑2628764                                                                                                                                            |                                                                                                                                               cTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGAc                                                                                                                                               >  2:1223679/1‑140 (MQ=255) cTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGAGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGAc                                                                                                                                               >  1:351815/1‑140 (MQ=255)   cGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAG‑ATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCt                                                                                                                                            >  2:415330/1‑140 (MQ=255)         tAATCCTACAAATGCAGCAAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATACGCGGGGAAATTAAACGGGTATCCCCAT‑GTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGT‑AAc                                                                                                                                      <  2:1844343/140‑1 (MQ=255)         aaaTCCTACAAATGCAGCAAACGGGAAAAGTATCAAAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGAGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AAc                                                                                                                                      <  2:351815/139‑1 (MQ=255)             ccaaaaaaTGCAGAAAACTGAAAAAGTATCAGAAACAGCGTTGAAGAAAACACAATCTCCAAGAATATACGCGGGAAAATTAAACGGGTACCCCCAT‑GTGAAACGTTGTTCCGCCGGGGAGAAAAGCCCCTGT‑AACCGcc                                                                                                                                  <  2:1144062/135‑1 (MQ=255)                aCAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCa                                                                                                                               >  2:438683/1‑140 (MQ=255)                   aaggcagaaAACTGAAAAAGTAAAAAAAACAGCTTGGAAGAAAACACAATCCGCAAAAATATCCGGGGGAAAATAAAACGGTTACCACCAT‑GTAAACCTTTGTTCCGCGGGGGAGAAAAGCCCCTGT‑AACCGCCCCAACa                                                                                                                            <  2:1079181/137‑1 (MQ=255)                   aacgcagcaAGATGGAAAAAAAAGAGAAACCGCGTTGGAAGAAACACAATCTGCCAGAATAGGCGCGGGGAAATAAAAGGGGTATCCCGTT‑GTGAACCGTTGTTCCGCCGGGGAATAAAGCCCTTTT‑AACCGCCCCAACa                                                                                                                            <  2:301291/137‑1 (MQ=255)                   aaTGCAGCAAACTTGAAAAGAATCAGAAACAGCGTTGCGGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACTCTGT‑AACCGCCCCAACa                                                                                                                            >  1:590024/1‑140 (MQ=255)                     tGCAGAAAACTGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAAAATATCCGCTGGAAAATTAAACGGTTACCCCCAT‑GTAAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGT‑AACCGCCCCAACaga                                                                                                                          <  2:260180/140‑1 (MQ=255)                       cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAgaga                                                                                                                        >  1:2004631/1‑140 (MQ=255)                       aagcaAAGTGGAAAAGAATAAAAACCAGCGTTGGAGAAAACACTATCGGCCAAAATACGCGTGGGAAAATTAAACGGGTCACCCCTT‑GGAAATCGTTGTAACGCCGGGGAGTAAAGCCCCTGT‑AACCGCCCCAACAgaga                                                                                                                        <  2:1099773/139‑1 (MQ=255)                         gcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGaat                                                                                                                      >  1:816740/1‑140 (MQ=255)                          caAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGaatg                                                                                                                     >  1:1133465/1‑140 (MQ=255)                            aaCTGAAAAAGTATCAAAAACAGCGTAGCAGAAACCACAATCTGTAAGAATATACGCGGGAAAATAAACCGGTTATCCCCAT‑GTGAAAGTTTGTTCCGCGGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGaatgat                                                                                                                   <  2:717353/140‑1 (MQ=255)                            aaCTGAAAAAGTATAAAAAACAGCGTGGCAGGAAACACAATCTGCCAAAATATCCGCTGGAAAATAAAACGGGTATCACCAT‑GTAAAACGTTGTTCCGCGGGGGAGAAAAGCCCCTGT‑AACCGCCCCAACAGAGaatgat                                                                                                                   <  2:1419095/140‑1 (MQ=255)                               tGGAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGT‑AACCGCCCCAACAGAGaatgataat                                                                                                                <  2:781759/140‑1 (MQ=255)                                ggAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGaatgataact                                                                                                               >  1:592138/1‑138 (MQ=255)                                   aaaGTATCAAAAACAGCGTTGCAGAAAACAAAATCTGCAAGAATATACGCTGGAAAATTAACCGGGTACCCCCAT‑GTAAAACTTTGTTACGCCGGGGAGAAAAGCCCCTGT‑AACCGCCCCAACAGAGAatgataatgata                                                                                                            <  2:816740/140‑1 (MQ=255)                                           aaaaaCAGTGTTGCACGAAACCAAATTTCCAAAAATATACGTGGAAAAATTAACGGGGTATCACCAG‑ATGAAAGGTTTTCACGCGGGGGAAAAAAGCCCCGGT‑AGCCGCCCCAACAGAGAATGATAATGATAATAAAAtt                                                                                                    <  2:1133465/138‑1 (MQ=255)                                             aaaCAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCGGGGGAGTAAAGCCCCTGC‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                  <  1:142307/140‑1 (MQ=255)                                             aaaCAGCGTTGCAGAAAACACAATCTGCCAGAATATCCGCTGGAAAATAAAACGGGTATCCCCAT‑GTAAAACTTTGTTACGCGGGGGAGTAAAGCCCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                  <  2:2054943/140‑1 (MQ=255)                                             aaaCAGCGTTGAAGAAAACACAATCTGCCAAAATATACGCTGGAAAATAAAACGGTTATCACCAT‑GTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                  <  2:116469/140‑1 (MQ=255)                                               aCAGCGTGGCAGAAAAAAAAATCTGCAAAAATATCCGCGGGAAAATAAAAGGGGTATCCCCAT‑GTAAACCTTTTTTCCGCGGGGGAGAAAAGCCCTTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGcc                                                                                                <  2:2075627/140‑1 (MQ=255)                                                cAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg                                                                                               >  1:970363/1‑140 (MQ=255)                                                cAGCGTTGCAGAAAACACAATCTGCAAGAATAAACGTTGGGAAATAAAAGGGGTATCACAAT‑GTGAAACGTTTTTACGCCGGGGAGAAAAGCCCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCg                                                                                               <  2:570868/140‑1 (MQ=255)                                                 acagTTGAAGAAAACAAAATCTGCAAAAAAATCCGCTGGAAAATAAAACGGGTATCCCCAT‑GTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                              <  2:1508794/137‑1 (MQ=255)                                                 aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                              >  1:935272/1‑140 (MQ=255)                                                   cGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCC‑ACAGAGAATGATAATGATAATAAAATTAGCCGTTaa                                                                                           >  1:1343263/1‑140 (MQ=255)                                                       gaagaaaaCAAAATCTGCAAAAATATACGCGGGAAAATAAACCGGGTACCACAAT‑GTAAACCTTTTTTCCGCGGGGGAGTAAAGCCCCTTA‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑a                                                                                       <  2:1572703/135‑1 (MQ=255)                                                           gAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATAATAATAAAATTAGCCGTTAATG‑ATAGc                                                                                   >  1:2015321/1‑140 (MQ=255)                                                             aaCACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGt                                                                                 >  2:1143478/1‑140 (MQ=255)                                                               cacaATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGttt                                                                               <  1:464250/140‑1 (MQ=255)                                                                acaATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCgctgt                                                                              >  1:758516/1‑136 (MQ=255)                                                                acaATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCgctgt                                                                              >  1:1508794/1‑136 (MQ=255)                                                                  aaTCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGt                                                                            >  1:1916379/1‑140 (MQ=255)                                                                    ccgaccATAATATTCGC‑GGGAAAATAAAGCGGTAATCCCATGGTGAAAC‑TTGTTGCCGCAGGGAATAAAACCCCCGTAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTaa                                                                          <  2:828571/136‑1 (MQ=255)                                                                      tGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTaaaa                                                                        >  1:1096275/1‑140 (MQ=255)                                                                         cAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAAcc                                                                     >  2:60142/1‑140 (MQ=255)                                                                                taCGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGctct                                                              >  2:1459038/1‑140 (MQ=255)                                                                                  cGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGc                                                            >  1:2138993/1‑140 (MQ=255)                                                                                     tGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCata                                                         >  1:618905/1‑140 (MQ=255)                                                                                     tGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCata                                                         >  1:1205378/1‑140 (MQ=255)                                                                                      gggAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                        >  1:492721/1‑140 (MQ=255)                                                                                       ggAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                       >  2:390820/1‑140 (MQ=255)                                                                                       ggAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                       >  1:1283838/1‑140 (MQ=255)                                                                                       ggAAATTAAACGGGGATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGACTGATAATGATAATAAAATTAGCAGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                       >  2:1486712/1‑140 (MQ=255)                                                                                                 ggggTATCCCCAT‑GAAAAATTTTTTTCCGCGGGGGAGAAAAGCCCCTGA‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                             <  2:761277/139‑1 (MQ=255)                                                                                                 ggggTATCCCAAT‑GTGAACCGTTTTTCCGCCGGGGAGTAAAGCCCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                             <  2:1341642/139‑1 (MQ=255)                                                                                                 cGGGTATCCCCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa                                             <  1:1443095/140‑1 (MQ=255)                                                                                                  ggttATCACCAT‑GTAAAAGGTTTTTCCGCGGGGGAGAAAAGCCCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaa                                            <  2:1283838/137‑1 (MQ=255)                                                                                                   ggTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaaa                                           >  2:688587/1‑140 (MQ=255)                                                                                                     tATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACC‑CCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATcc                                        >  1:894020/1‑140 (MQ=255)                                                                                                      aTCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATcc                                        >  2:1814836/1‑140 (MQ=255)                                                                                                      aTCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATcc                                        >  1:363676/1‑140 (MQ=255)                                                                                                       tCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATccc                                       >  1:1685995/1‑140 (MQ=255)                                                                                                            aT‑GTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                  <  2:935272/140‑1 (MQ=255)                                                                                                            aT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg                                  <  1:60142/140‑1 (MQ=255)                                                                                                             t‑gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATGATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑c                                >  1:2035583/1‑140 (MQ=255)                                                                                                             t‑gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑c                                >  1:761277/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGTGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑Cg                               >  1:64126/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑Cg                               >  1:1088144/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑Cg                               >  2:941861/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑Cg                               >  1:709326/1‑140 (MQ=255)                                                                                                               gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑Cg                               >  2:1178625/1‑140 (MQ=255)                                                                                                                tgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGa                              >  2:736191/1‑140 (MQ=255)                                                                                                                   aaCGTTGTTACGCGGGGGAGAAAAGCCCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATcc                           <  2:1685995/140‑1 (MQ=255)                                                                                                                      gttgttACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTg                        >  1:1384840/1‑140 (MQ=255)                                                                                                                         tttccgcgGGGGAGAAAAGCCCCGGA‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCa                     <  2:492721/136‑1 (MQ=255)                                                                                                                         gttACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCa                     >  1:1870119/1‑140 (MQ=255)                                                                                                                         gttACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCa                     >  1:363753/1‑140 (MQ=255)                                                                                                                          ttACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCaa                    >  1:81519/1‑140 (MQ=255)                                                                                                                          ttACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCaa                    >  1:57647/1‑140 (MQ=255)                                                                                                                                  gggAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCAATGCCGccc            >  2:474748/1‑140 (MQ=255)                                                                                                                                  gggAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCAATGCCGccc            >  1:61308/1‑140 (MQ=255)                                                                                                                                  gggAGAAAAGCCCCTGA‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCAATGCCGccc            <  2:709326/140‑1 (MQ=255)                                                                                                                                     agaaaaGCCCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCAATGCCGCCCAGc         <  1:114253/137‑1 (MQ=255)                                                                                                                                     aGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCAATGCCGCCCAGc         <  1:1459038/140‑1 (MQ=255)                                                                                                                                       tAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCAATGCCGCCCAGCTg       <  2:1096275/140‑1 (MQ=255)                                                                                                                                        aaaGACCCTGT‑AACC‑CCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCAATGCCGCCCAGCTGaa     <  2:894020/140‑1 (MQ=255)                                                                                                                                         aaGACCCTGT‑AACCGCCCCAACAGAGAATGATAA‑GGTAATAAAAATAACCGGTAATGAAAACCGTTTGGAAAAAACCCAGCTTGCCAAATTGCCGGAC‑AAAATCCCCAAGGCCAACCCCTGCCAATGCGCCCCAGCTGaa     >  2:762528/1‑140 (MQ=255)                                                                                                                                            aCCCTGT‑AACCGCCCCAACAGAGAATGATATTGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCAATGCCGCCCAGCTGaaaaa  >  2:2153156/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                               CTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCAT‑GTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGT‑AACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATG‑ATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGG‑CGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628487‑2628764

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