Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,074,744:1 +T 5.1% coding (685/1308 nt) mdtA ← multidrug transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,074,7441.T 5.1% 678.0 / 28.2 156coding (685/1308 nt)mdtAmultidrug transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base . (79/69);  new base T (6/2);  total (85/71)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

GGTGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGCCCCTGGAACT  >  NC_002947/4074595‑4074893
                                                                                                                                                      |                                                                                                                                                     
ggTGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAgcgc                                                                                                                                                        >  1:205487/1‑150 (MQ=255)
 gTGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGcgcg                                                                                                                                                      >  2:4217024/1‑150 (MQ=255)
 gTGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGcgcg                                                                                                                                                      <  1:2713599/150‑1 (MQ=255)
  tGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCgg                                                                                                                                                     >  1:4553413/1‑150 (MQ=255)
   ggATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCggt                                                                                                                                                    <  1:1425442/150‑1 (MQ=255)
    gATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCggtg                                                                                                                                                   <  2:4660007/150‑1 (MQ=255)
       tGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGgtggtg                                                                                                                                                >  1:4606808/1‑150 (MQ=255)
        gCAGCTTGCTGTCAACGCGGTCCCAGGCCTCGACCGTCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATAGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGgtgt                                                                                                                                               <  1:1021447/150‑1 (MQ=255)
         cAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTc                                                                                                                                              <  2:3711806/150‑1 (MQ=255)
          aGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCg                                                                                                                                             >  2:3465081/1‑150 (MQ=255)
              tGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTgg                                                                                                                                         <  2:945365/150‑1 (MQ=255)
                cTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTggcg                                                                                                                                       <  1:3148382/150‑1 (MQ=255)
                cTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTggcg                                                                                                                                       <  1:5212219/150‑1 (MQ=255)
                 tGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTggcgg                                                                                                                                      >  2:180189/1‑150 (MQ=255)
                     aCCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGccac                                                                                                                                  >  2:4329762/1‑150 (MQ=255)
                     aCCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGccac                                                                                                                                  >  2:796392/1‑150 (MQ=255)
                      ccGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGccacc                                                                                                                                 >  2:4763393/1‑150 (MQ=255)
                      ccGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGccacc                                                                                                                                 >  2:2210387/1‑150 (MQ=255)
                       cgcgGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGccacca                                                                                                                                >  1:4087484/1‑150 (MQ=255)
                       cgcgGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGccacca                                                                                                                                >  1:1696544/1‑150 (MQ=255)
                        gcgGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAg                                                                                                                               >  1:989881/1‑150 (MQ=255)
                         cgGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAgg                                                                                                                              >  2:4627708/1‑150 (MQ=255)
                            tCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTg                                                                                                                           <  1:1304884/150‑1 (MQ=255)
                                     tCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGccacca                                                                                                                                <  1:1442928/136‑1 (MQ=255)
                                      cGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGa                                                                                                                 <  1:5386074/150‑1 (MQ=255)
                                          cGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTg                                                                                                             <  1:2124482/150‑1 (MQ=255)
                                          cGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGATTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGAGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTg                                                                                                             >  2:4392781/1‑150 (MQ=255)
                                           ggCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACGGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTgg                                                                                                            >  2:1556400/1‑150 (MQ=255)
                                             cAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGgcg                                                                                                          <  1:1098947/150‑1 (MQ=255)
                                                gCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTgg                                                                                                            <  2:4089976/145‑1 (MQ=255)
                                                gCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTgg                                                                                                            >  1:4089976/1‑145 (MQ=255)
                                                gCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAg                                                                                                       <  1:4320832/150‑1 (MQ=255)
                                                 cTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAgg                                                                                                      >  2:2560992/1‑150 (MQ=255)
                                                     cgcCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCca                                                                                                  <  2:1696544/150‑1 (MQ=255)
                                                      gcCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCcac                                                                                                 >  1:1941313/1‑150 (MQ=255)
                                                      gcCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCcac                                                                                                 >  2:5409295/1‑150 (MQ=255)
                                                       ccACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCcaca                                                                                                >  2:3327391/1‑150 (MQ=255)
                                                         aCTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACg                                                                                              >  2:1298807/1‑150 (MQ=255)
                                                          cTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCGCGCTGgtggtg                                                                                                                                                <  1:2649438/99‑1 (MQ=255)
                                                          cTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCGCGCTGgtggtg                                                                                                                                                >  2:2649438/1‑99 (MQ=255)
                                                           tGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGAGGTGGTGACGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGc                                                                                            >  2:644316/1‑150 (MQ=255)
                                                              ggTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCAc                                                                                         <  1:3939330/150‑1 (MQ=255)
                                                              ggTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCAc                                                                                         <  2:4948287/150‑1 (MQ=255)
                                                               gTATCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACt                                                                                        >  1:3892264/1‑150 (MQ=255)
                                                                  gcgcTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGAt                                                                                     <  1:1538396/150‑1 (MQ=255)
                                                                   cgcTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCAACAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATg                                                                                    >  1:224363/1‑150 (MQ=255)
                                                                    gcTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATgg                                                                                   <  2:3040647/150‑1 (MQ=255)
                                                                    gcTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATgg                                                                                   <  2:246358/150‑1 (MQ=255)
                                                                      tCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGgc                                                                                 >  2:2601322/1‑150 (MQ=255)
                                                                        cAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGgcgc                                                                               >  1:4940501/1‑150 (MQ=255)
                                                                             gaacgCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGa                                                                          >  2:1232956/5‑150 (MQ=255)
                                                                              aGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGAt                                                                         >  2:5080882/1‑150 (MQ=255)
                                                                               gTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATc                                                                        >  2:892388/1‑150 (MQ=255)
                                                                                  cTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTgg                                                                     <  2:3121036/150‑1 (MQ=255)
                                                                                    cAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGt                                                                   <  2:5443280/150‑1 (MQ=255)
                                                                                    cAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGt                                                                   <  2:740343/150‑1 (MQ=255)
                                                                                    cAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGt                                                                   <  2:5443170/150‑1 (MQ=255)
                                                                                      gTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGa                                                                 >  1:2900855/1‑150 (MQ=255)
                                                                                       ttCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGaa                                                                >  1:4775908/1‑150 (MQ=255)
                                                                                         cGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAAt                                                              >  2:133434/1‑150 (MQ=255)
                                                                                           gTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCg                                                            >  2:3565738/1‑150 (MQ=255)
                                                                                            tttCCGGGAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGAGCGGTTGTGTCGTTGGCGGCCACCAGGTTGCACAGATCGAGCTGGCGCAGGTCCTCACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGa                                                           >  2:518882/1‑150 (MQ=255)
                                                                                            tttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTgg                                                                     >  2:2265384/1‑140 (MQ=255)
                                                                                            tttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGa                                                           <  1:2492974/150‑1 (MQ=255)
                                                                                            tttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGa                                                           <  1:2560163/150‑1 (MQ=255)
                                                                                            tttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGa                                                           >  2:1675820/1‑150 (MQ=255)
                                                                                             ttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAg                                                          >  1:1381343/1‑150 (MQ=255)
                                                                                             ttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAg                                                          <  2:1429898/150‑1 (MQ=255)
                                                                                             ttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAg                                                          <  2:1155983/150‑1 (MQ=255)
                                                                                                        tGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAggcggg                                               >  2:1450164/1‑150 (MQ=255)
                                                                                                               gACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGtt                                        >  2:2101890/1‑150 (MQ=255)
                                                                                                               gACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGtt                                        <  2:4581240/150‑1 (MQ=255)
                                                                                                               gACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGGTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGtt                                        <  2:4606808/150‑1 (MQ=255)
                                                                                                                 cGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGa                                      <  2:5207824/150‑1 (MQ=255)
                                                                                                                   ttGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGAcc                                    >  1:2804106/1‑150 (MQ=255)
                                                                                                                    tGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCt                                   <  1:4627997/150‑1 (MQ=255)
                                                                                                                           cTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTg                            >  1:2353899/1‑150 (MQ=255)
                                                                                                                            tCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGa                           >  2:3359903/1‑150 (MQ=255)
                                                                                                                             cGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGAt                          >  1:2276482/1‑150 (MQ=255)
                                                                                                                               gTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTg                        <  2:4295327/150‑1 (MQ=255)
                                                                                                                                  tGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTc                     >  1:4205048/1‑150 (MQ=255)
                                                                                                                                    ggTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCtt                   >  1:1236886/1‑150 (MQ=255)
                                                                                                                                    ggTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCtt                   >  1:5276115/1‑150 (MQ=255)
                                                                                                                                    ggTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCtt                   >  1:5275941/1‑150 (MQ=255)
                                                                                                                                    caacgacaccacCGCGCTGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCtt                   >  1:1794628/7‑150 (MQ=255)
                                                                                                                                      tGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGa                 <  2:4013606/150‑1 (MQ=255)
                                                                                                                                      tGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGa                 <  2:2353899/150‑1 (MQ=255)
                                                                                                                                      tGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGGTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGa                 <  2:3959347/150‑1 (MQ=255)
                                                                                                                                        gacaccacCGCGCTGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGAcc               >  2:174595/3‑150 (MQ=255)
                                                                                                                                          caccacCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAg             >  1:5042223/1‑150 (MQ=255)
                                                                                                                                          caccacCGCGCTGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAg             >  2:1479729/1‑150 (MQ=255)
                                                                                                                                           accaccAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGc            >  1:177050/1‑150 (MQ=255)
                                                                                                                                            ccaccaGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGcc           >  1:3213145/1‑150 (MQ=255)
                                                                                                                                             caccaGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGccc          >  2:2216034/1‑150 (MQ=255)
                                                                                                                                             caccaGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGccc          >  1:35183/1‑150 (MQ=255)
                                                                                                                                                caccgcgCTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGCCCCTg       >  1:1261759/4‑150 (MQ=255)
                                                                                                                                                caGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGCCCCTg       >  1:2859542/1‑150 (MQ=255)
                                                                                                                                                   cgcgGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGCCCCTGGaa    <  1:963585/150‑1 (MQ=255)
                                                                                                                                                   cgcgGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGCCCCTGGaa    <  1:3705699/150‑1 (MQ=255)
                                                                                                                                                     cgGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGCCCCTGGAACt  >  1:3334370/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGTGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGCCCCTGGAACT  >  NC_002947/4074595‑4074893

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