Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,040,075	0	C	G	42.6%	85.6 / 20.5	62	L200V (CTG→GTG)	yedY	membrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base C (5/30);  new base G (26/0);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAG  >  NC_000913/2039926‑2040222                                                                                                                                                      |                                                                                                                                                    ggTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGc                                                                                                                                                     <  2:209480/150‑1 (MQ=255)    cGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGaa                                                                                                                                                 <  2:448402/151‑1 (MQ=255)        ggtggtGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAata                                                                                                                                              >  2:975505/1‑150 (MQ=255)        ggtggtGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGgaaaa                                                                                                                                              >  1:99020/1‑148 (MQ=255)              gCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTAt                                                                                                                                       <  1:739011/151‑1 (MQ=255)                    gATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCtttttt                                                                                                                                     >  1:220431/1‑143 (MQ=255)                    gATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGaa                                                                                                                                 <  1:242388/151‑1 (MQ=255)                    gATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGaa                                                                                                                                 <  1:528272/151‑1 (MQ=255)                     aTTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCCTATGTCg                                                                                                                                   >  2:249511/1‑148 (MQ=255)                     aTTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAg                                                                                                                                <  1:442750/151‑1 (MQ=255)                        ggTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGAt                                                                                                                             <  1:828892/151‑1 (MQ=255)                        ggTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGAt                                                                                                                             <  1:345454/151‑1 (MQ=255)                         gTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTATGTCGAAGGAtt                                                                                                                            >  2:781823/1‑151 (MQ=255)                         gTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGAtt                                                                                                                            <  2:471592/151‑1 (MQ=255)                           tttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTg                                                                                                                           <  2:308329/150‑1 (MQ=255)                            ttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAAAATTCTTATTTCGAAGGGtt                                                                                                                            >  2:411123/1‑148 (MQ=255)                             tCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCTTATTTTGAAgg                                                                                                                               >  2:302414/1‑144 (MQ=255)                               cGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCtttttttgaa                                                                                                                                 >  1:970816/1‑132 (MQ=255)                                 cTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATTTCGAAGGATTGCg                                                                                                                         >  1:301086/1‑146 (MQ=255)                                   gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTcgac                                                                                                                  <  1:145558/151‑1 (MQ=255)                                   gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTcgac                                                                                                                  <  1:190458/151‑1 (MQ=255)                                   gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGAAAAATCCTTATTTTGGAGGGATGCGtgt                                                                                                                      >  1:701525/1‑145 (MQ=255)                                       aaaTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGAAAAATCCTTATTTTggaggaat                                                                                                                            >  1:134814/1‑135 (MQ=255)                                       aaaTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAg                                                                                                              <  1:513808/151‑1 (MQ=255)                                         aTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCa                                                                                                            <  2:942219/151‑1 (MQ=255)                                         aTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCa                                                                                                            <  2:726800/151‑1 (MQ=255)                                            gCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAAt                                                                                                          <  1:712004/150‑1 (MQ=255)                                                     tGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAATATTCTTTTTACGAAGGATTTCGTCCCCACGAAGCAATGCATcctc                                                                                                  >  1:589202/1‑147 (MQ=255)                                                       ccGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGGAAAATTCTTTTGTCGCAGGAGTGGGTCTCGCCGAAAGAATGGATCCGCTc                                                                                                >  1:809033/1‑149 (MQ=255)                                                        cGATCCTACCAGCAATGCGAAGTATGTCGGTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGACAGCAGGACGGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTc                                                                                                <  1:768284/148‑1 (MQ=255)                                                        cGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTCGAAGGGTTGCGTCTCGACGAAACAATGCATCCg                                                                                                   >  2:649856/1‑145 (MQ=255)                                                            ccTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACtga                                                                                         <  1:750888/151‑1 (MQ=255)                                                              tACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTTTGTCGAAGGGTTGCGTCTCGGCGAAGGAAAGCATCCCCTCACACtga                                                                                         >  2:117776/1‑149 (MQ=255)                                                                    cAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTa                                                                                 <  1:302351/151‑1 (MQ=255)                                                                    cAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTa                                                                                 <  2:826784/151‑1 (MQ=255)                                                                      aTGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCTTTTTTCGGAGGAATGCGTCTCGACGAAACAAAGCATCCCCTCACACCGATGACCGTAgg                                                                               >  2:773179/1‑151 (MQ=255)                                                                        gCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATAACCTTATGTCGAAGGATTGCGTCTCGACGAAGCAAAGCATCCGCCCACACTGATGACCGTAgg                                                                               >  2:293566/1‑149 (MQ=255)                                                                             gTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATTCTTTTGTCGAAGGATTGCGTCTCGACGAAGCAAAGCATCCGCTCACACTGATGACCGTAGGGGGTTAt                                                                        >  2:510150/1‑151 (MQ=255)                                                                              tATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTGGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACCCTAATGACCGTAGGTTTTTAt                                                                        <  1:181730/150‑1 (MQ=255)                                                                                    gCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAgg                                                                 <  1:87384/151‑1 (MQ=255)                                                                                        tCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGt                                                              <  1:74461/150‑1 (MQ=255)                                                                                         cGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTa                                                            <  2:134976/151‑1 (MQ=255)                                                                                               aaTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTTAAAAATCCTTTTGTTGAAGGGATGGGTCTCGACGAAAGAAAGCATCCGCTCACACTGATGGCCGTAGGGGGTTTTTGGAAGGCGTTAACGcc                                                       >  1:547581/1‑150 (MQ=255)                                                                                                aTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCTTATGTCGAAGGATTTCGTCTCGACGAAACAAAGCATCCCCTCACACAGATGACCGTAGGGGGTTATGGGAAGGGGGTACCGcc                                                       >  2:370649/1‑149 (MQ=255)                                                                                                aTTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATTCTTATTTCGAAGGGTTTCGTCTCCACGAAACAAAGCATCCCCCCACACCGATGACCGTAGGGGTTTTTGGCAAGGCGTTACCGcccc                                                     >  2:457440/1‑149 (MQ=255)                                                                                                 tttATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAATTCCTTATTTCGAAGGATTTGGTCTCGGCGGAGCAATGCATCCCCTCACCCTTATGACCGGAGGGGTTTTTTGGAAGGCGTTAACGcc                                                       >  1:64171/1‑148 (MQ=255)                                                                                                   tATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCTTTTGTCGAAGGGGTGTGTCTCGACGAAACAAAGCCTCCGCCCACACTGATGACCGGAGGGGTTTTTTGGAAGGCGGTACCCCCCCaa                                                   >  1:351465/1‑150 (MQ=255)                                                                                                     tGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTTTTTCGAAGGATTGCGGCTCGACGAAGCAATGCATCCCCCCACACTGATGACCGTAGGGGGTTTTGGCAAAGCGGTACCCCCCCaaaa                                                 >  2:967433/1‑150 (MQ=255)                                                                                                     tGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTCGAAAGGATGCGTCTCGACGGAGCAATGCATCCCCCCCCCCTGATGACCCTAGGGGTTTTTGGGAAGGGGGTACCGCCCCcaaa                                                 >  1:587000/1‑150 (MQ=255)                                                                                                       cACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGGAAAATCCTTATGTCGAAGGGTTGCGTCTCCACGGAACAAAGCAACCGCCCACACTGATGACCGTAGGGGTTTTTTGCAAGGCGGTACCCCCCCaaaa                                                 >  1:233005/1‑148 (MQ=255)                                                                                                           ggACCAGATGCCAGGCCAGCAGGACCGCTTTATCGGGGGTGGGCTGAAATATCCTTATGGCGAAGGATTGCGTCTCGCCGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTTTTTATGGCAGGGCTTTACCGCCACAAAATGgcgc                                           <  1:982772/150‑1 (MQ=255)                                                                                                             aaCAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGAAGGGCGTTACCGCCACAAAATGgcgc                                           <  1:484810/148‑1 (MQ=255)                                                                                                              aCAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATGTTGAAGGATTGCGTCTCGACGAAGCAAAGCATCCGCCCACACTGATGACCGTAGGGGTTTTTTGCAAGGCGGTACCCCCCCAAAAAGGCGCGcc                                        >  2:234581/1‑150 (MQ=255)                                                                                                                   tGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGGTTGCGTCTCGACGGAGCAAAGCATCCGCCCACACTGATGAACGTAGGGGGTTTTGGGAAGGCGGTACCCCCACAAAAAGGCGCGCCGGGGCg                                  >  1:659137/1‑151 (MQ=255)                                                                                                                            agcagGACCGCTTTATCGGCGGGGGGGGGGAATATCCTTTTTTCGCAGGATTGCGTCTCGACGGAGCAATGGATTCGCTCACACTGGTGAACCGGGGTGTTTTTGGCAAGGGGGTAACGCCACAAAAAGGCGCGGCCGTGGGGCTGGtt                           >  1:605149/1‑149 (MQ=255)                                                                                                                                ggACCGTTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCg                     <  2:840203/151‑1 (MQ=255)                                                                                                                                ggACCGCTTTATCGGGGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCg                     <  2:75561/151‑1 (MQ=255)                                                                                                                                ggACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTCGAAAGATTGCGTCTCGACGAAGCAATGCAACCGCCCACACTGATGACCGTAGGGGTTTTTGGCAAGGCGGTACCGCCACAAAAAGGCGCGCCGGGGCGGCTGATTGTGCCg                     >  2:14724/1‑151 (MQ=255)                                                                                                                                   ccGCTTTATCGGCGGTGGGTTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGAAATGCACCCGCTCACACTGATGACCGTAGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTg                   <  1:62851/150‑1 (MQ=255)                                                                                                                                              gcggTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGGGCCGTGGAAATATGGCtt       >  1:213873/1‑151 (MQ=255)                                                                                                                                              gcggTGGGCTGAAATATCCTTATGTCAAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTGGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCtt       <  1:463523/151‑1 (MQ=255)                                                                                                                                                   gggCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGCCCGTAGGTGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAg  <  2:716224/151‑1 (MQ=255)                                                                                                                                                   gggCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACGGATGACCGTGGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaaa   <  1:869296/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAG  >  NC_000913/2039926‑2040222

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