Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 722,987 C→T 12.4% 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,9870CT12.4% 646.7 / 61.0 226T42I (ACA→ATA) feuBIron‑uptake system permease protein FeuB
Reads supporting (aligned to +/- strand):  ref base C (93/105);  new base T (9/19);  total (102/124)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 232 total aligned reads displayed.

TTTTATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGT  >  USA300TCH1516_ALE/722897‑723077
                                                                                          |                                                                                          
ttttATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTAc                                                                                            <  2:2686384/91‑1 (MQ=255)
  ttATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACaa                                                                                          >  1:513953/1‑91 (MQ=255)
  ttATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACaa                                                                                          >  1:737836/1‑91 (MQ=255)
  ttATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACaa                                                                                          >  1:900765/1‑91 (MQ=255)
  ttATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACaa                                                                                          <  1:2727759/91‑1 (MQ=255)
  ttATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACaa                                                                                          <  1:1846989/91‑1 (MQ=255)
   tATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAt                                                                                         <  1:1171850/91‑1 (MQ=255)
   tATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAt                                                                                         <  2:2117643/91‑1 (MQ=255)
   tATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAt                                                                                         <  1:1617097/91‑1 (MQ=255)
   tATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGGGACGCCAAAATTCAAGCCTCTATAAt                                                                                         <  2:2247755/91‑1 (MQ=255)
     tCCTATATTATAGGTTTAAGTATGATACTACTAGTCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttg                                                                                       <  1:3503650/91‑2 (MQ=255)
     tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttg                                                                                       <  2:8650/91‑2 (MQ=255)
     tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttg                                                                                       <  2:733533/91‑2 (MQ=255)
     tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttg                                                                                       <  2:3212249/91‑2 (MQ=255)
     tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATta                                                                                       >  1:1025363/1‑91 (MQ=255)
     tCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATta                                                                                       >  1:1870737/1‑91 (MQ=255)
     tCCTATATTATAGGTTTAAGTATGAGACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAAttg                                                                                       <  2:2731864/91‑2 (MQ=255)
       cTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATtata                                                                                     >  1:2057766/1‑91 (MQ=255)
        tataTTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTag                                                                                    <  2:1842356/91‑1 (MQ=255)
        tataTTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATag                                                                                    <  1:2533685/91‑1 (MQ=255)
         atatTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATaga                                                                                   <  1:65608/91‑1 (MQ=255)
          tattatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTagag                                                                                  <  1:2641797/91‑1 (MQ=255)
          tattatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATagag                                                                                  >  2:3904478/1‑91 (MQ=255)
           attatAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAgg                                                                                 <  2:2402380/91‑1 (MQ=255)
             tataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTATAATTGTAGAGGCt                                                                               <  1:28123/91‑1 (MQ=255)
             tataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTATACAATTATAGAGGCt                                                                               <  2:3170997/91‑1 (MQ=255)
              ataGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTa                                                                              >  1:1541044/1‑91 (MQ=255)
                aGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAtt                                                                            <  1:83737/91‑1 (MQ=255)
                 ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                           >  1:3859346/1‑91 (MQ=255)
                 ggTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTAttt                                                                           >  2:2905026/1‑91 (MQ=255)
                    ttAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTCCAATTATAGAGGCTATTTTTa                                                                        >  2:2789230/1‑91 (MQ=255)
                    ttAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTa                                                                        <  1:2642467/91‑1 (MQ=255)
                    ttAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTa                                                                        <  1:3829484/91‑1 (MQ=255)
                      aaGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAt                                                                      <  1:2763588/91‑1 (MQ=255)
                      aaGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAt                                                                      <  2:3276369/91‑1 (MQ=255)
                       aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                     >  2:1659174/1‑91 (MQ=255)
                       aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                     <  2:347391/91‑1 (MQ=255)
                       aGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAAtt                                                                     <  1:469216/91‑1 (MQ=255)
                        gTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATta                                                                    >  2:1549826/1‑91 (MQ=255)
                        gTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATta                                                                    >  1:2406062/1‑91 (MQ=255)
                           tGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACCATTATAGAGGCTATTTTTAATTATaa                                                                 >  2:2417359/1‑91 (MQ=255)
                            gATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAAt                                                                <  2:1634656/91‑1 (MQ=255)
                               actactAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCt                                                             <  1:3700188/91‑1 (MQ=255)
                                ctactaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTa                                                            >  1:106404/1‑91 (MQ=255)
                                 tactaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAg                                                           <  2:3847977/91‑1 (MQ=255)
                                  actaATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGc                                                          <  1:2008872/91‑1 (MQ=255)
                                     aaTCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGCGGCTATTTTTAATTATAATCCTAGCaac                                                       >  2:3929477/1‑90 (MQ=255)
                                        cacaCTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCaa                                                    >  2:2897544/1‑91 (MQ=255)
                                        cacaCTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCaa                                                    <  2:1293723/91‑1 (MQ=255)
                                         acacTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGcaaccaac                                                   >  1:3298034/1‑91 (MQ=255)
                                           acTATTTATTTCCACATTAATAGGTGACGCCAAAATTCAAGCCTCTCCAATTATAGAGGCTCTTCTCCATTATAATCCTAGCCATCAACAg                                                 >  2:109676/1‑91 (MQ=255)
                                            cTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGc                                                <  1:1738587/91‑1 (MQ=255)
                                              atttatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCaa                                              <  2:513060/91‑1 (MQ=255)
                                                 tatttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAc                                           <  1:325012/91‑1 (MQ=255)
                                                  atttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAca                                          <  2:2227921/91‑1 (MQ=255)
                                                  atttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAca                                          >  1:3937178/1‑91 (MQ=255)
                                                    ttCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatc                                        >  2:1361522/1‑91 (MQ=255)
                                                     tCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatca                                       >  1:773151/1‑91 (MQ=255)
                                                      cTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAAcatcat                                      <  2:947629/91‑1 (MQ=255)
                                                       tACATCAAGAGGTGACGCCAAAATTCAAGCCTCTACAATTAGAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACatcatc                                     <  2:3643927/91‑1 (MQ=255)
                                                         cATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCaa                                   >  1:1299903/1‑91 (MQ=255)
                                                             aaTAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATgag                               <  1:527372/91‑1 (MQ=255)
                                                             aaTAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATgag                               <  1:2345136/91‑1 (MQ=255)
                                                               tAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGAt                             <  2:3076804/91‑1 (MQ=255)
                                                               tAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAAGCAACAGCAAAACATCATCAATGAGAt                             >  2:1155076/1‑91 (MQ=255)
                                                                   tGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAgg                         <  1:1212112/91‑1 (MQ=255)
                                                                   tGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAgg                         <  1:3765800/91‑1 (MQ=255)
                                                                    gACGCCAAAATTCAAGCCTCTACAATTATAGAGTCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGa                        >  1:2095736/1‑91 (MQ=255)
                                                                    gACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGa                        <  1:2118822/91‑1 (MQ=255)
                                                                    gACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGa                        <  1:2118506/91‑1 (MQ=255)
                                                                     aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                       >  2:1613958/1‑91 (MQ=255)
                                                                     aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                       >  1:1349147/1‑91 (MQ=255)
                                                                     aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                       <  2:1031201/91‑1 (MQ=255)
                                                                     aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                       >  1:3342146/1‑91 (MQ=255)
                                                                     aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                       >  1:1647939/1‑91 (MQ=255)
                                                                     aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                       <  2:786351/91‑1 (MQ=255)
                                                                     aCGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAt                       <  1:598929/91‑1 (MQ=255)
                                                                      cGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGAtt                      <  1:1748911/91‑1 (MQ=255)
                                                                       gCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTc                     >  1:1257456/1‑91 (MQ=255)
                                                                          aaaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCa                  <  1:2758242/91‑1 (MQ=255)
                                                                           aaaTTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAg                 >  2:259160/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)
                                                                                 aaGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAgtat           >  1:3117656/1‑91 (MQ=255)
                                                                                 aaGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAatat           <  2:505832/91‑1 (MQ=255)
                                                                                  aGCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAAtata          >  1:2690694/1‑91 (MQ=255)
                                                                                   gCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag         >  2:2838971/1‑91 (MQ=255)
                                                                                   gCCTCTATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag         >  1:3534439/1‑91 (MQ=255)
                                                                                   gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag         >  1:1449427/1‑91 (MQ=255)
                                                                                   gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag         >  1:1292529/1‑91 (MQ=255)
                                                                                   gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag         >  2:616015/1‑91 (MQ=255)
                                                                                   gCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATag         >  2:88070/1‑91 (MQ=255)
                                                                                    ccTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagc        <  1:2557407/91‑1 (MQ=255)
                                                                                     ctctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagca       >  1:2312655/1‑91 (MQ=255)
                                                                                      tctATAATTGTAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAACATagcag      >  2:2453224/1‑91 (MQ=255)
                                                                                       ctACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATagcagc     <  1:1613958/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)
                                                                                          |                                                                                          
TTTTATCCTATATTATAGGTTTAAGTATGATACTACTAATCACACTATTTATTTCTACATTAATAGGTGACGCCAAAATTCAAGCCTCTACAATTATAGAGGCTATTTTTAATTATAATCCTAGCAATCAACAGCAAAACATCATCAATGAGATTAGGATTCCCAGAAATATAGCAGCAGT  >  USA300TCH1516_ALE/722897‑723077

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: ‑