Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,431	0	A	C	48.3%	81.0 / 22.7	60	G34G (GGT→GGG)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (25/6);  new base C (0/29);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 1.46e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGAC  >  NC_000913/2248284‑2248578                                                                                                                                                    |                                                                                                                                                    ggCAAACGAACCGGAAACCGGCATATAAGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCt                                                                                                                                                  >  1:933052/1‑151 (MQ=255)      aCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTAcc                                                                                                                                             >  2:55934/1‑151 (MQ=255)      aCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTAcc                                                                                                                                             <  1:567312/151‑1 (MQ=255)         aaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGAt                                                                                                                                          <  2:679483/151‑1 (MQ=255)              gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAc                                                                                                                                     >  2:649184/1‑151 (MQ=255)                  ccGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGcc                                                                                                                                 <  1:248396/151‑1 (MQ=255)                  ccGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGcc                                                                                                                                 >  1:904182/1‑151 (MQ=255)                  ccGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGcc                                                                                                                                 >  2:224998/1‑151 (MQ=255)                   cGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCa                                                                                                                                <  2:864003/151‑1 (MQ=255)                               aGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAatcatc                                                                                                                    >  2:626027/1‑151 (MQ=255)                               aGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAatcatc                                                                                                                    >  2:149968/1‑151 (MQ=255)                                   agtcccccGAGACTGGCCATCAGGAAGAAAACCACCAGCCCAACCAGAATATACGAGAGCAACCCCCCCCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTACCGATGAAACCGCCAATGGCAATCATCGTc                                                                                                                 <  1:283151/144‑1 (MQ=255)                                       cACGGAGACGGGTCATAAGGAAGAAAACCATCAGGCAAACAAAAAAAAAGGAGAGCAACCCCCCCCCCGCCCCTGCCGGAAAAATCGTGGCGCGGGAGGAAAAAAAAAGCCCTGTACGGAGGGAACCGCCAATGGCAATCATCGTCAg                                                                                                               <  1:344846/148‑1 (MQ=255)                                       cACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTg                                                                                                            >  1:315272/1‑151 (MQ=255)                                           gagaCTGGTCATCAGGAATAAACCCATCGGGCCAATAAGAAAATAGGAGAGCAACCCCCCCCCCGGCCCTGCCGGAGAAATCGTGGCGCCAGAGGCAAAAAAAGGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAc                                                                                                          <  1:742903/149‑1 (MQ=255)                                            agaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgcg                                                                                                       >  1:767183/1‑151 (MQ=255)                                            agaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgcg                                                                                                       >  2:151645/1‑151 (MQ=255)                                              aCTGGTCATAAGGAAGTAAACCCTCAGGCCAATCAGAATAGCCGAGAGCAACCCCCCGCCCGGACCTGCCTGAGAAATCGTGGCGCCAGAGGCAAAAAAAAGCCCTGCACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGAcgcg                                                                                                       <  1:31736/149‑1 (MQ=255)                                                 ggTCATCAGGAAGAAACCCATCAGGCCAATCAGCATATACGAGAGCAATCCCCCGCCGGGACCTGCCTGAAAAATCGTTGCGCAAGAGGCAACAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCtt                                                                                                   <  2:1025144/150‑1 (MQ=255)                                                 gatcataaGGAAGTAAACCACCGGGCATACCCGAATATACAAGCGAAATGCCCCCCCCGGCCCTGCGTGAAAAATCTTTGCGCCAGATGCAACAAAAACCCCGGCCCGGAAGGACCCGCCAATGGCAATCATCGTCAGGTGAcgcg                                                                                                       <  2:502064/144‑1 (MQ=255)                                                         gaaaGTAAACCACAAGGCCAACAAGCATATACGAGAGCAATCCCCCGCCCGGCCCTGCCTGAGAAATCGTGGCGCCAGAGGCAAAAAAAAGACCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCa                                                                                            <  1:966889/147‑1 (MQ=255)                                                           aaGTAACCCACCGGGCAAATCAGCAGATACGAGACCACTGCCCCGCCCGGCCCTGCCTGAGAAACCGTGGCGCCGGGGGCAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGc                                                                                        <  1:250229/151‑1 (MQ=255)                                                             gTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGt                                                                                      >  2:143237/1‑151 (MQ=255)                                                               aacccaccaGGCAAAACAGAATATAGGAGAGCAATCCCCCGCCGGGCCCGGCCTGAAAAATCGTTGCGCAGGAGGCAAAAAAAAGCCCTGTCCCGTGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTa                                                                                     <  1:336986/147‑1 (MQ=255)                                                               aacccATCAGCCCAACCAGCATATACGAGAGCAATCCCCCGCCCGGCCCTCCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTACCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTa                                                                                     <  2:738397/147‑1 (MQ=255)                                                                        gcccAACCAGCATATCCGAGAGAAATGCCCCGCCCGGCCCTGCCTGAAAAATCGTGGCGCCAGAGGCAACAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgcg                                                                            <  2:577001/148‑1 (MQ=255)                                                                        gccaaaTCAGCATATACGAGAGAAATCCCCCGCCCGGACCTGCCTGAAAAATCGTTGCGCAAGAGGCAAAAAAAAGCCTTGTACCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgc                                                                             <  2:145384/145‑1 (MQ=255)                                                                              tCAGCATATACGAGAGCAGCGCCCCGCCCGGCCCTCCCAAAAAAATCGTTGCGCCAGAGCAAACAAAAACCCCGGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCg                                                                               <  2:124197/141‑1 (MQ=255)                                                                               cAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgtg                                                                    >  2:997278/1‑151 (MQ=255)                                                                               cAGCATAAACGAGAGCAACCCCCCCCCCGGCCCTGCCGGAGAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                                     <  1:792483/150‑1 (MQ=255)                                                                                  cATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtg                                                                      >  1:792238/1‑146 (MQ=255)                                                                                   atatACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGttt                                                                >  1:873062/1‑151 (MQ=255)                                                                                   atatACGAGAGCAATCCCCCGCCCGGCCCTGCCGGAGAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGttt                                                                <  2:967659/151‑1 (MQ=255)                                                                                    tataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                               >  1:28013/1‑151 (MQ=255)                                                                                     ataCGAGAGCAATGCCCCGCCCGGCCCTGCCTGAGAAATCGTTGCGCAAGAGGCAAAAAAAAGCCCTGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                               <  2:794684/150‑1 (MQ=255)                                                                                       aCGAGAGCAAGCCCCCGCCCGGCCCTGCCTGAGAAATCGTTGCGCGAGAGGCAAAAAAAAGCCCTGTCCCGGTGGACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtttt                                                               <  1:882943/148‑1 (MQ=255)                                                                                         gagaGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGttt                                                          >  1:415322/1‑151 (MQ=255)                                                                                            agCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCgg                                                       >  1:98802/1‑151 (MQ=255)                                                                                            agCAACCCCCCGCCCGGCCCTCCCTGAAAAACCGTTGCGCCAGAGGCAAAAAAAAGCCTTTTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCg                                                        <  1:140574/150‑1 (MQ=255)                                                                                              cAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGa                                                      >  2:582027/1‑150 (MQ=255)                                                                                               actcccccGCCCGGCCCTCCCGGAAAAATCGTTGCGCCAGGGGAAACAAAAAGCCCTGTACCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCg                                                        <  2:279293/143‑1 (MQ=255)                                                                                               aaTGCCCCGCCCGGACCTGCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                                     <  2:789574/150‑1 (MQ=255)                                                                                               aaTCCCCCGCCGGGACCTGCCTGAAAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                                     <  2:691711/150‑1 (MQ=255)                                                                                                aTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAc                                                   <  2:480176/151‑1 (MQ=255)                                                                                                 tGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAcc                                                  >  2:545046/1‑151 (MQ=255)                                                                                                   ccccccccGGCCCTGCCTGAAAAATCGTTGCGCAAGAGGCAACAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAc                                                   <  1:56830/148‑1 (MQ=255)                                                                                                      cgcccgGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATAc                                              >  1:336129/1‑150 (MQ=255)                                                                                                        ccgggcccggccgaagaAACCGTACCGCCAGGGGCAAAAAAAAGCCCTGTCACAATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAc                                                   <  1:806937/129‑1 (MQ=255)                                                                                                          cgGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaaaa                                          >  2:587167/1‑150 (MQ=255)                                                                                                              cctgcctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCt                                     >  2:696818/1‑151 (MQ=255)                                                                                                               ctgcctgAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCt                                     >  2:452196/1‑150 (MQ=255)                                                                                                                   cggaaaaaTCGTGGCGCAAGAGGCAACAAAAAGCCCTGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCa                                <  2:785153/146‑1 (MQ=255)                                                                                                                        aaaTCGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATc                              <  2:128987/148‑1 (MQ=255)                                                                                                                           tCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaa                         >  2:30392/1‑150 (MQ=255)                                                                                                                                 cgcCAGAGCCAAAAAAAAGCCCGGACCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAgg                  <  2:506480/151‑1 (MQ=255)                                                                                                                                   ccggagGAAAAAAAAATCCCTGTACCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgc                 <  2:933160/147‑1 (MQ=255)                                                                                                                                           aaaaaaaaGCCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGAtt           <  1:871934/145‑1 (MQ=255)                                                                                                                                              aaaaaGCCCTGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCa     <  2:172494/151‑1 (MQ=255)                                                                                                                                                aaggACCTGTCCCGTTGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGa   <  1:597251/148‑1 (MQ=255)                                                                                                                                                 aaGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGAc  >  2:658566/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGAC  >  NC_000913/2248284‑2248578

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