Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	639,519	T→C	19.9%	intergenic (+11/‑234)	ahpC → / → ahpF	alkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	639,519	0	T	C	19.9%	169.8 / 38.7	86	intergenic (+11/‑234)	ahpC/ahpF	alkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (36/32);  new base C (15/2);  total (52/34)
Fisher's exact test for biased strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGC  >  NC_000913/639421‑639617                                                                                                   |                                                                                                   ttCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTc                                                                                                   >  1:31653/1‑100 (MQ=255)  tctcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCt                                                                                                  >  1:1541586/1‑100 (MQ=255)  tctcACCCAGGTGAAGTTTGCCAGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCt                                                                                                  >  2:1309546/1‑100 (MQ=255)    tcACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCttt                                                                                                >  2:1950712/1‑100 (MQ=255)     cACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTc                                                                                               >  1:1616940/1‑100 (MQ=255)     cACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTc                                                                                               >  2:1442010/1‑100 (MQ=255)     cACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTc                                                                                               >  1:77332/1‑100 (MQ=255)       cccAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCAc                                                                                             >  1:733588/1‑100 (MQ=255)        ccAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACg                                                                                            >  1:1680820/1‑100 (MQ=255)         cAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGc                                                                                           >  2:1464202/1‑100 (MQ=255)          aGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGcc                                                                                          >  2:1647134/1‑100 (MQ=255)           ggTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCa                                                                                         >  1:35492/1‑100 (MQ=255)             tGAAGTTTGTCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATa                                                                                       >  1:1182990/1‑100 (MQ=255)               aaGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAgc                                                                                     <  2:1616929/100‑1 (MQ=255)                aGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAgcg                                                                                    >  1:2159785/1‑100 (MQ=255)                 gTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAgcgg                                                                                   >  1:1768342/1‑100 (MQ=255)                   ttGCCCGTCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCTGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAgcggcg                                                                                 <  2:232656/100‑1 (MQ=255)                    tGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGt                                                                                >  2:931251/1‑100 (MQ=255)                      cccGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTg                                                                              <  2:1687604/100‑1 (MQ=255)                      cccGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTg                                                                              <  1:1194909/100‑1 (MQ=255)                      cccGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTg                                                                              >  1:345939/1‑100 (MQ=255)                       ccGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTgg                                                                             <  1:1759420/100‑1 (MQ=255)                         ggCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCg                                                                           >  2:867628/1‑100 (MQ=255)                            tAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCg                                                                        <  2:733588/100‑1 (MQ=255)                             aaaTGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGc                                                                       <  1:1136380/100‑1 (MQ=255)                              aaTGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGcc                                                                      <  2:919147/100‑1 (MQ=255)                               aTGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGccc                                                                     >  2:950360/1‑100 (MQ=255)                               aTGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCGTCGccc                                                                     >  2:1551300/1‑100 (MQ=255)                               aTGGAAAGAAGGTGAAGAAACTCTGGCTCCGACTCTGGAACTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCATAGccc                                                                     >  2:841748/1‑100 (MQ=255)                                tGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCg                                                                    >  2:269582/1‑100 (MQ=255)                                tGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCg                                                                    >  1:353167/1‑100 (MQ=255)                                tGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCGTCGccc                                                                     <  2:370830/99‑1 (MQ=255)                                  gAAAGAAGGTGAAGCAACTCTGGCTCCGTCGCTGGACCTGGTTGGTAAAATCTAAATTTTCTTCGCCTTTCACGCCATAGCTGCGTTGGCGTCGCCCGCt                                                                  <  2:1981246/100‑1 (MQ=255)                                   aaaGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTc                                                                 >  2:1076939/1‑100 (MQ=255)                                    aagaagGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCa                                                                >  1:2205894/1‑100 (MQ=255)                                         ggTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCg                                                           >  1:1794422/1‑100 (MQ=255)                                             atgCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCa                                                       >  1:1513813/3‑100 (MQ=255)                                             aaGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTc                                                        <  1:322366/99‑1 (MQ=255)                                              aGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCATCGCCCGCTCACCCCGGTCac                                                      >  2:23246/1‑100 (MQ=255)                                               gCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCact                                                     <  2:4869/100‑1 (MQ=255)                                               gCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCact                                                     <  1:1727534/100‑1 (MQ=255)                                                cAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCactt                                                    >  1:674263/1‑100 (MQ=255)                                                   ctctGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttact                                                 <  1:1184133/100‑1 (MQ=255)                                                   ctctGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttact                                                 >  2:1735179/1‑100 (MQ=255)                                                    tctGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCacttactt                                                >  1:912720/1‑100 (MQ=255)                                                     ctGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTtg                                               >  2:1206148/1‑100 (MQ=255)                                                     ctGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTtg                                               >  2:2119596/1‑100 (MQ=255)                                                       ggCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTtgtg                                             >  2:1546473/1‑100 (MQ=255)                                                        gCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTgtgt                                            <  1:1531802/100‑1 (MQ=255)                                                         cTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTa                                           >  1:1094751/1‑100 (MQ=255)                                                          tCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCTCTCACCCCGGTCACTTACTTGTGTa                                           <  2:989957/99‑1 (MQ=255)                                                          tCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTaa                                          <  2:142955/100‑1 (MQ=255)                                                          tCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTa                                           <  1:803773/99‑1 (MQ=255)                                                             gTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCt                                       >  1:1881841/1‑100 (MQ=255)                                                             gTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTATTTGTGTAAGCt                                       >  1:1534489/1‑100 (MQ=255)                                                               ctctGGACCTGGTTGGTAAAATCTAAATTTCCTTCGGCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTcc                                     >  1:1583506/1‑100 (MQ=255)                                                                  tGGATCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCg                                   <  1:1206835/99‑1 (MQ=255)                                                                  tGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCgg                                  <  1:1941183/100‑1 (MQ=255)                                                                   ggACGTGGTGGGTAAAATCTAAATTTCGTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCggg                                 <  2:593596/100‑1 (MQ=255)                                                                   ggACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCggg                                 >  1:1952088/1‑100 (MQ=255)                                                                     aCCTGGTTGGTAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGATCACTTACTTGTGTAAGCTCCCGGGGa                               >  1:33929/1‑100 (MQ=255)                                                                      ccTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGAt                              >  2:725305/1‑100 (MQ=255)                                                                        tggttggtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTc                            >  1:1321896/1‑100 (MQ=255)                                                                         ggttggtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTca                           >  2:1154276/1‑100 (MQ=255)                                                                         ggttggtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTca                           <  1:2162666/100‑1 (MQ=255)                                                                         ggttggtAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTca                           >  2:629081/1‑100 (MQ=255)                                                                           ttggtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTcaca                         <  1:1980006/100‑1 (MQ=255)                                                                             ggtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAgg                       >  1:1951623/1‑100 (MQ=255)                                                                             ggtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAgg                       >  2:2186244/1‑100 (MQ=255)                                                                              gtAAAATCTAAATTTCCTTCGTCTTTCACGCCATAACGACGTTGGCGTAGCCCGCTCACCCCGGTCACTAACTTGCGTAACCTCCCGGGGGTTCACAgg                       >  2:286308/1‑99 (MQ=255)                                                                               tAAAATCTAAATTTCCTTCTTCTTTCACGCCATAGCAGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGc                      <  2:1881910/99‑1 (MQ=255)                                                                               tAAAATCTAAATTTCCTTCGCCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCt                     >  1:1672584/1‑100 (MQ=255)                                                                                   aTCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAgcc                 >  1:1352560/1‑100 (MQ=255)                                                                                   aTCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAgcc                 >  1:1351966/1‑100 (MQ=255)                                                                                    tCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAgccg                <  1:1633063/100‑1 (MQ=255)                                                                                     cTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAgccgc               >  1:994986/1‑100 (MQ=255)                                                                                      tAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAgccgcc              <  2:2065095/100‑1 (MQ=255)                                                                                          tttCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGc          <  2:297041/100‑1 (MQ=255)                                                                                          tttCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGc          <  2:595257/100‑1 (MQ=255)                                                                                            tCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGctn        <  2:1513803/100‑2 (MQ=255)                                                                                            tCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGctc        <  2:1230589/100‑1 (MQ=255)                                                                                            tCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGctc        <  2:2069587/100‑1 (MQ=255)                                                                                            tCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGctc        <  1:2186259/100‑1 (MQ=255)                                                                                             ccTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGctct       <  2:114864/100‑1 (MQ=255)                                                                                                 cGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTACTCTGAcg   <  1:1442020/100‑1 (MQ=255)                                                                                                  gTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGAcgc  >  1:2181226/1‑100 (MQ=255)                                                                                                   tCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGAcg   <  1:421787/98‑1 (MQ=255)                                                                                                   |                                                                                                   TTCTCACCCAGGTGAAGTTTGCCCGGCTAAATGGAAAGAAGGTGAAGCAACTCTGGCTCCGTCTCTGGACCTGGTTGGTAAAATCTAAATTTCCTTCGTCTTTCACGCCATAGCGGCGTTGGCGTCGCCCGCTCACCCCGGTCACTTACTTGTGTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGC  >  NC_000913/639421‑639617

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑