Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,649,308	0	G	C	45.5%	128.9 / 32.2	78	intergenic (+267/‑301)	ydfC/dicB	uncharacterized protein, Qin prophage/Qin prophage; cell division inhibition protein
Reads supporting (aligned to +/- strand):  ref base G (6/36);  new base C (35/0);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 4.52e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGT  >  NC_000913/1649158‑1649458                                                                                                                                                         |                                                                                                                                                            ccTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCg                                                                                                                                                             <  1:1446024/151‑1 (MQ=255)    gAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGtt                                                                                                                                                           <  2:424339/150‑1 (MQ=255)     aaaaGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGa                                                                                                                                                         <  2:565805/151‑1 (MQ=255)      aaaGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAg                                                                                                                                                        <  1:1179168/151‑1 (MQ=255)           cAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTTCCCCCCTGAgggaac                                                                                                                                                   >  2:105420/1‑147 (MQ=255)               ccATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCggg                                                                                                                                               <  1:383266/151‑1 (MQ=255)               ccATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGGG‑TTTCTTTCCCCTTGCGGaagcccgg                                                                                                                                               >  2:1147792/1‑144 (MQ=255)                 aTAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGtt                                                                                                                                             <  2:959766/151‑1 (MQ=255)                   aaCGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTg                                                                                                                                            <  2:718793/150‑1 (MQ=255)                        gAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGt                                                                                                                                     <  1:610814/151‑1 (MQ=255)                        gAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTTCCCCCTTGAGGAGACCGGGGTGT‑GAggg                                                                                                                                     >  2:854303/1‑150 (MQ=255)                         aaaGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGtt                                                                                                                                    <  2:505213/151‑1 (MQ=255)                          aaGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTCCCCCCTGAGGAGACCCGGgtgttcgggtg                                                                                                                                    >  1:490096/1‑144 (MQ=255)                               cATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTT‑TTCTTTCCCCCTTTGGGCACCCGGGTGTT‑CAGGttacca                                                                                                                                >  1:1395595/1‑146 (MQ=255)                               cATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTTTCCCCTTGAGGGCACCCGGGTGT‑CAGGGTG‑A‑CCAt                                                                                                                            >  2:748540/1‑151 (MQ=255)                                 tGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTT‑TTCTTTCCCCCTTGGGGAACCCGGGTTgctgggtgagcctac                                                                                                                              >  1:773349/1‑137 (MQ=255)                                 tGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTTTCCCCTTTAGGACACCCGGGTGG‑CAGGGTG‑G‑CCATAc                                                                                                                          >  1:408354/1‑151 (MQ=255)                                 tGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTTCCCCCTTGGGGACACCCGGGTGT‑TAgggttgcccta                                                                                                                             >  2:700497/1‑141 (MQ=255)                                     gAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTCCCCCTTTAGGACACCCGGGTGT‑CAGGGTG‑A‑CCATACGCtt                                                                                                                      >  2:1496035/1‑151 (MQ=255)                                      aGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTa                                                                                                                     <  1:374274/151‑1 (MQ=255)                                       ggTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTaa                                                                                                                    <  1:918120/151‑1 (MQ=255)                                             tGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTTCCCCGTGAGGGCAACGGGGTGT‑TAGGGTG‑A‑ACATACGCCTaaaagaa                                                                                                               >  2:924884/1‑145 (MQ=255)                                                 tCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTCCCCCTTGGGGACACCCGGGTGT‑CAGGGTG‑A‑CCATACGCTTTAGTGAAAAacc                                                                                                          >  2:621899/1‑151 (MQ=255)                                                  cATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAAcccc                                                                                                         <  1:1371470/151‑1 (MQ=255)                                                   atagatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTTTCCCCCTTGGGGGACCCGGGGT‑T‑CTGGTTT‑T‑CCCTACGCTTTAGTgggaaacccc                                                                                                        >  1:598823/1‑142 (MQ=255)                                                    tagatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTTCCCCTTGTGGGCACCGGGGTGT‑CAGGTTGTC‑CCATA‑GCTTTAGTGAGAAACCCGc                                                                                                       >  1:778665/1‑151 (MQ=255)                                                    tagatagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTCCCCCTTGGGGACACCCGGGTGG‑CAGGTTG‑G‑CCCCACCCTTTAGTGGCAAaccccc                                                                                                       >  1:1168273/1‑149 (MQ=255)                                                        tagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTTTCCCCCTGAGGGCACCCGGGTGT‑CCGGGTG‑T‑CCAAAACCTTTAATGACAAACCCCCTGc                                                                                                    >  2:1070993/1‑150 (MQ=255)                                                        tagCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTAGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTCCCCCTTGAGGGCACCCGGGTGT‑TAGGTTG‑A‑CCATACGCTTAAATGACAAaccccccgcc                                                                                                   >  2:483305/1‑145 (MQ=255)                                                           ccTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTA‑TTCTTTTCCCCTT‑TGGACACACGGTTGT‑CTCGGTG‑A‑CCCAACAGTTAAGTGACAAACCCGcccgcacc                                                                                                >  1:1190711/1‑144 (MQ=255)                                                           ccTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACg                                                                                                <  2:1426906/151‑1 (MQ=255)                                                             tGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTTCCCCTTGAGGGCACCCGGGTGT‑TAGGGTG‑A‑CCATACGCTTAAATGAAAAACCCCCTTCAAAGcc                                                                                              >  2:1497602/1‑151 (MQ=255)                                                               tCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTCCCCCTTGAGGACACCCGGGTGT‑CAGGGTG‑A‑ACATACGCTTTAATGACAAACCCCCCGCCAAGcccc                                                                                            >  2:983295/1‑150 (MQ=255)                                                                cGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCct                                                                                            <  1:22138/150‑1 (MQ=255)                                                                  ttAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTT‑TTCTTTCCCCG‑TGAGGGCACCCGGGTGT‑CAGGGTT‑G‑CCCTACGCCTTAGTGACAAACCCCCCGCCACGCCccctt                                                                                         >  2:1284963/1‑150 (MQ=255)                                                                    aaaTTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑tt                                                                                      <  1:870156/151‑1 (MQ=255)                                                                     aaTTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCTAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTa                                                                                     <  1:1025289/151‑1 (MQ=255)                                                                          aattcgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTCCCCCCTTAGGAAACCCAGGTGT‑TAGGGTG‑A‑CCATACGCCTAAATGACAACCCCCCCGCAACGCCCTCTT‑TTTTCAAt                                                                                >  2:1465888/4‑151 (MQ=255)                                                                             tcgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTCCCCCCTGGGGGTAACCGGTGGT‑GAGGGTG‑A‑CCATACCCTTAACGGAAAAACCCCCCGCGACACCCCCTC‑CTATTAAtttt                                                                             >  2:131723/1‑151 (MQ=255)                                                                              cgACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTc                                                                            <  1:627402/151‑1 (MQ=255)                                                                                aCCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTg                                                                          <  2:1264814/151‑1 (MQ=255)                                                                                  cGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTTTCCCCCTGGGGGCACCCGGGTGT‑CAGGTTG‑A‑CCATACCCTTTAGTGACAAACCCCCCGCAACGCCCCCTG‑TTTTTAAtttttttgt                                                                        >  2:1116045/1‑146 (MQ=255)                                                                                   gTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTTCCCCCCTGAGGACACCCGGGTGT‑TAGGGTG‑G‑CCATACGCCTTAATGAGAAACCCCCCGCCACCCCCCCCT‑TTTTTAATTTTTTtgtg                                                                       >  1:385688/1‑151 (MQ=255)                                                                                    tGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTTCCCCGTGAGGGCACCCGGGTGT‑TAGGGTG‑A‑CCCTACGCTTTAGTGACAAACCCCCCGCAACGCCCCCTT‑TTTTCAATTTTTTtgtgg                                                                      >  2:122776/1‑150 (MQ=255)                                                                                      cgcTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑g                                                                   <  2:1063299/151‑1 (MQ=255)                                                                                      cgcTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCTTCCCCCCTGGGGGCACCCGGGTGTCCGGGTTA‑A‑CCA‑ACCCCTTAGTGACAAACCCCCCGCAACGCCCCCTG‑TTTTCAAATTTTTTGTGGC‑g                                                                   >  2:1448694/1‑151 (MQ=255)                                                                                        cTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTTCCCCCTGAGGGCACCCGGGTGT‑TAGGGTG‑A‑CCATACGCTTAAATGACAAAACCCCCGCAACGCCCCCTC‑TTTTTAATTTTTTTGTGGC‑ggt                                                                 >  2:1083082/1‑149 (MQ=255)                                                                                           ccGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTg                                                               <  2:66766/150‑1 (MQ=255)                                                                                            cGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTT‑TTTCCTCCCCCTTGAGGGCACCCGGGTGT‑CAGGGTT‑G‑CCCTACCCTTTAGTGACAAACCCCCTGCAAAGCCCTCTG‑TTTTCCATTTTTTTGTGGC‑Gttttgg                                                             >  2:1110883/1‑148 (MQ=255)                                                                                              gTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTTCCCCCTGTGGGCACCCGGGTGT‑TAGGGTG‑A‑CCCTACCCTTTAGTGAGAAACCCCCCGCCACGCCCCCTT‑TTTTTTATTTTTTTGTGGCGGTTTgggg                                                            >  1:70115/1‑149 (MQ=255)                                                                                                    gCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTATTTCTTTCCCCCTTG‑GGGCACCCGGGTGT‑CAGGGTG‑A‑CCATACGCTTAAGTGACAACCCCCCCGCCACGCCCCCTG‑TTTTCAATTTTTTTGGGGC‑GTTTGGGGGGATc                                                       >  2:373198/1‑149 (MQ=255)                                                                                                     cAACCCGCGAAATGGCGCGGCGGTAAGTATGGCGGGGGTTT‑TTCTTTCCCCCTT‑AGGGCACCGGGGTGT‑GAGGGTG‑G‑GCATAAGCTTTAATGAGAAACCCCCCGCCACGCCCCCTG‑TTTTCTAATTTTTTGTGAC‑GGTTGggggggtccgg                                                    >  1:901534/1‑144 (MQ=255)                                                                                                     cAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTCCTTCCCCCCTGAGGGCACCGGGGTGT‑CAGGTTG‑A‑CCCTACGCCTTAGTGACAAACCCCCCGCCACGCCCCCTG‑TTATTAATTTTTTTGTGGC‑GGTTTGCGGGATCAGt                                                    >  2:433600/1‑151 (MQ=255)                                                                                                       aCCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGttt                                                  <  2:1463457/151‑1 (MQ=255)                                                                                                             gAAATGGCGCGGCGGTAAGTATGGCGGGGGTTT‑TTCTTTCCCCTTTG‑GGGCACCGGGGGGT‑CCGGGTG‑AGCCTAACCCTTAA‑TGACAACCCCCCCGCCAAGCCCCCTG‑TTTTTAATTTTTTTGTGGCGGTTTGgggggttcgtttttcttc                                             >  1:70160/1‑137 (MQ=255)                                                                                                              aaaTGGCGCGGCGGTAAGTATGGCGGGGGGTT‑TTCTTTCCCCCTT‑AGGGGACCGGGGTGT‑GAGGGTT‑T‑CCCAAAGCCTTAGTGACAAACCCGCCGCCAAGCCCCCCGTTTTTCAATTTTTTGG‑GAC‑GTTTTGGGGGGTCAGttttttttcc                                           >  1:420696/1‑145 (MQ=255)                                                                                                              aaaTGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCg                                           <  2:136291/151‑1 (MQ=255)                                                                                                                aTGGCGCGGCGGTAAGTATGGC‑GGGGGTA‑TTTCTTCCCCCTTTAGGACACCGGGGTGT‑CAGGGTG‑A‑CCATACGCCTTAGTGACAAACCCCCCGCCACGCCCCCTG‑TTTTTAATTTTTTTGTGGA‑GGTTTGCGGGATCAGtttttatccctt                                         >  2:1490395/1‑143 (MQ=255)                                                                                                                  ggCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGa                                        <  2:1257808/150‑1 (MQ=255)                                                                                                                      cggcggTAAGTATGGC‑GGGGGTA‑TTTCTTTCCCCTTGAGGACACCCGGGTGT‑TAGGGTG‑A‑CCAAACCCCTTAATGACAAACCCCCCGCCACCCCCCCTG‑TTATTAATTTTTTTGGGGC‑GTTTTGCGGGATCAGTTTTTTTCCCTGACTGct                                   >  2:503701/1‑151 (MQ=255)                                                                                                                        gcggTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGctct                                 <  1:1201262/151‑1 (MQ=255)                                                                                                                        gcggTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGctct                                 <  1:1566458/151‑1 (MQ=255)                                                                                                                           gTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTgcc                              <  2:1195522/151‑1 (MQ=255)                                                                                                                             aaGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTgccgc                            <  1:34183/151‑1 (MQ=255)                                                                                                                             aaGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTgccg                             <  2:1227873/150‑1 (MQ=255)                                                                                                                               gTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGGGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGccc                          >  1:1136092/1‑151 (MQ=255)                                                                                                                                 aTGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCt                         <  1:1405686/150‑1 (MQ=255)                                                                                                                                     c‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTa                    <  2:323732/151‑1 (MQ=255)                                                                                                                                       ggggTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTaa                   <  1:1268725/151‑1 (MQ=255)                                                                                                                                             a‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGaa             <  2:56047/151‑1 (MQ=255)                                                                                                                                                 ccttccCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAAttt          >  2:78918/1‑151 (MQ=255)                                                                                                                                                 ccttccCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAAttt          <  1:1290032/151‑1 (MQ=255)                                                                                                                                                   ttccCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTtg        >  1:851985/1‑151 (MQ=255)                                                                                                                                                    tccCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTtg        <  1:845393/150‑1 (MQ=255)                                                                                                                                                      cccGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTtgtg      >  2:484108/1‑150 (MQ=255)                                                                                                                                                      cccGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGa     <  2:1395850/151‑1 (MQ=255)                                                                                                                                                      cccGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGa     <  2:1558803/151‑1 (MQ=255)                                                                                                                                                       ccGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGAt    >  1:258698/1‑151 (MQ=255)                                                                                                                                                       ccGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGAt    <  2:490120/151‑1 (MQ=255)                                                                                                                                                         gTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGAtgt  >  1:420822/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                            CCTGAAAAGTCAGACCATAACGCGAAAGCGCATGGCGAGGTAGCTGGTTCATAGATAGCCTGTCGTTAAATTTTCGTCGACCGTGCGCTTCCGGTTGTGGCAACCCGCGAAATGGCGCGGCGGTAAGTATGGC‑GGGGTTA‑TTCCTTCCCCGTTGAGGACACCGGGTTGT‑CAGGTTG‑A‑CCATACGCTTAAGTGACAACCCCGCTGCAACGCCCTCTG‑TTATCAATTTTCTGGTGAC‑GTTTGGCGGTATCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGT  >  NC_000913/1649158‑1649458

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