Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	44.0%	120.0 / 34.5	92	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (41/10);  new base C (0/40);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 1.15e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACC  >  NC_000913/3556057‑3556357                                                                                                                                                              |                                                                                                                                                       gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACt                                                                                                                                                        <  2:110324/151‑1 (MQ=255) gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACt                                                                                                                                                        <  1:217201/151‑1 (MQ=255) gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACt                                                                                                                                                        >  1:1312740/1‑151 (MQ=255)   gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTcc                                                                                                                                                      >  2:1262653/1‑151 (MQ=255)   gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTc                                                                                                                                                       <  1:711049/150‑1 (MQ=255)    ccGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTccc                                                                                                                                                     <  1:220912/151‑1 (MQ=255)       cAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGc                                                                                                                                                  >  1:1082582/1‑151 (MQ=255)         acccccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGcc                                                                                                                                                 <  2:314239/150‑1 (MQ=255)            accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGCGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAg                                                                                                                                               >  2:1226314/1‑149 (MQ=255)               tCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAgtgt                                                                                                                                            >  2:246688/1‑149 (MQ=255)                           cttttttCCCCGCCAAAAAAAAAGCGTTG‑GGTGATATTTTCTCTTTCAACTCAAAGCGAGGCGGTATCCCCC‑GGCCC‑CTGG‑CGGCGCGGAGGG‑TTAAGAATTTT‑TTGTT‑CATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACg                                                                                                                                <  1:1236951/146‑1 (MQ=255)                           cgcTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCATGc                                                                                                                               >  1:1300234/1‑150 (MQ=255)                            gcTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                             >  2:194475/1‑151 (MQ=255)                              tttcccccGACAAAAAAAAGGGTTTC‑GTTAATTTTTTCCTTTTCGACTTAAGGCGAGACGTTATTCCCC‑GGCCC‑CTGT‑TCGCGCGGCAGT‑TTATGAATTTT‑CTGTT‑CATGCGGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                             <  1:736380/143‑1 (MQ=255)                               ttttccccaagaaaaaaaaGTGTTCGGGTGATATTTCCCCTTTCGCCTTCAACCGAGCCGTTATTCCCG‑GCCCC‑CTGG‑CCGCGCGGCAGG‑TTAAAAAAATT‑TTTTT‑CAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                                            <  1:146980/137‑1 (MQ=255)                               tttCCCCGACGAAAAAAAAGCGTTG‑GGTAATTTTTTCCTTTTCTACTCAAAGCGAGAGGGTTTTCCCC‑GGCCC‑CTGG‑TCGCGGGGCAGT‑TTATGAATATT‑TTTTT‑CATGAAGGGAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAg                                                                                                                           <  1:59341/150‑1 (MQ=255)                               tttCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                                          >  2:868161/1‑151 (MQ=255)                                  cACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                                       >  2:854480/1‑151 (MQ=255)                                    ccGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                                     >  1:191761/1‑151 (MQ=255)                                    ccGACAAAAAAAAGGCGTGC‑GGTGATTTTTCCACTTCCGACTTCAAGCGAGCCGTTATTCCCC‑GCCCC‑TTGG‑TCGCGGGGCAGT‑TTATAAATATT‑CTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACa                                                                                                                      <  1:869360/150‑1 (MQ=255)                                     cgacaaAAAAAAAGGGTTG‑GTTGATTTTTTCGCTTCCGCCTAAAACAGAGAGGTTATTCCCC‑GCCCG‑CGGG‑GCGCGCGGAAGT‑TAATAAATATT‑TTGTT‑CAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                                    <  2:434108/151‑1 (MQ=255)                                         gacAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                                >  1:648531/1‑151 (MQ=255)                                          acAAAAAGGCGTTG‑GGTAATTTTTCCCCTTCCGCCTTCAAGCGAGCCGGTATCCCCC‑GCCCC‑CGGG‑TCGCGGGGCAGT‑TTATGAATTTT‑CTGTT‑CAGGAAGAGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                                   <  2:151643/147‑1 (MQ=255)                                           cAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                              >  1:17437/1‑151 (MQ=255)                                            aaaaaaGTGTCG‑GGTAATATTTTCCCTTCCGACTCAAAGGGAGCGGGTATTCCCC‑GGCCC‑CGGG‑TCGCCCGGCAGG‑TTATGAATTTT‑CTGTT‑CAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                             <  1:255222/151‑1 (MQ=255)                                            aaaaaaGCTTTC‑GGTATTTTTTTCACTTGCGGCTGAAGAGGGGACGGTATTCCCC‑GCCCC‑CGGG‑CCGCG‑GGCGGGTTTTTGATTTTT‑CTGTT‑CACGGGGAAAAAAACACCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                               <  1:539026/149‑1 (MQ=255)                                            aaaaaaGCGTTC‑GTTGTTTTTTCCCCTTTCGACTAAAAGGGAGACGGTATT‑CCC‑CGCCC‑CTGGTTCCCGCGGCAGG‑TTAGAAATATT‑CGGTT‑CAGGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                             <  2:1147321/151‑1 (MQ=255)                                            aaaaaaGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                             >  2:140756/1‑151 (MQ=255)                                            aaaaaaGCGTTC‑GGGAATATTTTCACTTTCGCCTCCAAGCGAGACGGTATTCCCC‑GCCCC‑CGGG‑TCGCGCGGCAGT‑TTATAAATATT‑CTGTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                             <  1:1064369/151‑1 (MQ=255)                                            aaaaaaGCGTCG‑GGTGATTTTTTCCCTTCCGTCTACAAGGGGGACGGTATTCCCC‑GCCCC‑CTGC‑CCGCGGGGCAGT‑TTATGAATTTT‑CTGTC‑CAGGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                             <  1:426650/151‑1 (MQ=255)                                             aaaaaGCTTTC‑GGTGATATTTTCTTTTCCGACTTCAAGCGAGCCGGTATCCCCC‑GCCCC‑CGGG‑TCGCGCGGCAGT‑TTATAAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                            <  1:271383/151‑1 (MQ=255)                                              aaaaGCGTTG‑GGTGATTTTTCCCCTTTCGCCTTAAAGCGAGACGGTATTCCCC‑GCCCC‑GTGG‑TCGCGCGGAAGT‑TTATAAATTTT‑CTGTT‑AATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                                           <  1:1494724/151‑1 (MQ=255)                                               aaaGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                                           >  1:684368/1‑150 (MQ=255)                                                aaGCGTTC‑GGTGATTTTTTCCCTTTCGACTTAAAGCGAGACGGTTTTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGT‑TTATAAATATT‑CTGTT‑AATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                                         <  1:480563/151‑1 (MQ=255)                                                     ttC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCAAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAg                                                                                                      >  2:84309/1‑149 (MQ=255)                                                         ggTGATTTTTTCATTTTCGACTTCAAGCGAGACGTTATTCCCC‑GGCCC‑CGGG‑TCGCGCGGCAGT‑TTATAAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                                    <  1:1501601/148‑1 (MQ=255)                                                          gTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                                 >  1:1278163/1‑150 (MQ=255)                                                            gATTTTTTCCCTTTCGCCTTCAAGCGAGACGGTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                                <  1:1159522/149‑1 (MQ=255)                                                             atatTTTCTCTTTCGACTTCAAGCGAGACGGTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCGGT‑TTATGAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                              <  1:1183612/150‑1 (MQ=255)                                                              ttttttCTCTTCCGCCTACAAGCGAGACGTTATTCCCC‑GCCCC‑CTGG‑TCGCGCGGCAGT‑TTATGAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                              <  1:70533/147‑1 (MQ=255)                                                              ttttttCCCTTTCGACTTAAAGGGAGACGTTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTTTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                               <  1:889041/146‑1 (MQ=255)                                                              ttttttCACTTCCGACTTCAAGCGAGACGGTTTTCCCC‑GGCCC‑CTGG‑CCGCGGGGCAGT‑TTAGGAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTa                                                                                                   <  2:698479/142‑1 (MQ=255)                                                              tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                            >  2:1380111/1‑151 (MQ=255)                                                                ttttCCTTTTCGCTTTCAAGCGGGACGTTTTCCCCC‑GGCCC‑CTGG‑CCGCGCGGCGGT‑TTATAAATTTT‑CTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                               <  2:1005951/146‑1 (MQ=255)                                                                ttttCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                          >  1:214451/1‑151 (MQ=255)                                                                 tttCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                                         >  2:964792/1‑151 (MQ=255)                                                                    cACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                                      >  1:1157525/1‑151 (MQ=255)                                                                     aCTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                                      >  1:463437/1‑150 (MQ=255)                                                                      cTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                                    >  1:867618/1‑151 (MQ=255)                                                                       tttCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                                   >  2:558435/1‑151 (MQ=255)                                                                        ttCGCCTTCAAGCGAGGCGGTATTCCCC‑GCCCC‑CGGG‑TCGCGGGGCAGG‑TTATGAATTTT‑CTTTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                                   <  1:457597/150‑1 (MQ=255)                                                                        ttCGACTTCAAGCGAGACGGTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGT‑TTATGAATTTT‑CTGTT‑AATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                                   <  1:1471712/150‑1 (MQ=255)                                                                          cGACTCCAAGCGAGACGGTTTTCCCC‑GCCCC‑CTGG‑CCGCGCGGCAGT‑TTATAAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                                 <  1:295877/150‑1 (MQ=255)                                                                            aCTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                              >  2:844646/1‑151 (MQ=255)                                                                            aCTTAAAGCGAGGCGGTATTCCCC‑GGCCC‑CTGG‑CCGCGCGGCGGT‑TTATGAATTTT‑CTGTT‑CAGGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                                    <  1:374326/145‑1 (MQ=255)                                                                             cTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                             >  1:598991/1‑151 (MQ=255)                                                                              ttaaaGCGAGCCGGTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGT‑TTATAAATTTT‑CTGTT‑CATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                              <  1:194664/146‑1 (MQ=255)                                                                              tcaaaGCGAG‑CGGTTTTTCCCGGGCCC‑CTGG‑TCGCGCGGCGGT‑TTATAAATTTT‑CTTTC‑CAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                                 <  1:798265/143‑1 (MQ=255)                                                                              taaaggCGAGACGGTATTCCCG‑GGCCC‑CTGG‑CCGCGGGCAAGT‑TTATAAATTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                              <  2:1356657/144‑1 (MQ=255)                                                                                cAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                          >  2:1424216/1‑151 (MQ=255)                                                                                 aggCGAGACTGTATTCCCC‑GCCCC‑CTGG‑TCGCGGGGCAGT‑TTATAAATTTT‑CTGTT‑CATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                          <  1:618276/148‑1 (MQ=255)                                                                                 aaGCGAGACGTTATTCCCC‑GGCCC‑CTGG‑TCGCGCGGCAGT‑TTATGAATATT‑CTTTT‑CATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                                         <  1:1000512/151‑1 (MQ=255)                                                                                 aaGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                          >  2:1258614/1‑150 (MQ=255)                                                                                     gagaCGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                                     >  1:225779/1‑151 (MQ=255)                                                                                      agaCGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                                    >  1:943573/1‑151 (MQ=255)                                                                                        gccgatttcccc‑GGCCCTCTGC‑GCGCGCGGCGGT‑TGTTAAATTTT‑GTGTT‑CATGCGGCAAAACCCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATg                                                                      <  2:531806/139‑1 (MQ=255)                                                                                          ggTATCCCCC‑GGCCC‑CGGG‑TCGCGCGGCAGT‑TTATAAATATT‑CTGTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                                <  2:744205/151‑1 (MQ=255)                                                                                             tttgcccgggCCC‑TTGG‑CCGCGCGTCGTG‑TTA‑GAATTTTCCTGTT‑CATGCAGGCAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                                 <  2:48184/140‑1 (MQ=255)                                                                                              tccccc‑GGCCC‑CGGG‑TCGCGCGGGAGT‑TTATAAATTTT‑CTTTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                                <  2:1494489/144‑1 (MQ=255)                                                                                                cccc‑GCCCC‑CGGG‑CCGCGGGGAGGT‑TTATAAATTTT‑CTTTT‑CAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                                <  2:1089043/144‑1 (MQ=255)                                                                                                 ccc‑GGCCC‑CTGG‑CCGCGCGGCAGG‑TTATAAATATT‑‑TGTTCCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                           <  2:889506/149‑1 (MQ=255)                                                                                                 ccc‑GGCCC‑CGGG‑TCGCGCGGAAGT‑TTATAAATATT‑TTGTT‑CAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                           <  2:273097/149‑1 (MQ=255)                                                                                                 ccc‑GCCCC‑CTGG‑TCGCGCGGCAGG‑TTATAAATATT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                           <  1:1415862/149‑1 (MQ=255)                                                                                                        cc‑cTGG‑CCGCGCGGCAGT‑TTATAAATATT‑CTGTT‑CATGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                                   <  1:1262671/151‑1 (MQ=255)                                                                                                                  gcgcgGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtg                                           >  2:1224937/1‑151 (MQ=255)                                                                                                                       ggcggtttATAAATATT‑CTGTT‑CATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                                       <  2:1503162/145‑1 (MQ=255)                                                                                                                       gCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                                       >  2:1482981/1‑150 (MQ=255)                                                                                                                       gCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                                       >  1:311862/1‑150 (MQ=255)                                                                                                                          gg‑TTATGAATATT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                     <  1:868177/149‑1 (MQ=255)                                                                                                                             ttATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                                 >  2:340545/1‑151 (MQ=255)                                                                                                                             ttATGAATATT‑CTGTT‑CATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                                   <  1:209819/149‑1 (MQ=255)                                                                                                                              tATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                                >  2:595155/1‑151 (MQ=255)                                                                                                                              tATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                                >  2:574375/1‑151 (MQ=255)                                                                                                                              tATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                                >  1:515550/1‑151 (MQ=255)                                                                                                                               aTGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                               >  1:609764/1‑151 (MQ=255)                                                                                                                                  aaTTTT‑CTGTT‑CATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                              <  2:463431/149‑1 (MQ=255)                                                                                                                                  aaTATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                            >  1:647125/1‑151 (MQ=255)                                                                                                                                     atT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc                         >  2:1061257/1‑151 (MQ=255)                                                                                                                                       t‑CTGTT‑CATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg                        <  1:1136124/150‑1 (MQ=255)                                                                                                                                          tGTT‑CAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag                      <  2:1275822/150‑1 (MQ=255)                                                                                                                                                  gCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc                <  1:194084/149‑1 (MQ=255)                                                                                                                                                       aaaaaaCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc                 <  1:889621/142‑1 (MQ=255)                                                                                                                                                        aaaaaCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccac        >  1:1567024/1‑151 (MQ=255)                                                                                                                                                         aaaaCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccacc       >  2:1232253/1‑151 (MQ=255)                                                                                                                                                           aaCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCcaccac     >  1:604919/1‑151 (MQ=255)                                                                                                                                                              tCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACAcc  >  2:293112/1‑151 (MQ=255)                                                                                                                                                              |                                                                                                                                                       GTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTC‑GGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCC‑GGCCC‑CTGG‑TCGCGCGGCAGG‑TTATGAATATT‑CTGTT‑CATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACC  >  NC_000913/3556057‑3556357

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: ‑