Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 722,992 A→G 12.0% 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,9920AG12.0% 574.5 / 63.7 216I44V (ATA→GTA) feuBIron‑uptake system permease protein FeuB
Reads supporting (aligned to +/- strand):  ref base A (83/107);  new base G (19/7);  total (102/114)
Fisher's exact test for biased strand distribution p-value = 6.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

Only 100 of 228 total aligned reads displayed.

TCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTAATTG  >  USA300TCH1516_ALE/722902‑723082
                                                                                          |                                                                                          
tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATta                                                                                            >  1:1870737/1‑91 (MQ=255)
tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATta                                                                                            >  1:3918868/1‑91 (MQ=255)
tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATta                                                                                            >  1:1025363/1‑91 (MQ=255)
  cTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATtata                                                                                          >  1:2057766/1‑91 (MQ=255)
   tataTTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTag                                                                                         <  1:1294794/91‑1 (MQ=255)
    atatTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATaga                                                                                        <  1:65608/91‑1 (MQ=255)
     tattatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTagag                                                                                       <  1:2641857/91‑1 (MQ=255)
     tattatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATagag                                                                                       <  1:1252494/91‑1 (MQ=255)
     tattatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATagag                                                                                       >  2:3904478/1‑91 (MQ=255)
       ttatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGc                                                                                     <  2:1086888/91‑1 (MQ=255)
        tataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTAGAGGCt                                                                                    <  1:1191451/91‑1 (MQ=255)
        tataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTAGAGGCt                                                                                    <  1:28123/91‑1 (MQ=255)
        tataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCt                                                                                    <  1:869494/91‑1 (MQ=255)
         ataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTAGAGGCtt                                                                                   <  1:3855707/91‑2 (MQ=255)
          taGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAt                                                                                  >  1:2510505/1‑91 (MQ=255)
            ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                                >  2:321061/1‑91 (MQ=255)
               ttAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTa                                                                             >  2:2787152/1‑91 (MQ=255)
                tAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTaa                                                                            <  2:2827652/91‑1 (MQ=255)
                 aaGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAt                                                                           <  2:3276369/91‑1 (MQ=255)
                 aaGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAt                                                                           <  1:1773867/91‑1 (MQ=255)
                 aaGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAt                                                                           <  1:2763588/91‑1 (MQ=255)
                   gTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATta                                                                         >  2:1549826/1‑91 (MQ=255)
                   gTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATta                                                                         >  1:2406062/1‑91 (MQ=255)
                      tGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAACATTCAAGCCTTTTCACTTATAAAGGCCATTTTTAATTATaa                                                                      >  1:2655373/1‑91 (MQ=255)
                      tGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACCATTATAGAGGCTATTTTTAATTATaa                                                                      >  2:2417359/1‑91 (MQ=255)
                      tGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATaa                                                                      >  1:3773929/1‑91 (MQ=255)
                      tGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATaa                                                                      >  1:3489509/1‑91 (MQ=255)
                        aTACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATc                                                                    <  1:506208/91‑1 (MQ=255)
                           ctactaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTa                                                                 >  1:106404/1‑91 (MQ=255)
                            tactaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAg                                                                <  2:3711156/91‑1 (MQ=255)
                            tactaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAg                                                                <  2:3847977/91‑1 (MQ=255)
                             actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                               >  1:400037/1‑91 (MQ=255)
                             actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                               <  1:2008872/91‑1 (MQ=255)
                             actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                               <  2:1133835/91‑1 (MQ=255)
                             actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                               >  2:122916/1‑91 (MQ=255)
                                aaTCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGCGGCTATTTTTAATTATAATCCTAGCaac                                                            >  2:3929477/1‑90 (MQ=255)
                                 aTCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATc                                                           >  1:1723169/1‑91 (MQ=255)
                                 aTCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATc                                                           >  1:3759238/1‑91 (MQ=255)
                                  tCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCa                                                          >  1:2742855/1‑91 (MQ=255)
                                   cacaCTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCaa                                                         >  2:2897544/1‑91 (MQ=255)
                                   cacaCTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCaa                                                         <  2:1293723/91‑1 (MQ=255)
                                       cTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGc                                                     <  1:1738769/91‑1 (MQ=255)
                                       cTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGc                                                     >  2:2594342/1‑91 (MQ=255)
                                           ttatttCTACATTAATAGGTGACGCCAAATTGCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaaaa                                                 <  1:3166113/91‑1 (MQ=255)
                                           ttatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaaaa                                                 >  2:836129/1‑91 (MQ=255)
                                            tatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAc                                                <  1:325012/91‑1 (MQ=255)
                                             atttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACGATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAca                                               <  2:3463036/91‑1 (MQ=255)
                                             atttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAca                                               <  2:2227921/91‑1 (MQ=255)
                                               ttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatc                                             <  2:1651922/91‑1 (MQ=255)
                                                tCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatca                                            <  1:1258237/91‑1 (MQ=255)
                                                 cTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatcat                                           <  2:947629/91‑1 (MQ=255)
                                                  tACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACatcatc                                          >  1:1300428/1‑91 (MQ=255)
                                                  tACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACatcatc                                          <  1:1598371/91‑1 (MQ=255)
                                                  tACATCAAGAGGTGACGCCAAAATTCAAGCCTCTACAATTAGAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACatcatc                                          <  2:3643927/91‑1 (MQ=255)
                                                        aaTAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATgag                                    <  1:2345136/91‑1 (MQ=255)
                                                        aaTAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATgag                                    >  1:1480127/1‑91 (MQ=255)
                                                        aaTAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATgag                                    <  2:528214/91‑1 (MQ=255)
                                                        aaTAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATgag                                    <  1:527317/91‑1 (MQ=255)
                                                        aaTAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATgag                                    <  2:1007182/91‑1 (MQ=255)
                                                          tAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGAt                                  <  2:3076804/91‑1 (MQ=255)
                                                           aGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGAtt                                 >  1:3572474/1‑91 (MQ=255)
                                                              tGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAgg                              >  2:771895/1‑91 (MQ=255)
                                                               gACGCCAAAATTCAAGCCTCTACAATTATAGAGTCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGa                             >  1:2095736/1‑91 (MQ=255)
                                                               gACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGa                             <  1:2118506/91‑1 (MQ=255)
                                                               gACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGa                             <  1:2118822/91‑1 (MQ=255)
                                                                aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                            >  1:1349147/1‑91 (MQ=255)
                                                                aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                            >  1:1647939/1‑91 (MQ=255)
                                                                aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                            >  1:700062/1‑91 (MQ=255)
                                                                aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                            >  1:3342146/1‑91 (MQ=255)
                                                                aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                            >  2:1614177/1‑91 (MQ=255)
                                                                 cGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAtt                           >  2:240007/1‑91 (MQ=255)
                                                                     aaaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCa                       <  1:2758498/91‑1 (MQ=255)
                                                                     aaaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCa                       <  2:2586905/91‑1 (MQ=255)
                                                                      aaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAg                      >  1:2276185/1‑91 (MQ=255)
                                                                        aTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGaa                    <  2:1649654/91‑1 (MQ=255)
                                                                          tCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAat                  <  2:1898479/91‑1 (MQ=255)
                                                                           cAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAata                 >  2:3303545/1‑91 (MQ=255)
                                                                           cAAGCCTCTACAATTATAGAGGCTATTTATAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAata                 >  2:2923807/1‑91 (MQ=255)
                                                                            aaGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAatat                <  2:505832/91‑1 (MQ=255)
                                                                              gCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag              >  1:3534439/1‑91 (MQ=255)
                                                                              gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag              >  2:88070/1‑91 (MQ=255)
                                                                              gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag              >  2:616015/1‑91 (MQ=255)
                                                                              gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag              <  2:3141816/91‑1 (MQ=255)
                                                                               ccTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagc             <  1:2557407/91‑1 (MQ=255)
                                                                                 tctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAACATagcag           >  2:2453224/1‑91 (MQ=255)
                                                                                  ctataaTTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAATACATCATCAATGAGATTAGGATTCACAGAAATATagcagc          >  2:199329/5‑91 (MQ=255)
                                                                                  ctataaTTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc          >  2:3792064/5‑91 (MQ=255)
                                                                                  ctataaTTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc          >  1:2592870/5‑91 (MQ=255)
                                                                                  ctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc          <  1:1613958/91‑1 (MQ=255)
                                                                                  ctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc          <  2:1724789/91‑1 (MQ=255)
                                                                                    aCAATTATCGAGGCTATTTTTAATTAGAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGcagcag        <  2:3431629/91‑1 (MQ=255)
                                                                                    aCAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGcagcag        >  2:3585731/1‑91 (MQ=255)
                                                                                     cAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGt       <  2:402819/91‑1 (MQ=255)
                                                                                      aaTTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTa      >  1:3278412/1‑91 (MQ=255)
                                                                                      aaTTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTa      >  1:59454/1‑91 (MQ=255)
                                                                                      aaTTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACACCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTa      <  1:665127/91‑1 (MQ=255)
                                                                                       aTTGTAGAGGCTATTTCGAATTATAATCCTAGCAATCCACAGCAACACATAATAGATGAGATAAGGATTCCTAGAAATATAGCAGGAGTaa     >  2:1427192/1‑91 (MQ=255)
                                                                                        ttATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTAAt    <  2:774483/91‑1 (MQ=255)
                                                                                         tataGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAAAATAGCAGCAGTAAtt   >  1:783507/1‑91 (MQ=255)
                                                                                          ataGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGTAATTg  >  2:447393/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: