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

CGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTG  >  NC_000913/1535636‑1535933                                                                                                                                                       |                                                                                                                                                      cgTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGa                                                                                                                                                       <  1:830642/150‑1 (MQ=255)       cGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAaacaac                                                                                                                                                <  1:482945/151‑1 (MQ=255)         cAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGaaaaaaaa                                                                                                                                               >  1:568168/1‑148 (MQ=255)          aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACAcc                                                                                                                                             <  1:198735/151‑1 (MQ=255)              cGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAaaaacaaccgc                                                                                                                                          >  1:415214/1‑145 (MQ=255)              cGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCa                                                                                                                                         <  2:396018/151‑1 (MQ=255)                  aCTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCggacgat                                                                                                                                      >  2:842116/1‑146 (MQ=255)                    ttCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAAcc                                                                                                                                             >  2:585330/1‑141 (MQ=255)                       aaCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCcg                                                                                                                                <  2:704383/151‑1 (MQ=255)                        aCCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCCGAATTTCcg                                                                                                                                >  1:406648/1‑150 (MQ=255)                              ggCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACCGCAGAATTGCCGCGgaa                                                                                                                           >  2:591408/1‑149 (MQ=255)                                  tAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAACACCCGCGGCATTGTCGCGcaaaag                                                                                                                        >  2:59615/1‑145 (MQ=255)                                    aaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAAACGGGAGGAATGTCGCGgaagaggtggg                                                                                                                   >  2:405977/1‑143 (MQ=255)                                             aTGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCCGGAATGCCGCCGAAAGGGTGTggggg                                                                                                               >  2:581747/1‑145 (MQ=255)                                             aTGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGt                                                                                                               <  1:581825/146‑1 (MQ=255)                                                cgcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCAt                                                                                                       <  2:764688/151‑1 (MQ=255)                                                    caTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACGGCAGAAATGCCGCGGAAAGGGTGTGGGGGaaaa                                                                                                           >  1:728761/1‑143 (MQ=255)                                                    caTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGAc                                                                                                   <  1:48837/151‑1 (MQ=255)                                                     aTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACg                                                                                                  <  1:28771/151‑1 (MQ=255)                                                       cGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAAAACCCGCAGAAATTCCGCGCAACGGGTGTGGGGGAAAAGCAacgacggg                                                                                                >  2:218501/1‑149 (MQ=255)                                                       cGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACGGCAGAATTGCCGCGGAACGGGTGTGGGGGAAAAAAATCGACggg                                                                                                >  2:554437/1‑149 (MQ=255)                                                          tGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAAACGGCCGAATTGCCGCGCAAAGGTTGTGAGGGAAAAGCAACGACGTggg                                                                                              >  2:608011/1‑149 (MQ=255)                                                           gtcgtcGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCCGAATTGCCGCGGAAAGGGTGTGGGGGAAAAACAACGACGCgg                                                                                               >  1:653468/1‑148 (MQ=255)                                                              gtcgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAATTGCCGCGAAACGGGTGTGAGGGAAAAACAACCACGGGGGACAcc                                                                                         >  1:172471/1‑151 (MQ=255)                                                                cgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTg                                                                                       <  2:675898/151‑1 (MQ=255)                                                                 gCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAATTGCCGCGGAACGGGTGTGGGGGAAAAAAAACGACGTGGGaca                                                                                           >  1:678283/1‑146 (MQ=255)                                                                  cTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAATTGCCGCGGAAAGGGTGTGGGGGAAAAACATCGACGTGGGACACCcgc‑g                                                                                    >  1:12678/1‑151 (MQ=255)                                                                    gtgCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCAGAATTGCCGCGGAACGGGTGTGGGGGAAAAAAAACGACGGGGGACACCTGC‑gg                                                                                   >  2:135112/1‑150 (MQ=255)                                                                      gcagcaGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCCGCCGAAATGCCGCGGAACGGGTGGGGGGGAAAAAAAACGACGGGGGACACCCGC‑GGgtt                                                                                >  1:552376/1‑151 (MQ=255)                                                                         gcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCCGCAGAATTTCCGCGGAAAGGGTGTGGGGGAAAAACATCGACGTGGGACACCCGC‑GGGTTTAt                                                                             >  1:225819/1‑151 (MQ=255)                                                                         gcagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCAt                                                                             <  1:358018/151‑1 (MQ=255)                                                                           agAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTa                                                                           <  1:101017/151‑1 (MQ=255)                                                                             aaaaGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCGGCAGAATTGCCGCGGAAAGGGTGTGGGGGAAAAACAACGACGGGGTACACCCGC‑GGGTTTAttta                                                                          >  2:877772/1‑148 (MQ=255)                                                                              aaaGCAGTTTCCGAATGGCCTAACTCGCGAGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACg                                                                        <  1:262092/151‑1 (MQ=255)                                                                                      ttCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAAAAACGGCAGAATTGCCGCGGAAAGGGTGTGGGGGAAAAAAATCGACGGGGGACACCCGC‑GGGTTTCTTAAAATTCTcc                                                                  >  1:893643/1‑149 (MQ=255)                                                                                        ccGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCtgtg                                                              <  2:366347/151‑1 (MQ=255)                                                                                              ggCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCCGCAGGAATTCCGCCGAACCGGTGTGGGGGAAAAACAACGACGGGGGACACCTGC‑GGGTTTATTAACGGTCTCCCGTGACTgg                                                         >  2:467429/1‑150 (MQ=255)                                                                                              ggCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGga                                                        <  2:658458/151‑1 (MQ=255)                                                                                                ccTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAACACCCGCAGAAATGCCGCGGAACGGGTGTGGGGTAAAAACATCGACGCGGGACACCCGC‑GGGTTTATTTACGTTCTCCTGTGACTgggg                                                       >  1:521025/1‑148 (MQ=255)                                                                                                ccTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGgaga                                                      <  2:905540/151‑1 (MQ=255)                                                                                                ccTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGgaga                                                      <  1:499649/151‑1 (MQ=255)                                                                                                 cTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGGTGTGGGGGAAAAGCAACGACGTGGTACACCCGC‑GGGTTTATTAACGTTCTCCTGTGACTGgag                                                       >  1:528664/1‑149 (MQ=255)                                                                                                 cTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAAAAACACCCGCAGAATTTCCGCGAAACGGTTGTGGGGGAAAAACAACGACGTGGGACACCTGC‑GGGTTCATTAAAGTTCTCCTGTGACTGgag                                                       >  1:78913/1‑149 (MQ=255)                                                                                                       cgcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATc                                                <  1:773819/150‑1 (MQ=255)                                                                                                          gggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAACAACACCGGCAGAATTGCCGCGAAAAGGGTGTGGGGGAAAAACAACGACGTGGGACACCTGCGGGTTTCTTTAA‑ATTTTCCTGTGAATTGGGGAATATCATa                                             >  1:65898/1‑150 (MQ=255)                                                                                                            gCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAACACACCCGTAGGATTTGCCCGAAACCGGTGTGGGGGAAAAAAACCAACGTGGGACACCTCA‑GGGTTCTTTTAAGTGTTTCTGTTACTGGGGAAATATCAAAgac                                          >  2:609953/1‑149 (MQ=255)                                                                                                                 aaGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                     <  2:607911/151‑1 (MQ=255)                                                                                                                  aGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACAACCGCAGAATTTGCGCGGAACGGGTGTGGGGGAAAAACAACGACCGGGTACAACCGC‑GGGTTCATTAAAGTTCTCCTGTGACTGGGGGAATAATATAGcc                                          >  2:218030/1‑145 (MQ=255)                                                                                                                      gCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTg                                  <  2:894473/149‑1 (MQ=255)                                                                                                                      gCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTg                                  <  1:578002/149‑1 (MQ=255)                                                                                                                       cATTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCGGCAGAATTGGCGCGAAAAGGGTGTGGGGGAAAAGCAACGACGGGGGACACCCGC‑GGGTTCATTAAAGTTCTCCTGTGACTGGAGAACAATCATAACCTGCAAGTGGcc                               >  2:297016/1‑151 (MQ=255)                                                                                                                        aTTTGCGCTTCGCTAAGATTATCGGC‑GGGGGAAAAAACACCCGCGGGATTGCCGCGAAACGGTTGTGGGGGAAAAAAATCGACGAGGGACACCCGC‑GGGTTCATTAAAGTTCTCCTGTTGATGGGGAACTATCATAACCTGCCAGTGGCCg                              >  2:239027/1‑151 (MQ=255)                                                                                                                            gcgcTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGGAACGGGTGTGGGGGAAAAACATTCACGTGGGACACCTGC‑GGGTTTATTAAAATTCTCCTCTGGATGGAGAACAATCATAACCCGCAAGGGGGCgggg                           >  1:858684/1‑148 (MQ=255)                                                                                                                             cgcTTCGCTAAGATTATCGGC‑GGGGGAAAAAAAACCGGCCGAATTGCCGCGGAACGGGTGTGGGGGAAAAACATCGACGGGGTACACCTGC‑GGGTTTATTAAAGTTCTCCTGTGACTGGGGAAATATCATAACCTGCCAGGGGGCgggg                           >  1:573394/1‑147 (MQ=255)                                                                                                                             cgcTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGagag                         <  2:673410/151‑1 (MQ=255)                                                                                                                                  cGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGaa                     <  1:179838/150‑1 (MQ=255)                                                                                                                                      aaGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCt                <  2:124645/151‑1 (MQ=255)                                                                                                                                                gc‑ggGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAg      <  2:441761/151‑1 (MQ=255)                                                                                                                                                 c‑ggGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg     <  1:591488/151‑1 (MQ=255)                                                                                                                                                 c‑ggGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg     <  2:396788/151‑1 (MQ=255)                                                                                                                                                    gggAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGt   <  2:830952/151‑1 (MQ=255)                                                                                                                                                     ggAAAACAACACCGGCAGAATTTCCGTTAAACGGTTGTGCGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTg  >  2:855120/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGC‑GGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGC‑GGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTG  >  NC_000913/1535636‑1535933

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