Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,810,873	A→G	16.2%	H35R (CAC→CGC) ‡	mprA →	transcriptional repressor of microcin B17 synthesis and multidrug efflux

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,810,873	0	A	G	16.2%	166.4 / 25.1	74	H35R (CAC→CGC)	mprA	transcriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base G (6/6);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

AGTTCGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGA  >  NC_000913/2810776‑2810972                                                                                                   |                                                                                                    aGTTCGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcaca                                                                                                   <  2:575836/100‑1 (MQ=255)  gTTCGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacat                                                                                                  >  1:813930/1‑100 (MQ=255)  gTTCGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacat                                                                                                  >  1:604088/1‑100 (MQ=255)   ttCGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatg                                                                                                 <  2:813928/100‑1 (MQ=255)     cGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgca                                                                                               <  1:2236516/100‑1 (MQ=255)     cGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgca                                                                                               <  2:802627/100‑1 (MQ=255)       tttACGCCCATTGAACAAATGCT‑AAAATTTAGCGCCGGCCGCCACGAAGATTTTCCTTCTCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCaaa                                                                                             >  2:1564455/1‑100 (MQ=255)        ttACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaag                                                                                            <  2:1888936/100‑1 (MQ=255)        ttACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaag                                                                                            >  2:410541/1‑100 (MQ=255)         tACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaaga                                                                                           <  1:558267/100‑2 (MQ=255)          aCGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagac                                                                                          <  1:2141137/100‑3 (MQ=255)           cGCCCATTGAACAAATGCTAAAAATTTCGCGCCGGCCCCCAGGAAAATTTTCCTAATCAGGAGATCCTTCTCACCCGTCTTTGCATGCACATGAAaagca                                                                                          >  2:1300657/1‑100 (MQ=255)           cGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaa                                                                                         >  1:2220730/1‑100 (MQ=255)           cGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaa                                                                                         >  1:241648/1‑100 (MQ=255)            gCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaag                                                                                        <  2:548958/100‑1 (MQ=255)             cccATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGTGCATGCAaagacgag                                                                                       <  2:1367535/100‑6 (MQ=255)             cccATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaagc                                                                                       <  1:604933/100‑1 (MQ=255)              ccATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGAATCGTCTTTGCATGTGCATGCAaagacgagt                                                                                      <  2:337604/100‑7 (MQ=38)               cATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCCTGTGCATGCAaagacgagtc                                                                                     <  2:1456252/100‑8 (MQ=255)               cATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCATTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATTCAaagcaaggt                                                                                      >  2:510999/1‑97 (MQ=255)                aTTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAgctgc                                                                                    >  2:137292/1‑100 (MQ=255)                  tGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGctgctg                                                                                  >  2:455068/1‑100 (MQ=255)                   gAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTgg                                                                                 <  1:1371242/100‑1 (MQ=255)                    aaCAAATGCT‑AACATTTCGCGCCAGCCGCCACGAAGATTGTCCTCATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGga                                                                                <  2:2169103/100‑1 (MQ=255)                    aaCAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGCCTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGga                                                                                <  2:1278606/100‑1 (MQ=255)                    aaCAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGga                                                                                <  1:1364077/100‑1 (MQ=255)                     aCAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGgag                                                                               <  2:553959/100‑1 (MQ=255)                       aaaTGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGaa                                                                             <  1:795695/100‑1 (MQ=255)                        aaTGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAAc                                                                            <  1:984383/100‑1 (MQ=255)                        aaTGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAAc                                                                            <  2:1512211/100‑1 (MQ=255)                           gCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCg                                                                          >  2:852820/1‑99 (MQ=255)                           gCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGc                                                                         <  2:1956256/100‑1 (MQ=255)                           gCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGc                                                                         >  1:1591874/1‑100 (MQ=255)                            cT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCa                                                                        <  2:1848739/100‑1 (MQ=255)                            cT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCa                                                                        <  2:1591862/100‑1 (MQ=255)                            cT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCa                                                                        >  1:829863/1‑100 (MQ=255)                             t‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCaa                                                                       >  1:1870409/1‑100 (MQ=255)                               aaaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAAt                                                                      >  2:1547796/1‑100 (MQ=255)                                aaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAAt                                                                      <  1:1901954/99‑1 (MQ=255)                                aaaTTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATa                                                                     >  1:1559519/1‑100 (MQ=255)                                    ttCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAAt                                                                 <  2:190610/100‑1 (MQ=255)                                      cgcgCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGc                                                               <  2:2152897/100‑1 (MQ=255)                                      cgcgCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGc                                                               >  1:229151/1‑100 (MQ=255)                                         gcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCGTCTGGCTCGTCTTTGCGTGCACATGGAAAGCAAGCTGCTGGAGATCCGCAATAAAATGc                                                               <  2:115617/97‑1 (MQ=255)                                             gccgccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGc                                                        >  1:701788/1‑100 (MQ=255)                                                gccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCa                                                     <  1:136936/100‑1 (MQ=255)                                                gccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCa                                                     >  1:33632/1‑100 (MQ=255)                                                    cGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAgggg                                                 <  1:1438044/100‑1 (MQ=255)                                                    cGAAGATTTTCCTTATCAGGAGATCCTTCTGACTAGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAAGCGCAATAAAATGCTGAAGGCTCAgggg                                                 >  2:2251758/1‑100 (MQ=255)                                                      aaGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGAt                                               <  2:2286491/100‑1 (MQ=255)                                                       aGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGAtt                                              <  2:986614/100‑1 (MQ=255)                                                           tttCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACg                                          <  1:1415739/100‑1 (MQ=255)                                                           tttCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACg                                          <  2:2253140/100‑1 (MQ=255)                                                               cTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGAc                                      >  1:148734/1‑100 (MQ=255)                                                               cTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGAc                                      >  2:2323431/1‑100 (MQ=255)                                                                  aTCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGTACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgtt                                   >  2:1529123/1‑100 (MQ=255)                                                                  aTCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgtt                                   <  1:1627927/100‑1 (MQ=255)                                                                   tCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                  >  2:1594073/1‑100 (MQ=255)                                                                   tCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttg                                  <  2:1707593/100‑1 (MQ=255)                                                                    cAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGCTTAACGAGACgttgt                                 >  2:1541628/1‑100 (MQ=255)                                                                    cAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgttgt                                 >  1:1164954/1‑100 (MQ=255)                                                                    cAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACgatgt                                 >  1:1509074/1‑100 (MQ=255)                                                                         gaTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATg                            >  2:1913210/1‑100 (MQ=255)                                                                          aTCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATgg                           <  2:1608269/100‑1 (MQ=255)                                                                               tCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTg                      <  1:1064622/100‑1 (MQ=255)                                                                                 tGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGAt                    >  1:1557189/1‑100 (MQ=255)                                                                                   agcttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGGGAACCGCAATAAGATGCTGCAGTCTCAGGCGATTAACGAGATTGTGTTTATGGCGTTGATTa                  >  2:1671484/7‑100 (MQ=25)                                                                                    gcttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTAc                 >  1:1503092/6‑100 (MQ=255)                                                                                    cTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTAc                 >  1:2356907/1‑100 (MQ=255)                                                                                     cttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACg                >  1:2109621/5‑100 (MQ=255)                                                                                      ttgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGc               >  2:2035104/4‑100 (MQ=255)                                                                                      cGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGc               >  2:1519289/1‑100 (MQ=255)                                                                                       tgcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCt              >  1:687178/3‑100 (MQ=255)                                                                                        gcTTTGCATGTGCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTg             >  2:617468/2‑100 (MQ=255)                                                                                         cTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTgg            <  2:432566/100‑1 (MQ=255)                                                                                          tttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa           >  2:1349986/1‑100 (MQ=255)                                                                                          tttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa           >  1:956517/1‑100 (MQ=255)                                                                                          tttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa           >  2:1790279/1‑100 (MQ=255)                                                                                           ttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAg          <  1:570041/100‑1 (MQ=255)                                                                                           ttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAg          <  1:1487688/100‑1 (MQ=255)                                                                                            tgcatgcaCATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGt         >  2:277638/1‑100 (MQ=255)                                                                                              catgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtct       >  1:2087291/1‑100 (MQ=255)                                                                                               atgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGtctc      <  1:527359/100‑1 (MQ=255)                                                                                                tgcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCa     >  1:1818875/1‑100 (MQ=255)                                                                                                 gcacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAg    <  1:1751281/100‑1 (MQ=255)                                                                                                  cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg   <  1:1994531/100‑1 (MQ=255)                                                                                                  cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAgg   >  1:134555/1‑100 (MQ=255)                                                                                                  cacatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGACTACGCTGGAGCCTCAgg   >  2:917694/1‑100 (MQ=255)                                                                                                   acatgcaAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGACTCAGGa  >  1:147632/1‑100 (MQ=255)                                                                                                   |                                                                                                    AGTTCGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGGA  >  NC_000913/2810776‑2810972

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: ‑