Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,040,075	0	C	G	43.8%	96.1 / 38.3	90	L200V (CTG→GTG)	yedY	membrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base C (11/39);  new base G (39/0);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 2.12e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGG  >  NC_000913/2039931‑2040224                                                                                                                                                 |                                                                                                                                                        aTGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGgaa                                                                                                                                                    >  2:462603/1‑150 (MQ=255)   ggtggtGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAatat                                                                                                                                                 <  2:476302/151‑1 (MQ=255)          ccGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTAt                                                                                                                                           >  1:83258/1‑150 (MQ=255)          ccGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATACtttt                                                                                                                                           >  2:993938/1‑148 (MQ=255)          ccGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATg                                                                                                                                          <  1:439283/151‑1 (MQ=255)          ccGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATg                                                                                                                                          <  2:175182/151‑1 (MQ=255)           cGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAataaaatta                                                                                                                                            >  2:635759/1‑142 (MQ=255)           cGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAAACCtt                                                                                                                                             >  1:902369/1‑147 (MQ=255)           cGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGt                                                                                                                                         <  1:145945/151‑1 (MQ=255)                  tGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGCGGAAAATCCTTATGTCGGAgg                                                                                                                                   >  1:406640/1‑150 (MQ=255)                   ggTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATACTTATTTTGAAg                                                                                                                                    >  2:1261085/1‑148 (MQ=255)                   ggTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGAt                                                                                                                                 <  2:1212112/151‑1 (MQ=255)                    gTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGAtt                                                                                                                                <  2:597790/151‑1 (MQ=255)                    gTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGAtt                                                                                                                                <  1:648072/151‑1 (MQ=255)                      tttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAAATCTTATGTTGa                                                                                                                                      >  2:144331/1‑143 (MQ=255)                        tCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTATGTCGAAGGAGTGCg                                                                                                                             >  2:1236442/1‑150 (MQ=255)                             tGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATTCTTTTTTCGaaaggatgggt                                                                                                                            >  2:663650/1‑137 (MQ=255)                              gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTcgac                                                                                                                      <  2:20319/151‑1 (MQ=255)                              gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGGCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGTCTGAAATATCCTTATGTCGAAGGATTGCGTCTcg                                                                                                                        <  2:198693/149‑1 (MQ=255)                               cacaAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAAAATCCTTATTTCGGAGGAATGCGTCTcgac                                                                                                                      >  1:346048/1‑150 (MQ=255)                                 caAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGaa                                                                                                                   <  2:1408317/151‑1 (MQ=255)                                    aTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCa                                                                                                                <  2:1174475/151‑1 (MQ=255)                                    aTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCa                                                                                                                <  1:542640/151‑1 (MQ=255)                                          ggCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCAt                                                                                                          <  1:177832/151‑1 (MQ=255)                                            cgcTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATcc                                                                                                        <  1:1330909/151‑1 (MQ=255)                                              cTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATTCTTATTTTGGAGGGTTGGGTCTCGACGAAACaa                                                                                                               >  1:695808/1‑142 (MQ=255)                                                tGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCt                                                                                                     <  1:583522/150‑1 (MQ=255)                                                   cGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAATATCCTTTTTTTGAAGGATTTCGTCTCGACGAAGCAATGCATcccccccc                                                                                                  >  1:565463/1‑144 (MQ=255)                                                     aaCCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTTTTTCGGAGGAGTGCGTCTCGGCGAAACAAAGCATcccccc                                                                                                    >  2:1271860/1‑142 (MQ=255)                                                     aaCCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTTGAAGGAGTGCGTCTCGACGGAGCAATGGATCCGCCcaca                                                                                                 >  2:686164/1‑149 (MQ=255)                                                      aCCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATTCTTTTTTCGGAGGAATTCGTCTCGACGGAGCAAAGCATCCCCTcac                                                                                                  >  1:65120/1‑147 (MQ=255)                                                        cTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAAAACTTTTTTTGAAGGGATTCGTCTCGCCGAAGGAAAGGATCCCCTCacac                                                                                                >  2:98538/1‑147 (MQ=255)                                                             aGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTTTTTCGAAGGATTTCGTCTCGACGAAGCAAAGCATCCCCCCACACtgatgaac                                                                                        >  1:510760/1‑148 (MQ=255)                                                               cAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATAATCTTTTTTCGAAGGAGTGCGGCTCGACGGAGCAAAGGATCCGCTCACACTGATGAACgg                                                                                      >  2:748265/1‑149 (MQ=255)                                                               cAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAAACCTTATGTCGGAGGATTTCGTCTCGACGAAGCAAAGCATCCCCTCACACTGATGAcc                                                                                        >  1:104801/1‑148 (MQ=255)                                                                    cGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAAAATCCTTTTGTCGGAGGAATTCGGCTCGGCGAAACAAAACATCCGCCCACACTTATGAACGTA‑gggg                                                                                >  1:273375/1‑148 (MQ=255)                                                                     gAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTTTGTCGGAAGGTTTCGTCTCGACGCAGCAAAGGATCCGCTCCCCCTGATGACCGTA‑gggg                                                                                >  1:646482/1‑147 (MQ=255)                                                                     gAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATTCTTTTTTCGGAAGGTTGCGGCTCGGCGGAACAAAGCATCCCCCCCCACCGATGACCGTA‑gggggt                                                                              >  1:529568/1‑147 (MQ=255)                                                                     gAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATTTCGAAGGGATGCGTCTCGACGAAGAAAAGCATCCGCTCCCACTGATGAACGTA‑GGGGtt                                                                              >  2:900369/1‑151 (MQ=255)                                                                     gAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCTACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGtt                                                                              <  1:1112379/151‑1 (MQ=255)                                                                      aaGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGAGGGGTGGAAAATCCTTATGTTGAAGAATTGCGGGTCGAGGAAGCAAAGCATCCCCTCACAACTATGATCGCA‑GGGGttt                                                                             >  1:1006266/1‑151 (MQ=255)                                                                        gTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTAt                                                                           <  1:375890/151‑1 (MQ=255)                                                                              cGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATCCTTATTTCGAAGGATTTCGTCTCGACGGAGCAAAGCATACGCTCACACTTATGACCCGA‑GGGGGTTATGG‑Caa                                                                     >  1:1252917/1‑150 (MQ=255)                                                                                  ttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGAAAAATACTTATTTTGAAGGGTTGCGTCTCCGCGAAGCAAAGCATCCGCTCACACTGATGAACCGA‑GGGGGTTATGG‑CAAGGCg                                                                 >  2:1284692/1‑150 (MQ=255)                                                                                  ttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGGCg                                                                 <  1:477752/150‑1 (MQ=255)                                                                                    cGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATTCTTTTTTCGAAGGAGTGCGTCTCGACGAAGCAAAGCATCCGCCCACACTGTTGACCGGA‑GGGGGTTTTGG‑CAAgggggt                                                               >  2:604738/1‑146 (MQ=255)                                                                                    cGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGAAATATCCTTTTTTCGGAAGATTTCGGCTCGACGGAACAAAACATTCCCCCCCACCGATGACCGGT‑GGTGGTTTTGG‑GCAGGCggt                                                               >  1:332804/1‑148 (MQ=255)                                                                                     gAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGAAATATTCTTATTTCGAAGGATTGCGGCTCCACGAAGCAATGCATCCGCTCACACTGATGAACGTA‑GGTGTTTTTTG‑GAAGGCg                                                                 >  2:809363/1‑147 (MQ=255)                                                                                     gAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGTCGTTAc                                                             <  1:978046/151‑1 (MQ=255)                                                                                      aaaCAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCATATTTCGGAGGGATGTGGCTCGCCGAAGCCAAGCATCCCCTCACACTGATGACCGGC‑GGTGTTTATTG‑GAAGGCGGTAc                                                             >  2:541329/1‑150 (MQ=255)                                                                                      aaaCAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGAGGGGTGAAAAATCCTTATTTCTAAGGATTGGGTCTCGACGAAGGAAAGAATACGCCCAACCCGATGACCGTA‑GGGGTTTTTGG‑CCAgggggt                                                               >  2:1209435/1‑144 (MQ=255)                                                                                        aCAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCCTTTGTCGAAGGATTGCGGCTCGGCGCAACAATGCATCCGCCCACACTGATGACCGCA‑GGTGGTTTTTG‑CAAgggggtatccc                                                          >  2:869396/1‑142 (MQ=255)                                                                                            tttATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAATCCTTTTGTCGAAGGGTTGCGTCTCGACGAAGCAAAGCATCCCCCCACACTGATGAACGTA‑GGGGTTTTTGG‑CCAGGGGGTACCGCcac                                                       >  1:520004/1‑150 (MQ=255)                                                                                               aTGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATAAACTTATGTCGAAGGGTTTCGTCTCCACGAAGCAAAGGAACCGCTCCCCCCGATGACCGTA‑GGGGGTTTTTG‑GAAGGCGGTAACGCCACaaa                                                    >  2:1104414/1‑150 (MQ=255)                                                                                               aTGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGAAAAAACCTTTTGTCGAAGGGGTGCGTCTCGACGAAACAAAGCATCCCCCCACACTGATGACCGTA‑GGGGGTTATTG‑GAAGGCGGTACCGCCCCaa                                                     >  2:254105/1‑149 (MQ=255)                                                                                                    ccGGAACAGTTGCGAGGCCAGCGGGCCCGTTTTATCGGGGGTGGGCTAAAATATCCTTATGTGGAAGGTTTGCGTCACGACGAAGCAATGCACCCGCTCACACTGATGACCGTA‑GGTGTTTAGGG‑CAAGGCTTTACCGCCACAAAATGgc                                               <  2:273372/150‑1 (MQ=255)                                                                                                    ccGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAAAACTTATGTCGAAGGGGTGCGTCTCGACGAAACAAAGCATCCGCTCACACTGATGACCGCG‑GGGGGTTTTGG‑GAAGGCGGTACCCCCCCaaaaaagcg                                              >  2:39372/1‑146 (MQ=255)                                                                                                     cGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAAAAACCCTATGTCGGAGGAATGCGTCTCGCCGAAACAAAGCATCCGCTCACACCGATGACCGTA‑GGGGTTTTTGG‑CAAGGCGGTACCGCCACcaaaaa                                                 >  2:960115/1‑145 (MQ=255)                                                                                                     cGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTTTGTTGAAGGATTGCGTCTCGACGAAGCAAAGCATCCCCCCCCACTTATGACCGGA‑GGGGGTTTTTG‑CAAGGCGGTACCGCCCCAAAAAGgcg                                              >  1:1103460/1‑150 (MQ=255)                                                                                                     cGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGAAAAATTCTTCTGTCGAAGGGTTTCGTCCCGGCGGAGAAAAGCATCCCCTCTAACACGTGACCGCA‑GGGGGTTTTGG‑CGAGGGGCTTCCACCCCAAAAAggg                                               >  1:308171/1‑148 (MQ=255)                                                                                                     cGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTTTGTCGCAGGATTTCGGCTCGACGAAGCAATGCATCCCCTCACACTGATGACCGGA‑GGGGTTTATGG‑CAAAGCGGTACCGCCCCCAAAAGgcg                                              >  1:1021523/1‑150 (MQ=255)                                                                                                      ggAGCAGATGCCGGGCCAGGGGGACCGCTTTATGGGGGGTGGGCTAAAATATCCTTATGTAGAAGGTTTGCGTCGCAACGAAGCATTGCATCGCCTCACACTGAGGACCGTG‑GGTGTTTATGA‑AAAGGCGTTACCGCCACAAAATGgcgcg                                            <  1:457655/151‑1 (MQ=255)                                                                                                         aCAGATGCCAGGCCAGCAGGCCCGTTTTATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑AAGGGCGTTACCGCCACAAAATGGCGCGCCg                                         <  1:1236474/151‑1 (MQ=255)                                                                                                          cAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATTCCTATTTCGAAGGATTGCGTTTCGACGGAGCAAAGCATACGCCCCCAATGATGACCGTA‑GGGGGTTATGG‑CAAGGCGTTACCGCCCCAAAAAAGAGCGcc                                          >  2:473074/1‑149 (MQ=255)                                                                                                           aGATGCCAGGCCAGCAGGCCCGTTTTATCGGCGGTGGGCTAAAATATCTTTATGTCGAAGGATTGCGTCTCGACAAAGCAATGCATCCGCTCACACTGATGCCCGTA‑GGTTTTTATGG‑CAAGGCTTTACCGCCACAAAATGGCGCGCCGGt                                       <  1:993963/151‑1 (MQ=255)                                                                                                            gATGCCAGGCCAGCAGGCCCGTTTTATGGGGGGGGGGCTAAAATATCCTTATGTGGAAGGTTTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGACCGTG‑GTTGTTTAGGG‑CAGGGCTTTACCGCCACAAAATGGCGCGCCGGt                                       <  2:1302404/150‑1 (MQ=255)                                                                                                              tGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGCTGAAAAAACCCTATGTCGAAAGATTGGGTCTCGGCGAAGCCAAGCATCCGCCCACACTGATGACCGCAGGGGGTTTTTGG‑C‑AGGGGGTATCCCCCCAAAAAAGCGCGCCGGGGCg                                    >  2:1262966/1‑151 (MQ=255)                                                                                                                ccaggccagCAGGACCGCTTTATCGGCGGTGGGGTGCAATATTCCTATTTTCAAGGGTTGCCGCTCGGCGCAGACAAGCATCCGCCCCAACTGTTGAACGCG‑GGGGGTTTTGGCCAAGGGGTACCCCCCCCAAAAGGGGGCGCgggg                                       >  2:57967/1‑146 (MQ=255)                                                                                                                   ggccagCAGGACCGCTTTATCGGGGGGGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGAAATGCATCCGCTCACACTGATGACCGTG‑GGTGTTTAGGG‑AAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACt                                 <  2:364958/149‑1 (MQ=255)                                                                                                                     ccagCAGGACCGCTTTATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACTGATGCCCGTA‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGAt                              <  1:866312/150‑1 (MQ=255)                                                                                                                        gcagGACCGTTTTATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATGGCGTCTCGGCAAAGCAATGCATCCGCTCACACTGATGCCCGTA‑GGTGTTTATGG‑CAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATtgt                           <  2:134927/150‑1 (MQ=255)                                                                                                                          agGCCCGTTTTATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTG‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGc                         <  2:346043/150‑1 (MQ=255)                                                                                                                           ggccccttttATCGGCGGGGGGCTAAAATATCCTTATGTCAAAGGATTGCGTCCCGACGAAGCAATGCACCCGCTCACCTTGATGACCGTA‑GGTTTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCg                       <  1:475460/144‑1 (MQ=255)                                                                                                                              ccGTTTTATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTg                     <  1:75365/150‑1 (MQ=255)                                                                                                                                  tttATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAGGACCGTG‑GTTGTTTATGG‑CAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGa                   <  2:1252885/148‑1 (MQ=255)                                                                                                                                    tATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTAGGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAata               <  1:229252/150‑1 (MQ=255)                                                                                                                                      tCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATgg            >  2:1278053/1‑151 (MQ=255)                                                                                                                                      tCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATgg            >  2:1063189/1‑151 (MQ=255)                                                                                                                                       cggcggTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATgg            >  2:174206/1‑150 (MQ=255)                                                                                                                                       cggcggTGGGCTAAAATATCCTTATGTCGAGGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTAGGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATgg            <  2:1078213/150‑1 (MQ=255)                                                                                                                                        gggggTGGGCTAAAATATCCTTATGTCGAAGGATTTCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTTTTTAGGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATg             <  1:1206400/145‑1 (MQ=255)                                                                                                                                         gcggTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCCCACTGATGACCGTA‑GGTGTTTATGG‑CAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCtt         <  1:1247454/151‑1 (MQ=255)                                                                                                                                          cggTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCttt        >  2:1382701/1‑151 (MQ=255)                                                                                                                                           ggggggCTAAAATACCCTTATGTCGAGGGATTGCGTCCCGAGGAAGCAATGCACCCGCTCACACGGAGGACCGTA‑GGTGTTTAGGG‑CAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATgg            <  2:1397980/143‑1 (MQ=255)                                                                                                                                           ggTGGGTTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGAGGACCGTA‑GGTTTTTAGGG‑CAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTa       <  2:895795/151‑1 (MQ=255)                                                                                                                                            gTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTa       >  1:865646/1‑150 (MQ=255)                                                                                                                                              gggCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACAAAGCAATGCACCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑AAGGGGGTTCCCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTaaa     <  1:1277539/150‑1 (MQ=255)                                                                                                                                              gggCTGAAATATCCTAATGTGGAAGGTTTGCGTCTCAGGGAAGAAATGCATCCGCTCACACGGATGACGGTA‑GGTTTTTAGGG‑CAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTACCGTGGAAATATg             <  2:340377/142‑1 (MQ=255)                                                                                                                                              gggCTAAAATATCTTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCCCACGGATGACCGTG‑GGTGTTTATGG‑CAGGGGTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAg    <  2:529560/151‑1 (MQ=255)                                                                                                                                                gCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGCAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAggg  >  1:670754/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                        ATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTA‑GGTGTTTATGG‑CAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGG  >  NC_000913/2039931‑2040224

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