Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 722,992 A→G 5.9% I44V (ATA→GTA)  feuB → Iron‑uptake system permease protein FeuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE722,9920AG 5.9% 1582.0 / 58.2 458I44V (ATA→GTA) feuBIron‑uptake system permease protein FeuB
Reads supporting (aligned to +/- strand):  ref base A (219/211);  new base G (17/10);  total (237/221)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 474 total aligned reads displayed.

TCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTAATTG  >  USA300TCH1516_ALE/722902‑723082
                                                                                          |                                                                                          
  cTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATtata                                                                                          <  2:824303/91‑1 (MQ=255)
  cTATATTATAGGTTTAAGTATGATAATACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATtata                                                                                          <  2:6491866/91‑1 (MQ=255)
    atatTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATaga                                                                                        >  2:6264045/1‑91 (MQ=255)
     tattatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTGCATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTagag                                                                                       <  1:6158190/91‑1 (MQ=255)
     tattatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATagtg                                                                                       <  2:4775275/91‑3 (MQ=255)
      attatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAgg                                                                                      <  2:2586422/91‑1 (MQ=255)
       ttatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGc                                                                                     <  2:5900358/91‑1 (MQ=255)
        tataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTAGAGGCt                                                                                    <  2:3209046/91‑1 (MQ=255)
        tataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCt                                                                                    <  1:468383/91‑1 (MQ=255)
         ataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTa                                                                                   <  1:6553527/91‑1 (MQ=255)
         ataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTa                                                                                   >  2:4221272/1‑91 (MQ=255)
         ataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTa                                                                                   >  1:3852155/1‑91 (MQ=255)
          taGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTAGAGGCttg                                                                                  <  1:3730240/91‑3 (MQ=39)
          taGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAt                                                                                  >  2:1909322/1‑91 (MQ=255)
            ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                                <  1:3621390/91‑1 (MQ=255)
            ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                                >  2:3088568/1‑91 (MQ=255)
            ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                                >  1:1825040/1‑91 (MQ=255)
               ttAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTa                                                                             <  2:5210836/91‑1 (MQ=255)
                 aaGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAt                                                                           >  2:5588979/1‑91 (MQ=255)
                 aaGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAt                                                                           <  1:5316338/91‑1 (MQ=255)
                  aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                          <  2:2612877/91‑1 (MQ=255)
                  aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                          >  1:4612553/1‑91 (MQ=255)
                  aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                          <  1:2460717/91‑1 (MQ=255)
                   gTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATta                                                                         <  1:5466212/91‑1 (MQ=255)
                   gTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATta                                                                         <  2:2552798/91‑1 (MQ=255)
                   gTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATta                                                                         <  2:90558/91‑1 (MQ=255)
                    tATGATACTACTTATCACACTATTTATTTATACGTTAATAGGTGACGCCAAAATTCAAGCCTCTACGATTATAGTGGCTATTTTTAATtat                                                                        <  2:4636673/91‑1 (MQ=255)
                      tGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATaa                                                                      >  1:3463240/1‑91 (MQ=255)
                      tGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATaa                                                                      >  2:4315534/1‑91 (MQ=255)
                      tGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATaa                                                                      >  1:2550443/1‑91 (MQ=255)
                       gATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAAt                                                                     >  1:4260764/1‑91 (MQ=255)
                         tactacTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATcc                                                                   >  1:1642398/1‑91 (MQ=255)
                          actactAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCt                                                                  <  2:5987954/91‑1 (MQ=255)
                          actactAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCt                                                                  >  2:6060619/1‑91 (MQ=255)
                           ctactaATCACACTATTTATTTCTACATTACTAGGTGACGCCAAAATTCAAGCCTCTACAATTAGAGAGGCTATTTTTAATTATAATCCTa                                                                 <  2:6920858/91‑1 (MQ=255)
                           ctactaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTa                                                                 <  2:447669/91‑1 (MQ=255)
                           ctactaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTa                                                                 <  2:2129182/91‑1 (MQ=255)
                             actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                               <  2:3619403/91‑1 (MQ=255)
                             actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                               >  2:1200673/1‑91 (MQ=255)
                             actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                               <  1:6329007/91‑1 (MQ=255)
                             actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                               >  1:5964384/1‑91 (MQ=255)
                             actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                               <  1:4979186/91‑1 (MQ=255)
                                aaTCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAAt                                                            >  1:977727/1‑91 (MQ=255)
                                 aTCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATc                                                           >  2:1707023/1‑91 (MQ=255)
                                   cacaCTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCaa                                                         >  2:6705217/1‑91 (MQ=255)
                                   cacaCTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCaa                                                         <  1:5909441/91‑1 (MQ=255)
                                   cacaCTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCaa                                                         >  2:2417385/1‑91 (MQ=255)
                                   cacaCTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCaa                                                         >  2:5139607/1‑91 (MQ=255)
                                     cacTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCaa                                                         <  2:2269385/89‑1 (MQ=255)
                                     cacTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACa                                                       >  2:1054791/1‑91 (MQ=255)
                                       cTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGc                                                     <  1:2953815/91‑1 (MQ=255)
                                          tttatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaaa                                                  >  1:2223220/1‑91 (MQ=255)
                                          tttatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaaa                                                  >  1:5218140/1‑91 (MQ=255)
                                            tatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAc                                                <  1:3150699/91‑1 (MQ=255)
                                            tatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAc                                                <  2:1491331/91‑1 (MQ=255)
                                             atttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCa                                                          >  1:1722784/1‑80 (MQ=255)
                                             atttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAca                                               >  1:6728560/1‑91 (MQ=255)
                                             atttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAca                                               >  1:5110025/1‑91 (MQ=255)
                                              tttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcat                                              >  1:3564823/1‑91 (MQ=255)
                                               ttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatc                                             >  2:1339952/1‑91 (MQ=255)
                                               ttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatc                                             >  2:3438064/1‑91 (MQ=255)
                                                tCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatca                                            >  2:2969290/1‑91 (MQ=255)
                                                 cTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatcat                                           >  1:749920/1‑91 (MQ=255)
                                                  tACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACatcatc                                          >  1:6796929/1‑91 (MQ=255)
                                                   aCATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACAtcatca                                         >  1:5664368/1‑91 (MQ=255)
                                                     aTTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAAt                                       >  1:6225831/1‑91 (MQ=255)
                                                      ttAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATg                                      >  2:386618/1‑91 (MQ=255)
                                                      ttAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATg                                      >  2:6815360/1‑91 (MQ=255)
                                                        aaTAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGGGGCTATTTTTACTTATAATCCTAGCAATCAACAGCAAAACATCATCAATgag                                    <  1:5827293/91‑1 (MQ=255)
                                                        aaTAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATgag                                    <  1:4803997/91‑1 (MQ=255)
                                                           aGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGAtt                                 >  2:4847467/1‑91 (MQ=255)
                                                            ggTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTa                                <  1:6380176/91‑1 (MQ=255)
                                                            ggTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTa                                <  1:4173626/91‑1 (MQ=255)
                                                             gTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAg                               <  2:1620118/91‑1 (MQ=255)
                                                             gTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAg                               <  1:2143519/91‑1 (MQ=255)
                                                              tGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAgg                              >  2:5350807/1‑91 (MQ=255)
                                                              tGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAgg                              >  1:2093311/1‑91 (MQ=255)
                                                               gACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGa                             <  1:4486786/91‑1 (MQ=255)
                                                                  gCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTc                          <  2:3736824/91‑1 (MQ=255)
                                                                   ccAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTcc                         <  1:1356737/91‑1 (MQ=255)
                                                                     aaaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCa                       <  1:5723739/91‑1 (MQ=255)
                                                                     aaaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCa                       <  2:270599/91‑1 (MQ=255)
                                                                     aaaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCa                       <  1:3385662/91‑1 (MQ=255)
                                                                      aaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAg                      >  1:4755/1‑91 (MQ=255)
                                                                       aaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGa                     <  1:1033806/91‑1 (MQ=255)
                                                                        aTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGaa                    <  2:763487/91‑1 (MQ=255)
                                                                        aTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGaa                    <  2:4064525/91‑1 (MQ=255)
                                                                        aTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGaa                    <  2:1574044/91‑1 (MQ=255)
                                                                             aGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAAtata               >  1:2016606/1‑91 (MQ=255)
                                                                              gCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag              >  2:2823810/1‑91 (MQ=255)
                                                                              gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag              >  1:1169205/1‑91 (MQ=255)
                                                                               ccTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagc             >  2:5438770/1‑91 (MQ=255)
                                                                                ctctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagca            <  1:3762798/91‑1 (MQ=255)
                                                                                  ctataaTTGTAGAGGCTATTTTTAATTATAATTCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc          >  2:5767019/5‑91 (MQ=255)
                                                                                   tACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAgcagca         <  2:3371552/91‑1 (MQ=255)
                                                                                    aCAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGcagcag        <  1:273026/91‑1 (MQ=255)
                                                                                    aCAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGcagcag        <  1:39554/91‑1 (MQ=255)
                                                                                     cAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGt       <  1:4020356/91‑1 (MQ=255)
                                                                                      aaTTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTa      >  2:6511639/1‑91 (MQ=255)
                                                                                         tataGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTAAtt   >  2:3902412/1‑91 (MQ=255)
                                                                                          |                                                                                          
TCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTAATTG  >  USA300TCH1516_ALE/722902‑723082

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