Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,810,872	C→T	16.2%	H35Y (CAC→TAC) ‡	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,872	0	C	T	16.2%	173.0 / 23.2	74	H35Y (CAC→TAC)	mprA	transcriptional repressor of microcin B17 synthesis and multidrug efflux
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base T (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 = 9.20e-01

GATAGTTCGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGG  >  NC_000913/2810773‑2810971                                                                                                     |                                                                                                    gATAGTTCGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCGGGAGATCCTTCTGACTCGTCTTtgcatgc                                                                                                     <  2:2002490/100‑1 (MQ=255)    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                                                                                              <  2:802627/100‑1 (MQ=255)        cGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGcatgcacatgca                                                                                              <  1:2236516/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:241648/1‑100 (MQ=255)              cGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAaagcaa                                                                                        >  1:2220730/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                                                              >  1:229151/1‑100 (MQ=255)                                         cgcgCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGc                                                              <  2:2152897/100‑1 (MQ=255)                                            gcCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCGTCTGGCTCGTCTTTGCGTGCACATGGAAAGCAAGCTGCTGGAGATCCGCAATAAAATGc                                                              <  2:115617/97‑1 (MQ=255)                                                gccgccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGc                                                       >  1:701788/1‑100 (MQ=255)                                                   gccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCa                                                    >  1:33632/1‑100 (MQ=255)                                                   gccACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCa                                                    <  1:136936/100‑1 (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          >  1:956517/1‑100 (MQ=255)                                                                                             tttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa          >  2:1790279/1‑100 (MQ=255)                                                                                             tttGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGa          >  2:1349986/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)                                                                                                     |                                                                                                    GATAGTTCGTTTACGCCCATTGAACAAATGCT‑AAAATTTCGCGCCAGCCGCCACGAAGATTTTCCTTATCAGGAGATCCTTCTGACTCGTCTTTGCATGCACATGCAAAGCAAGCTGCTGGAGAACCGCAATAAAATGCTGAAGGCTCAGGGGATTAACGAGACGTTGTTTATGGCGTTGATTACGCTGGAGTCTCAGG  >  NC_000913/2810773‑2810971

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