Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,511,923	0	A	C	42.2%	108.3 / 11.2	64	V141G (GTT→GGT)	yhfZ	putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (27/10);  new base C (0/27);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 4.39e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCA  >  NC_000913/3511775‑3512064                                                                                                                                                       |                                                                                                                                              aCGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAacc                                                                                                                                             <  2:921450/151‑1 (MQ=255)        tCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGcc                                                                                                                                       >  2:557197/1‑150 (MQ=255)         ccAGCCCCCCGCGCTCCACGTTTGCGGCCTCGCCTTTCCGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAC‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAACCCACCGCCAt                                                                                                                                     <  1:957898/151‑1 (MQ=255)                cACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATac                                                                                                                              <  2:629739/151‑1 (MQ=255)                    cgcTTCGCGTTTGCGGCCTCCCCTTTACGGCAAATCAGCTGGTCCTCGCCACCGTAGGTGTGCGCCCCCAC‑CCCCAGCGCGAGGCATAAGCTTTTGTGCGTGAGATAACTTT‑CCCCCGCTAGTCGCGAAACCACCGCCATGTCATacac                                                                                                                            <  1:737023/149‑1 (MQ=255)                           cGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTa                                                                                                                     >  1:742616/1‑149 (MQ=255)                             tttGCGGCCTCGCTTTTACGGAAAATCACCTGGTGCTCGCCAACGTGGGTGTGCGGCCCCAA‑CCCCAGCGCGGGGCATAAGCCTTTTTGCGTGAGATAACTTC‑CCGCCGCTAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAg                                                                                                                   <  2:403021/149‑1 (MQ=255)                                ggggccccGCTTTTCCGCCAAATAACCTGGTGCTCCCCACCGTGGGTG‑GGGGCCCCAACCCCCGGCGCGGGGAATAAGCTTTTTTGGGTGAAAAACCTTC‑CCGCCGCTAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAg                                                                                                                <  1:951179/142‑1 (MQ=255)                                ggggccTCGCTTTTCCGGCAAACCACCTGGTCCCCCCAAACGTGGGTGGGCGCCCCAAC‑CCCCAGCGCAAGGCAAAACCTTTTTTGCGTGAAAAACCTTC‑CCGCCGCTGGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCa                                                                                                                 <  1:743407/143‑1 (MQ=255)                                 cGGACTCGCTTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTC‑CCGCCGCTAGCCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAg                                                                                                                <  2:645618/148‑1 (MQ=255)                                  ggccTCGCCTTTACGGCAAATCAACGGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CCCCAGCGCGAGCCATAAGCCTTTTTGCGTGAGATAATTTT‑CCGCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGaca                                                                                                             <  2:791967/147‑1 (MQ=255)                                   gACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACt                                                                                                           >  1:810076/1‑151 (MQ=255)                                        gCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCAcg                                                                                                      >  1:505068/1‑151 (MQ=255)                                         cttttCCGGAAAATCAACTGGTGCTCGCCACCGTGGTTGTGCGGCCCCAC‑CCCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATACCTTC‑CCGCCGTTAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCa                                                                                                        <  2:872857/146‑1 (MQ=255)                                         ccTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCAcgc                                                                                                     >  1:286099/1‑151 (MQ=255)                                             tACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAat                                                                                                 >  2:539991/1‑151 (MQ=255)                                                ggaaaaTCAACGGGTGCTCGCCAACGTGGGGGGGCGGCCCCAA‑CCCCAGCGCGGGGCATAACCTTTTTTGCGTGAGATAACTTC‑CCGCCGCTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGaa                                                                                                  <  2:156853/144‑1 (MQ=255)                                                ggaaaaCCAACGGGTGCCCCCAACCGTAGGTGGCCGCCCCCAC‑CCCCAGCGCGGGGCAAAAGCTTTTTTGCGGAAGATACCTTC‑CCGCCGCTAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAatat                                                                                               <  1:433664/147‑1 (MQ=255)                                                   aaagcacgtGGTCCTCCCCAACGTAGGTGTCGGCCCCCAC‑CCCCACCGAGAGGCTAAAGCTTTTTCGCGTAAGATATCTTC‑CCGCCGCTGGCGGCAAACCCGACGCCGTGTCAGACACGCCATTAAGCAGACACTCCACGCGAata                                                                                                <  1:568039/138‑1 (MQ=255)                                                    aaTCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCcgcg                                                                                          >  1:1144919/1‑151 (MQ=255)                                                        aaCTGGTGCTCGCCAACGTGGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTC‑CCGCCGCTAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCa                                                                                       <  2:793756/150‑1 (MQ=255)                                                            ggTGCTCGCCAACGTAGGGGTGCGGCCCCAC‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAAATAACTTC‑CCGCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGc                                                                                    <  2:907218/149‑1 (MQ=255)                                                             gTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtg                                                                                 >  2:890537/1‑151 (MQ=255)                                                                cccaccAACGTAGGTG‑GCGGCCCCAACCCCCCCCGGGGGGCATACGCTTTTTTGCGTAAAAAAACTTC‑CCGCCGCTAGCCGCAAAACCGCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtgtg                                                                               <  1:1066823/146‑1 (MQ=255)                                                                  cGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTgcgc                                                                            >  1:369055/1‑151 (MQ=255)                                                                      aaCGTAGGTGTGCGGCCCCAA‑CCCCAGCGCGGGGCATAAGCTTTTTTGCGTGAGATAACTTC‑CCGCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAg                                                                        <  2:771458/151‑1 (MQ=255)                                                                      aaCGTAGGGGGGCGCCCCAAC‑CCCCAGCGGGAGGAAAAACCTTTTTTGGGGGAGATACCTTC‑CCGCCGCTGGTCGCAAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAg                                                                        <  1:206469/151‑1 (MQ=255)                                                                          tgggTGTGCGGCCCCAA‑CCCCAGCGCGAGGCATAAGCTTTTTTGCGTGAGATAACTTC‑CCGCCGCTAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGt                                                                       <  2:1077590/146‑1 (MQ=255)                                                                          tAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTC‑CCGCCGCTAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaa                                                                     <  2:43836/150‑1 (MQ=255)                                                                            ggTGTGCGCCCCCAC‑CCCCAGCGCGGGGAAAAACCTTTTTTGGGGGAAATACTTTC‑CCGCCGCTAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaaa                                                                   <  1:421496/150‑1 (MQ=255)                                                                            ggTGTCCGCCCCCAC‑CCCCAGCGCGAGGAAAAAGCTTTTTTGCGGAAAAAAATTTC‑CCGCCGCTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTa                                                                      <  2:732329/147‑1 (MQ=255)                                                                                  cGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATg                                                            >  1:513962/1‑151 (MQ=255)                                                                                   ggCCCCAA‑CTCCAGCGCGAGGCATAAGCTTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAgg                                                                <  2:522411/146‑1 (MQ=255)                                                                                        ccacccccacggcgGGGCAAAACCTTTTTTGCGTGAAAAACCTTC‑CCGCCGCTGGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATc                                                       <  1:15796/140‑1 (MQ=255)                                                                                         aa‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTC‑CCGCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCaa                                                     <  2:807946/151‑1 (MQ=255)                                                                                           ccacctccgcgGGGCAAAGGCTTTTTTGCGTAAGATAGCTTC‑CGGCCGTTCTCCGCGAAACCCCCGCCATGTCATCCACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCaaa                                                    <  2:2939/144‑1 (MQ=255)                                                                                              ccAGCGCGAGGCATAAGCTTTTTTGCGTGAGATACCTTC‑CCGCCGCTAGTCGCGAAACCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTg                                                 <  2:852483/151‑1 (MQ=255)                                                                                                   gcgAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCtt                                            >  2:453636/1‑151 (MQ=255)                                                                                                    cgAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCttt                                           >  1:62615/1‑151 (MQ=255)                                                                                                       ggCATAAGCTTTTTTGCGTGAGATACCTTT‑CCGCCGCTAGTCGCAAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCttt                                           <  2:923621/148‑1 (MQ=255)                                                                                                          aTAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATcc                                     >  2:296350/1‑151 (MQ=255)                                                                                                           tAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCg                                    >  2:983507/1‑151 (MQ=255)                                                                                                           tAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCg                                    >  2:45303/1‑151 (MQ=255)                                                                                                            aaGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGc                                   >  1:1009279/1‑151 (MQ=255)                                                                                                               cttttttGCGGGAGAAATCTCC‑CCGCCTCGGGCCGCAACCCCCCCCCCATGTCATCCACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTc                                 <  1:629286/148‑1 (MQ=255)                                                                                                                 tttttGCGGGAGATAACTTC‑CCGCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGc                               <  2:997302/150‑1 (MQ=255)                                                                                                                    ttGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCa                             >  1:1121778/1‑149 (MQ=255)                                                                                                                       cGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCt                        >  2:400514/1‑151 (MQ=255)                                                                                                                           agaTACCTTC‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTc                      <  2:983641/149‑1 (MQ=255)                                                                                                                           agaTAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGt                    >  2:289755/1‑151 (MQ=255)                                                                                                                           aaaTAACTTT‑CCGCCGCTAGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCt                        <  1:83675/145‑1 (MQ=255)                                                                                                                               accTTC‑CCGCCGCTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTa                   <  2:751870/146‑1 (MQ=255)                                                                                                                               aaCTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACa                 >  1:163091/1‑150 (MQ=255)                                                                                                                                aCTTTCCCGCCGCTAGCCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCg                     <  1:364090/146‑1 (MQ=255)                                                                                                                                 cTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAAc              >  2:766573/1‑151 (MQ=255)                                                                                                                                 cTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAAc              >  1:835101/1‑151 (MQ=255)                                                                                                                                   tt‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACgg            >  1:31985/1‑151 (MQ=255)                                                                                                                                   tt‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACgg            >  1:1085492/1‑151 (MQ=255)                                                                                                                                     cccgccgCTGGCCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACggg           <  1:153083/150‑1 (MQ=255)                                                                                                                                      ccgccgCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACggg           >  2:275185/1‑150 (MQ=255)                                                                                                                                       cgccgcTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGt          >  2:934684/1‑150 (MQ=255)                                                                                                                                        gccgcTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtat        >  2:37067/1‑151 (MQ=255)                                                                                                                                         ccgcTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtat        <  2:500248/150‑1 (MQ=255)                                                                                                                                         ccgcTAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtata       <  2:1070404/151‑1 (MQ=255)                                                                                                                                                   cgAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagggca  <  2:332727/146‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                              ACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAGCCTTTTTGCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCA  >  NC_000913/3511775‑3512064

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