Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,535,785	0	A	G	41.4%	83.5 / 15.9	58	F25S (TTC→TCC)	yddE	PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base A (6/28);  new base G (24/0);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 4.79e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGT  >  NC_000913/1535636‑1535932                                                                                                                                                       |                                                                                                                                                       cgTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGaa                                                                                                                                                       <  2:294757/151‑1 (MQ=255)    gACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGgaaa                                                                                                                                                     >  1:304530/1‑150 (MQ=255)          aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAcc                                                                                                                                              >  2:988793/1‑151 (MQ=255)           aaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGaaaa                                                                                                                                                     >  2:527402/1‑143 (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                                    <  2:56233/151‑1 (MQ=255)                                                                                                                 aaGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                    <  1:187813/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)                                                                                                                                                       |                                                                                                                                                       CGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT‑GAGGTAAAAGCATCGACGTGGTACACCTGC‑GG‑TTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGT  >  NC_000913/1535636‑1535932

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