Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,535,789	0	C	A	42.4%	99.3 / 34.9	59	V24F (GTT→TTT)	yddE	PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base C (5/29);  new base A (25/0);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 5.49e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAA  >  NC_000913/1535645‑1535936                                                                                                                                                  |                                                                                                                                                       aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAcc                                                                                                                                                  >  2:988793/1‑151 (MQ=255)  aaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCg                                                                                                                                                 <  2:1143687/151‑1 (MQ=255)       gCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGa                                                                                                                                            <  2:421231/151‑1 (MQ=255)         aCTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAAt                                                                                                                                          <  1:21870/151‑1 (MQ=255)              aaCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACAACCGCAGAAATGcc                                                                                                                                      >  2:522576/1‑150 (MQ=255)               aCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCcgc                                                                                                                                    <  1:936198/151‑1 (MQ=255)                ccGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACGCCGGCAGAATTGCcgcg                                                                                                                                   <  1:148445/151‑1 (MQ=255)                          aaaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACAACCGCAGGAttagcgccgaacaggt                                                                                                                           >  1:308218/1‑135 (MQ=255)                          aaaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTtgt                                                                                                                         <  1:145705/151‑1 (MQ=255)                            aGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCCGAAATTCCGCGCAACgggt                                                                                                                           >  2:215541/1‑145 (MQ=255)                             gTAGCGGACGCGCAACTCGCTTTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATTGGC‑GGGGAAAACATCACCGGCAGAATTGCCGCGGAACGGTtgtta                                                                                                                       <  1:229920/150‑3 (MQ=255)                               aGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAACCCGCAGAAATGCCGCGGAACGGGTGT‑ggggg                                                                                                                   >  1:604617/1‑150 (MQ=255)                               aGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGt                                                                                                                   <  2:52321/151‑1 (MQ=255)                               aGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGt                                                                                                                   <  1:542178/151‑1 (MQ=255)                                 cGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTaa                                                                                                                 <  2:958917/151‑1 (MQ=255)                                  ggATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTaaa                                                                                                                <  1:407718/151‑1 (MQ=255)                                          acaTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACCGCAGAAATTCCCCGGAACCGGTGG‑GAGGGaaaa                                                                                                               >  2:1035986/1‑144 (MQ=255)                                           caTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGCATTGCCGCGGAACGGGTGG‑GGGGGAAAAACaac                                                                                                          >  2:536331/1‑146 (MQ=255)                                            aTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACGGCAGGATTGGCGCGAAACGGGTGT‑ggggggaaaacaaaga                                                                                                        >  1:1194121/1‑137 (MQ=255)                                             tCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACCGCAGAAATGGCGCGGAAAGGGTGT‑GGGGGAAAAGCaaaga                                                                                                        >  1:162309/1‑144 (MQ=255)                                              cGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACGGCAGAAATGCCGCGGAACGGGTGT‑GGGGGAAAAAAAacgacggg                                                                                                    >  2:969837/1‑149 (MQ=255)                                                cTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCGGCCGGAATGCCGCGGAAAGGGTGG‑GGGGGAAAAGCAACGACGTgg                                                                                                   >  2:268519/1‑150 (MQ=255)                                                    cgtcgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTaca                                                                                               <  2:1098432/150‑1 (MQ=255)                                                        gCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGc                                                                                          <  1:741296/151‑1 (MQ=255)                                                         cTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAACAACCGCCGAAATGCCGCGAAACGGGTGT‑GAGGGAAAAGCATCGACGTGGTACAcc                                                                                             >  2:1190178/1‑147 (MQ=255)                                                             gcagcaGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAAATTGCGCGAAAAGGGTGT‑GGGGGGAAAACAACGACGGGGGACACCCGC‑gg                                                                                       >  2:880739/1‑148 (MQ=255)                                                               agcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACAACCGCAGAAATGCCGCGGAAAGGGTGT‑GGGGGAAAAAAAACGACGGGGGACACCCGC‑GG‑gtttt                                                                                 >  2:536450/1‑149 (MQ=255)                                                                gcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAACCCGCAGAATTGCCGCGGAAAGGGTGT‑GGGGGGAAAAAAACGACGGGGGACACCCGC‑GG‑GTTcct                                                                                >  1:838255/1‑149 (MQ=255)                                                                     aaaGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAACACCCGCAGAAATGTCGCGGAAAGGGTGG‑GGGGGAAAAACATCGACCTGGGACACCCGC‑GG‑GTTCttt                                                                               >  1:38912/1‑147 (MQ=255)                                                                       aGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGtt                                                                         <  2:805954/151‑1 (MQ=255)                                                                        gCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAATTGCCGCGGAAAGGGTGT‑GGGGGGAAAAAATCGACGAGGGACAACTAC‑GG‑GTTTATTaa                                                                             >  1:973487/1‑146 (MQ=255)                                                                         cAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTtct                                                                       <  1:235494/151‑1 (MQ=255)                                                                          aGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACGGCAGGAATGCCCCGGAAAGGGGGT‑GTGGGGAAAGGAACGACGTGGGGCACCCCC‑GG‑TTTTCTTTAAGTtct                                                                       >  1:515489/1‑150 (MQ=255)                                                                                  aaTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACt                                                              <  1:732117/151‑1 (MQ=255)                                                                                           aCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTa                                                     <  1:789990/151‑1 (MQ=255)                                                                                           aCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAAAAACGGCAGAATTGCCGCGGAACGGTTGT‑GAGGGAAAAACATCGCCGGGGGACACCTGC‑GG‑GTTTATTAACGTTCTCCCGTGACt                                                              >  2:370226/1‑142 (MQ=255)                                                                                            cTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACGGCAGAGAGTTCGCGCAACGGTTGTGGGGGTAAAAG‑AGCCACGGGGGGGACCCGG‑GG‑TTTCTTTAACTCTTTCCTGTGaactgggaaca                                                      >  2:86226/1‑139 (MQ=255)                                                                                             tCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCCCCGAATTTTCGCGGAAAGGGTGT‑GGGGGAAAAAAAACGACGTGGGAAACCCGC‑GG‑GTTCCTTTACGTTTTCCTGTGGCTTGGGAAcca                                                     >  1:220630/1‑147 (MQ=255)                                                                                               gcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCCGCAGAATTGCCGCGGAACGGGTGT‑GGGGGAAAAGAATCGAAGGGGGACACCCGC‑GG‑GTTTATTAAAGTTCTCCTTGGACTGGGGAAATa                                                     >  1:311524/1‑147 (MQ=255)                                                                                                cgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATa                                                <  1:1154067/151‑1 (MQ=255)                                                                                                cgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAAAACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGGAAAAGCATCGACGTGGTACACCTGC‑GG‑GTTTATTAAAGTTCTCCTTTGACTGGAGAACTATCa                                                  >  2:925759/1‑149 (MQ=255)                                                                                                 gggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAAATTGCGCGGAAAGGGTGT‑GGGGGAAAAACATCCACGGGGGACACCCGCGGG‑TTTATTTAAAGTT‑TCCTTTTGCTGGAGAAATATc                                                   >  2:78969/1‑147 (MQ=255)                                                                                                    cGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGG‑AGAAATTCCGCGAAACGGGTGT‑GGGGGAAAAACATCGACGTGGGACACCCGC‑GGTTTTCTTTAAC‑TTCTTCTTTGACTGGGGAAATa                                                     >  2:616731/1‑142 (MQ=255)                                                                                                    cGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCt                                            <  2:830213/151‑1 (MQ=255)                                                                                                      aTAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGc                                          <  1:584584/151‑1 (MQ=255)                                                                                                        aaGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                        <  1:187813/151‑1 (MQ=255)                                                                                                        aaGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                        <  2:56233/151‑1 (MQ=255)                                                                                                        aaGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCa                                         <  2:774022/150‑1 (MQ=255)                                                                                                             gCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTgg                                    <  2:834699/150‑1 (MQ=255)                                                                                                              cATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAAAACCCGCAGAAATGCCGCGGAAAGGGTGT‑GAGGGAAAAAAAACGACGGGGGACACCCGC‑GG‑GTTCATTAAAGTTCTCCTGTGACTGGGGAAATATCATAACCTGCCAGAgg                                    >  2:79217/1‑149 (MQ=255)                                                                                                                tttGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAACACCCCCCGAAATGCCGCCAAACGGATGT‑GGGGTTAAAAAAACCACCTGGTACACCCCC‑CG‑TTTCCTTAACCTTCTCCCCTGACTTGAGGACTATTATACCCTCCTAGAGGc                                   >  1:505154/1‑148 (MQ=255)                                                                                                                tttGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCCGAAATTCCGCGGAACGGGTGT‑GGGGGAAAAACAACCACGTGGGAAACCCGC‑GG‑GTTTATTTAAATTTTCCCGTGACTGGGGAAATATCATAAACCGCCAGTGGcc                                  >  2:951214/1‑149 (MQ=255)                                                                                                                    cgcTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGagag                            <  2:308219/151‑1 (MQ=255)                                                                                                                    cgcTTCGCTAAGATTATCGGC‑GGGGAAAAAAAAACCCGCCGAATTTCCGCGAAAAGGGTGT‑GGGGGAAAAAAAACGACGTGGGACAACCGC‑GG‑GTTTATTAAAGTTCTCCCTTGGCTGGGGGAATATCATAACCCGCAAGTGGGCgggg                              >  2:1059530/1‑147 (MQ=255)                                                                                                                       ttCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGGAGAATTGCCGCGAAAAGGGTGT‑GGGGGGAAAACATCGACGTGGGACACCTGC‑GG‑GTTTTTTTAAGTTCTCCTGTGGATGGGGAAACATCATAACCTGCCAGGGGGCgggg                              >  2:15497/1‑144 (MQ=255)                                                                                                                        tCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCg                          <  1:887318/149‑1 (MQ=255)                                                                                                                                          ggggAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACGCCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCGAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg        <  1:537733/150‑1 (MQ=255)                                                                                                                                           gggAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGt      <  2:760729/151‑1 (MQ=255)                                                                                                                                               aaaCAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaaa  <  2:304531/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                       AAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAA  >  NC_000913/1535645‑1535936

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