Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,711,515 T→C 5.2% M383V (ATG→GTG)  mntH_2 ← Divalent metal cation transporter MntH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,711,5150TC 5.2% 737.9 / 21.7 232M383V (ATG→GTG) mntH_2Divalent metal cation transporter MntH
Reads supporting (aligned to +/- strand):  ref base T (129/91);  new base C (2/10);  total (131/101)
Fisher's exact test for biased strand distribution p-value = 5.72e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

Only 100 of 234 total aligned reads displayed.

ATGTCTTTTAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTAAAATCCATCCAT  >  USA300TCH1516_ALE/1711425‑1711604
                                                                                          |                                                                                         
aTGTCTTTTAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACAt                                                                                           <  2:82463/91‑1 (MQ=255)
  gTCTTTTAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATcc                                                                                         <  1:1709635/91‑1 (MQ=255)
   tCTTTTAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCa                                                                                        >  1:1325375/1‑91 (MQ=255)
    cTTTTAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCaa                                                                                       >  1:1136572/1‑91 (MQ=255)
      tttAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAAtt                                                                                     <  1:204925/91‑1 (MQ=255)
       ttAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACACCCaaatt                                                                                    <  2:2599395/91‑4 (MQ=255)
       ttAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACACCCaaatt                                                                                    <  1:3536767/91‑4 (MQ=255)
       ttAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACACCCaaatt                                                                                    <  1:884988/91‑4 (MQ=255)
        tAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTg                                                                                   <  2:2465386/91‑1 (MQ=255)
         aaaaTAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTgg                                                                                  <  1:2479875/91‑1 (MQ=255)
            aTAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGgtg                                                                               <  1:284512/91‑1 (MQ=255)
            aTAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGgtg                                                                               <  2:266261/91‑1 (MQ=255)
            aTAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACACCCAAATTggatg                                                                               <  1:3515466/91‑4 (MQ=255)
            aTAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTAATGCAACGATACCAAAGATGAACATCCAATTTGGgtg                                                                               >  2:3057975/1‑91 (MQ=255)
              aaCTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGtt                                                                             >  1:665231/1‑91 (MQ=255)
                cTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTta                                                                           >  1:1690211/1‑91 (MQ=255)
                cTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTta                                                                           >  1:3071756/1‑91 (MQ=255)
                  gAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTtata                                                                         >  1:953312/1‑91 (MQ=255)
                  gAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTtata                                                                         >  1:2497638/1‑91 (MQ=255)
                   aaGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATaa                                                                        >  2:1619096/1‑91 (MQ=255)
                   aaGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATaa                                                                        >  1:1302119/1‑91 (MQ=255)
                   aaGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATaa                                                                        <  1:2594175/91‑1 (MQ=255)
                    aGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAAt                                                                       >  1:507191/1‑91 (MQ=255)
                    aGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAAt                                                                       >  1:1053205/1‑91 (MQ=255)
                     gTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATc                                                                      <  1:3847601/91‑1 (MQ=255)
                     gTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATc                                                                      <  1:3848476/91‑1 (MQ=255)
                     gTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATc                                                                      <  1:2016592/91‑1 (MQ=255)
                        cttcttTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCAcc                                                                   >  2:1111126/1‑91 (MQ=255)
                        cttcttTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCAcc                                                                   <  2:654969/91‑1 (MQ=255)
                          tcttTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACcaa                                                                 >  2:442133/1‑91 (MQ=255)
                               gAATGAGAATATACCTGTCAATATTGGAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAAtt                                                            <  2:3540019/91‑1 (MQ=255)
                                aaTGAGAATATACCTGTCAATATTGTAACAAGTGCTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTg                                                           <  1:3176992/91‑1 (MQ=255)
                                aaTGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTg                                                           <  2:3807054/91‑1 (MQ=255)
                                aaTGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTg                                                           <  1:1544692/91‑1 (MQ=255)
                                aaTGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTg                                                           <  1:1785491/91‑1 (MQ=255)
                                aaTGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTg                                                           <  1:396124/91‑1 (MQ=255)
                                 aTGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGa                                                          >  2:1235126/1‑91 (MQ=255)
                                  tGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGAt                                                         >  1:990009/1‑91 (MQ=255)
                                    agaATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGAttt                                                       <  2:3409570/91‑1 (MQ=255)
                                    agaATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGAttt                                                       <  1:624908/91‑1 (MQ=255)
                                     gaATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTc                                                      <  1:797998/91‑1 (MQ=255)
                                     gaATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTc                                                      >  2:1297033/1‑91 (MQ=255)
                                     gaATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTc                                                      <  1:1132781/91‑1 (MQ=255)
                                     gaATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTc                                                      <  1:3675254/91‑1 (MQ=255)
                                     gaATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGACTTc                                                      >  1:3000757/1‑91 (MQ=255)
                                      aaTATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCt                                                     <  2:1025522/91‑1 (MQ=255)
                                      aaTATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCt                                                     <  2:1626993/91‑1 (MQ=255)
                                      aaTATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCt                                                     <  1:2060820/91‑1 (MQ=255)
                                      aaTATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCt                                                     <  2:1690211/91‑1 (MQ=255)
                                       atatACCTGTCAATACTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTAGGGTGTTTATAATCACCAACAATTGATTTCtt                                                    >  1:3313507/1‑91 (MQ=255)
                                         ataCCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCtttt                                                  >  1:843358/1‑91 (MQ=255)
                                         ataCCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCtttt                                                  >  1:486067/1‑91 (MQ=255)
                                         ataCCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCtttt                                                  >  2:2064029/1‑91 (MQ=255)
                                          taCCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCttttt                                                 >  2:2871953/1‑91 (MQ=255)
                                          taCCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCttttt                                                 >  2:810842/1‑91 (MQ=255)
                                          taCCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCttttt                                                 >  1:279233/1‑91 (MQ=255)
                                           aCCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTa                                                >  1:26718/1‑91 (MQ=255)
                                             cTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACt                                              >  2:2815092/1‑91 (MQ=255)
                                              tGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACGTCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACtt                                             >  2:2992136/1‑91 (MQ=255)
                                                   atatTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAAt                                        <  2:3513908/91‑1 (MQ=255)
                                                     atTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTATTTTTACTTGCAATCa                                      >  2:3704063/1‑91 (MQ=255)
                                                      ttGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGTTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCaa                                     >  1:3080523/1‑91 (MQ=255)
                                                      ttGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCaa                                     >  1:2741238/1‑91 (MQ=255)
                                                      ttGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCaa                                     >  2:514178/1‑91 (MQ=255)
                                                        gTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCaaaa                                   >  1:1665926/1‑91 (MQ=255)
                                                         tAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAAt                                  >  1:3697911/1‑91 (MQ=255)
                                                         tAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAAt                                  >  2:2989146/1‑91 (MQ=255)
                                                          aaCAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAAtt                                 <  1:1031754/91‑1 (MQ=255)
                                                           aCAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTg                                >  1:3532973/1‑91 (MQ=255)
                                                            cAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGc                               <  2:1741342/91‑1 (MQ=255)
                                                             aaGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCa                              >  1:2827655/1‑91 (MQ=255)
                                                                tACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCt                           <  2:2825235/91‑1 (MQ=255)
                                                                tACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCt                           >  1:1422516/1‑91 (MQ=255)
                                                                  cTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTaa                         >  1:2990367/1‑91 (MQ=255)
                                                                  cTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTaa                         >  1:3431103/1‑91 (MQ=255)
                                                                  cTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTaa                         >  1:1064148/1‑91 (MQ=255)
                                                                  cTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGAGTTCTTTTTACTTGCAATCAAAATTGCACCTaa                         <  2:3814043/91‑1 (MQ=255)
                                                                    gCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACg                       >  2:682698/1‑91 (MQ=255)
                                                                    gCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACg                       >  1:1130658/1‑91 (MQ=255)
                                                                      aaCGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTa                     >  2:2494672/1‑91 (MQ=255)
                                                                      aaCGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTa                     >  2:2725370/1‑91 (MQ=255)
                                                                      aaCGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTa                     >  2:197514/1‑91 (MQ=255)
                                                                       aCGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTaa                    >  1:3825650/1‑91 (MQ=255)
                                                                        cGATACCAAAGATGAACATCCAGTTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAAt                   <  2:3133971/91‑1 (MQ=255)
                                                                          aTACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAg                 <  1:1235126/91‑1 (MQ=255)
                                                                           tACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAgg                <  2:2018791/91‑1 (MQ=255)
                                                                            aCCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATCGCACCTAACGTAATAGGt               >  2:3694748/1‑91 (MQ=255)
                                                                             ccAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTa              <  2:3313507/91‑1 (MQ=255)
                                                                              cAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTaa             <  2:3429568/91‑1 (MQ=255)
                                                                              cAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTaa             >  2:1084700/1‑91 (MQ=255)
                                                                              cAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTaa             >  2:1066018/1‑91 (MQ=255)
                                                                              cAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTaa             >  2:3176610/1‑91 (MQ=255)
                                                                               aaaGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTaaa            >  1:3063443/1‑91 (MQ=255)
                                                                                aaGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAAt                   <  2:1337086/83‑1 (MQ=255)
                                                                                aaGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAAt                   >  1:1337086/1‑83 (MQ=255)
                                                                                  gATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTAAAatc         >  2:2313653/1‑91 (MQ=255)
                                                                                    tGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTAAAatcca       <  2:2110680/91‑1 (MQ=255)
                                                                                     gAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTAAAatccat      <  1:3809809/91‑1 (MQ=255)
                                                                                      aaCACCCAAATTGGATGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTAAAatccatc     >  2:185854/1‑91 (MQ=255)
                                                                                        cATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTAAAAtccatcca   <  1:3557250/91‑1 (MQ=255)
                                                                                          |                                                                                         
ATGTCTTTTAAAATAACTGAAGTACTTCTTTGAATGAGAATATACCTGTCAATATTGTAACAAGTACTGCAACGATACCAAAGATGAACATCCAATTTGGGTGTTTATAATCACCAACAATTGATTTCTTTTTACTTGCAATCAAAATTGCACCTAACGTAATAGGTAAAATCCATCCAT  >  USA300TCH1516_ALE/1711425‑1711604

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