Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 722,987 C→T 19.1% T42I (ACA→ATA)  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,9870CT19.1% 309.3 / 61.9 131T42I (ACA→ATA) feuBIron‑uptake system permease protein FeuB
Reads supporting (aligned to +/- strand):  ref base C (47/59);  new base T (10/15);  total (57/74)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 132 total aligned reads displayed.

ATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGT  >  USA300TCH1516_ALE/722901‑723077
                                                                                      |                                                                                          
aTCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAtt                                                                                        <  1:922210/91‑1 (MQ=255)
aTCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAtt                                                                                        <  1:3053040/91‑1 (MQ=255)
 tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttg                                                                                       <  1:4077732/91‑2 (MQ=255)
 tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttg                                                                                       <  1:2100057/91‑2 (MQ=255)
 tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATta                                                                                       >  1:2187967/1‑91 (MQ=255)
  ccTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttgt                                                                                      <  1:826147/91‑3 (MQ=255)
  ccTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttgt                                                                                      <  1:770324/91‑3 (MQ=255)
  ccTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttgt                                                                                      <  1:769954/91‑3 (MQ=255)
  ccTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttgt                                                                                      <  1:1308530/91‑3 (MQ=255)
    tataTTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATag                                                                                    <  1:1068627/91‑1 (MQ=255)
    tataTTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATag                                                                                    <  1:4050987/91‑1 (MQ=255)
    tataTTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATag                                                                                    <  1:3723812/91‑1 (MQ=255)
     atatTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTaga                                                                                   <  1:2214291/91‑1 (MQ=255)
     atatTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTaga                                                                                   <  1:1524201/91‑1 (MQ=255)
      tattatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATagag                                                                                  >  1:2997564/1‑91 (MQ=255)
       attatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTAGAgg                                                                                 <  1:4095201/91‑1 (MQ=255)
        ttatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGc                                                                                <  1:1607687/91‑1 (MQ=255)
         tataGGTTTTAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCt                                                                               <  1:3722443/91‑1 (MQ=255)
         tataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTAGAGGCt                                                                               <  1:3056313/91‑1 (MQ=255)
          ataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTCCATTAATAGGTGACGCCAAAATCCAAGCCTCTACAATTATAGAGGCTa                                                                              <  1:3658214/91‑1 (MQ=255)
             ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                           >  1:1206431/1‑91 (MQ=255)
             ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                           >  1:1977750/1‑91 (MQ=255)
             ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                           >  1:3950486/1‑91 (MQ=255)
              gTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttgt                                                                          >  1:2471770/1‑89 (MQ=255)
                 tAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTaa                                                                       >  1:3669530/1‑91 (MQ=255)
                 tAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTaa                                                                       <  1:544637/91‑1 (MQ=255)
                 tAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTaa                                                                       <  1:544903/91‑1 (MQ=255)
                 tAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTaa                                                                       >  1:3669548/1‑91 (MQ=255)
                   aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                     <  1:361338/91‑1 (MQ=255)
                   aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                     <  1:3474233/91‑1 (MQ=255)
                   aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                     <  1:3474104/91‑1 (MQ=255)
                   aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                     <  1:4081495/91‑1 (MQ=255)
                     tATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATtat                                                                   <  1:2217625/91‑1 (MQ=255)
                      aTGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATtata                                                                  >  1:94209/1‑91 (MQ=255)
                       tGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATaa                                                                 >  1:2605223/1‑91 (MQ=255)
                        gATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAAt                                                                >  1:780885/1‑91 (MQ=255)
                        gATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAAt                                                                >  1:605594/1‑91 (MQ=255)
                            ctactaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTa                                                            <  1:804816/91‑1 (MQ=255)
                            ctactaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTa                                                            >  1:3108259/1‑91 (MQ=255)
                              actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                          >  1:1223124/1‑91 (MQ=255)
                               ctaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCa                                                         >  1:4004359/1‑91 (MQ=255)
                                taATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCaa                                                        >  1:499115/1‑91 (MQ=255)
                                   tCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCa                                                     >  1:2472130/1‑91 (MQ=255)
                                      cacTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACa                                                  <  1:2142210/91‑1 (MQ=255)
                                       acTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAg                                                 <  1:2739560/91‑1 (MQ=255)
                                        cTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGc                                                <  1:3894535/91‑1 (MQ=255)
                                        cTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGc                                                >  1:4016089/1‑91 (MQ=255)
                                         tatttattTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCa                                               <  1:3641660/91‑1 (MQ=255)
                                          atttatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaa                                              <  1:1022090/91‑1 (MQ=255)
                                          atttatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaa                                              <  1:3614078/91‑1 (MQ=255)
                                           tttatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaaa                                             >  1:3011349/1‑91 (MQ=255)
                                             tatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAc                                           <  1:2165456/91‑1 (MQ=255)
                                              atttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAca                                          >  1:2674953/1‑91 (MQ=255)
                                                  cTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatcat                                      >  1:1665022/1‑91 (MQ=255)
                                                    aCATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACAtcatca                                    >  1:3257152/1‑91 (MQ=255)
                                                    aCATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACAtcatca                                    >  1:3257217/1‑91 (MQ=255)
                                                      aTTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAAt                                  <  1:1588837/91‑1 (MQ=255)
                                                      aTTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAAt                                  <  1:1588852/91‑1 (MQ=255)
                                                           tAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGAt                             <  1:148752/91‑1 (MQ=255)
                                                            aGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATGAATGAGAtt                            >  1:3977769/1‑91 (MQ=255)
                                                            aGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGAtt                            <  1:1356985/91‑1 (MQ=255)
                                                             ggTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTa                           <  1:2735966/91‑1 (MQ=255)
                                                             ggTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTa                           >  1:3816363/1‑91 (MQ=255)
                                                                gACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGa                        <  1:2134699/91‑1 (MQ=255)
                                                                 aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                       <  1:478139/91‑1 (MQ=255)
                                                                  cGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAtt                      <  1:3153898/91‑1 (MQ=255)
                                                                   gCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATAAGGATTc                     >  1:2076668/1‑91 (MQ=255)
                                                                      aaaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCa                  <  1:760785/91‑1 (MQ=255)
                                                                       aaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAg                 <  1:443170/91‑1 (MQ=255)
                                                                       aaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAg                 <  1:103323/91‑1 (MQ=255)
                                                                       aaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAg                 >  1:2485333/1‑91 (MQ=255)
                                                                        aaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGa                >  1:656945/1‑91 (MQ=255)
                                                                        aaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGa                >  1:1608987/1‑91 (MQ=255)
                                                                        aaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGa                >  1:2387433/1‑91 (MQ=255)
                                                                        aaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGa                <  1:3671470/91‑1 (MQ=255)
                                                                        aaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGa                >  1:3173006/1‑91 (MQ=255)
                                                                         aTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGaa               <  1:1534489/91‑1 (MQ=255)
                                                                            cAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAata            >  1:2281139/1‑91 (MQ=255)
                                                                             aaGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAatat           <  1:4063541/91‑1 (MQ=255)
                                                                              aGCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAAtata          >  1:2943787/1‑91 (MQ=255)
                                                                               gCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag         >  1:1121480/1‑91 (MQ=255)
                                                                               gCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag         >  1:584323/1‑91 (MQ=255)
                                                                               gCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag         >  1:1830515/1‑91 (MQ=255)
                                                                               gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag         >  1:398674/1‑91 (MQ=255)
                                                                                ccTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagc        <  1:3357659/91‑1 (MQ=255)
                                                                                ccTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagc        <  1:1189519/91‑1 (MQ=255)
                                                                                ccTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagc        >  1:940797/1‑91 (MQ=255)
                                                                                 ctctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagca       >  1:600054/1‑91 (MQ=255)
                                                                                 ctctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagca       >  1:1585694/1‑91 (MQ=255)
                                                                                  tctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcag      >  1:1343654/1‑91 (MQ=255)
                                                                                  tctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcag      >  1:2380271/1‑91 (MQ=255)
                                                                                  tctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcag      >  1:2638698/1‑91 (MQ=255)
                                                                                  tctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcag      >  1:2638603/1‑91 (MQ=255)
                                                                                  tctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcag      <  1:2396670/91‑1 (MQ=255)
                                                                                   ctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc     <  1:1635603/91‑1 (MQ=255)
                                                                                   ctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc     <  1:159744/91‑1 (MQ=255)
                                                                                   ctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc     <  1:3820384/91‑1 (MQ=255)
                                                                                     aCAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGcagcag   <  1:2682855/91‑1 (MQ=255)
                                                                                     aCAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGcagcag   <  1:2682296/91‑1 (MQ=255)
                                                                                      cAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGt  >  1:3474544/1‑91 (MQ=255)
                                                                                      |                                                                                          
ATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGT  >  USA300TCH1516_ALE/722901‑723077

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