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

CGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTG‑CAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTG  >  NC_000913/1535636‑1535933                                                                                                                                                      |                                                                                                                                                        cgTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGaa                                                                                                                                                        <  1:790545/151‑1 (MQ=255)  gTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGaaa                                                                                                                                                       <  2:339643/151‑1 (MQ=255)     aCCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAaaca                                                                                                                                                    <  2:114336/151‑1 (MQ=255)      ccGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGgaaaaa                                                                                                                                                    >  2:136588/1‑148 (MQ=255)        gCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGgaaaa                                                                                                                                                     >  2:247888/1‑146 (MQ=255)          aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAGaacaa                                                                                                                                                   >  2:1046067/1‑147 (MQ=255)          aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACAcc                                                                                                                                               <  1:791340/151‑1 (MQ=255)           aaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGgaaaaa                                                                                                                                                    >  2:257686/1‑143 (MQ=255)             tCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCgg                                                                                                                                             <  2:481593/150‑1 (MQ=255)              cGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAGAACAAAAcc                                                                                                                                               >  1:1074291/1‑147 (MQ=255)               ggCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑g                                                                                                                                         <  2:399223/151‑1 (MQ=255)                     tCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCA‑GAAtt                                                                                                                                     >  1:832020/1‑149 (MQ=255)                      cAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGcc                                                                                                                                  <  2:737075/151‑1 (MQ=255)                          cGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCA‑GAAATGcc                                                                                                                                  >  2:233758/1‑147 (MQ=255)                           gTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAGACAAAAACAGCA‑GCAATGCcgc‑g                                                                                                                              >  1:861841/1‑149 (MQ=255)                            ttGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAAAAACCGCAGGAATTGcg                                                                                                                                  >  1:1075783/1‑145 (MQ=255)                               gCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCGGCA‑GAATTTCCGC‑GGAAAgg                                                                                                                        >  2:794851/1‑151 (MQ=255)                               gCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACgg                                                                                                                        <  2:922080/151‑1 (MQ=255)                                 gTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCA‑GAATTGGCGC‑GGAACgg                                                                                                                        >  1:593181/1‑149 (MQ=255)                                    aaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCC‑GAAATGCCGC‑GAAACAGTtg                                                                                                                     >  2:1049511/1‑149 (MQ=255)                                    aaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGCA‑GGATTTTCGC‑GGAACgggt                                                                                                                      >  2:939681/1‑146 (MQ=255)                                    aaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTtgtg                                                                                                                   <  2:846698/151‑1 (MQ=255)                                       tAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAACAACGGCA‑GAATTTTCGC‑GCAAAGGGTGTgggg                                                                                                                >  1:961513/1‑151 (MQ=255)                                          cGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTaa                                                                                                             <  1:19245/151‑1 (MQ=255)                                           ggATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACGGCA‑GAAATGTCGC‑GGAACGGGtgggggggaa                                                                                                             >  2:98772/1‑142 (MQ=255)                                           ggATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGCC‑GGATTGCCGCGGAAACGGtgtgggg                                                                                                                 >  2:960245/1‑141 (MQ=255)                                           ggATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTaaa                                                                                                            <  1:603045/151‑1 (MQ=255)                                             aTGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAg                                                                                                          <  1:443541/151‑1 (MQ=255)                                                 gcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATc                                                                                                      <  1:1103605/151‑1 (MQ=255)                                                   acaTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGa                                                                                                    <  2:129036/151‑1 (MQ=255)                                                          tGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAACAAAACCCGCA‑GAAATGGCGC‑GGAAAGGGTGTGAGGGAAAAACATCGACgggg                                                                                               >  1:857597/1‑149 (MQ=255)                                                           gtcgtcGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTac                                                                                            <  2:732652/151‑1 (MQ=255)                                                            tcgtcgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCGGCA‑GAAATTACGC‑GAAAAGGCTGTggggggaaaaaaacga                                                                                                    >  2:448884/1‑130 (MQ=255)                                                                   tgtgCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCA‑GAATTGCCGC‑GGAAAGGGTGTGGGGGGAAAGAAACCACGGGGGACACCCGCgg                                                                                    >  1:1017703/1‑151 (MQ=255)                                                                          cagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAACACCCGCA‑GAATTGTCCC‑GAAACGGGTGTGGGGGAAAAGCAACCACGTGGGACACCTCCGGttttct                                                                              >  1:523642/1‑147 (MQ=255)                                                                           agAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTa                                                                            <  1:79766/151‑1 (MQ=255)                                                                             aaaaGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACGGCA‑GAATTGCCGC‑GAAAAGGGTGTGAGGGAAAAAAAACGACGTGGGACACCCGCGGGTTCAtt                                                                             >  2:355698/1‑148 (MQ=255)                                                                                   aGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCA‑GAAATGCCGC‑CGAAAGGGTGTGGGGGAAAAAAAAAGAAGAGGGACACCTGCGGGTTCTTTAAAGtt                                                                       >  2:1089747/1‑148 (MQ=255)                                                                                        ccGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAACAACGGCC‑GAAATGCCGC‑GCAACGGGTGTGGGGGAAAAGAAACGACGTGGGACACCCGCGGGTTCATTAACGTTCTCCtgt                                                                >  1:533840/1‑150 (MQ=255)                                                                                        ccGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAACAACCGCA‑GAAATTCCGC‑GAAAAGGGTGTGAGGGAAAAACAACGACGGGGGACACCCGCGGGTTTATTAAAGTTCTcc                                                                   >  1:1076670/1‑147 (MQ=255)                                                                                           aaTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACGGCA‑GAAATGCCGC‑CGAAAGGGTGTGGGGGGAAAACATCGACGTGGGAAACCCGCGGGTTTATTTAAGTTCTCCTGTGACt                                                            >  1:450692/1‑151 (MQ=255)                                                                                              ggCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTgg                                                          <  1:939825/150‑1 (MQ=255)                                                                                                   aaCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCA‑GGAATGCCGC‑GGAAAGGGTGGGGGGGAAAAAAAACGACGTGGGACACCCGCGGGTTTATTTAAGTTTTCCttttgctgggg                                                        >  1:780770/1‑137 (MQ=255)                                                                                                   aaCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAAc                                                     <  1:742303/150‑1 (MQ=255)                                                                                                    aCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAACCCGCA‑GAATTGCCGC‑GGAAAGGGTGTGGGGGAAAAACAACGCCGCGGGACACCCGCGGGTTTTTTAAAATTCTCCTGTTGCTGGGGaa                                                      >  1:180177/1‑148 (MQ=255)                                                                                                    aCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGGAAAAAAAAACCGCA‑GAAATGCCGC‑GGAACGGGTGTGGGGGAAAAACAACGACGGGGGACACCCCCGGGTTTCTTAAAATTCTCCTGTGACTGGGGaa                                                      >  1:447480/1‑148 (MQ=255)                                                                                                       cgcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTAt                                                  <  1:877641/149‑1 (MQ=255)                                                                                                        gcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCAt                                               <  1:770680/151‑1 (MQ=255)                                                                                                              gATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCt                                          <  2:257703/150‑1 (MQ=255)                                                                                                              gATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTg                                         <  1:897207/151‑1 (MQ=255)                                                                                                                 aaGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTG‑Caa                                     <  2:582295/151‑1 (MQ=255)                                                                                                                         tttGCGCTTCGCTAAGATTATCGGCGGGGGGAACAAAACCCGCA‑GGAATTCCGC‑GGAACCGGTGTGAGGTGAAAAAAATTACGCGGGAAATCTGCGGTTTTGTTAAAGTTCTCCCGTGATTCGGGGAATATCATAGCCAGCCAAGTGGcgg                              >  2:335081/1‑149 (MQ=255)                                                                                                                            gcgcTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTG‑CAAGTGGCCGgag                           <  1:724770/150‑1 (MQ=255)                                                                                                                                  cGCTAAGATTATCGGCGGGGGGAAAAACAACGGGA‑GAATTGGCGC‑CGAAAGGGTGTGTGGGAAAAAAAACCAAGGGGTAAACCTTCGGGTTTATTAACAATGTCCTGCGACTCGGGAAGTATCATACACTG‑CAAGGGGTCgg                             >  2:413383/1‑142 (MQ=255)                                                                                                                                      aaGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTG‑CAAGTGGCCGGAGAGCGAAGGGCt                <  1:754874/151‑1 (MQ=255)                                                                                                                                          ttATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTG‑CAAGTGGCCGGAGAGCGAAGGGCTATc             <  2:557345/150‑1 (MQ=255)                                                                                                                                            aTCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTG‑CAAGTGGCCGGAGAGCGAAGGGCTATCCgg          <  2:716411/151‑1 (MQ=255)                                                                                                                                              cggcggGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTG‑CAAGTGGCCGGAGAGCGAAGGGCTATCCGGcc        <  1:1074162/151‑1 (MQ=255)                                                                                                                                                 cggGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTG‑CAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAgg     <  1:559130/151‑1 (MQ=255)                                                                                                                                                    ggAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTG‑CAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTg  <  1:1085988/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                        CGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCA‑GAATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTG‑CAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTG  >  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 < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑