Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,245,337 A→T 8.8% intergenic (+142/‑101) fabG → / → acpP 3‑oxoacyl‑[acyl‑carrier‑protein] reductase FabG/Acyl carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,245,3370AT 8.8% 388.8 / 21.8 114intergenic (+142/‑101)fabG/acpP3‑oxoacyl‑[acyl‑carrier‑protein] reductase FabG/Acyl carrier protein
Reads supporting (aligned to +/- strand):  ref base A (57/47);  new base T (3/7);  total (60/54)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 118 total aligned reads displayed.

GAATTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAAGGA  >  USA300TCH1516_ALE/1245248‑1245426
                                                                                         |                                                                                         
gaATTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATaa                                                                                          >  1:2191964/1‑91 (MQ=255)
 aaTTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATaaa                                                                                         >  1:1266876/1‑91 (MQ=255)
 aaTTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATaaa                                                                                         >  1:1497934/1‑91 (MQ=255)
 aaTTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATaaa                                                                                         >  1:1526238/1‑91 (MQ=255)
   ttGTTTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAAtt                                                                                       >  1:1657039/1‑91 (MQ=255)
    tGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGAGAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTg                                                                                      >  1:3303576/1‑91 (MQ=255)
    tGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGAGAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTg                                                                                      >  1:3303596/1‑91 (MQ=255)
     gTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGa                                                                                     <  1:2265422/91‑1 (MQ=255)
     gTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGa                                                                                     >  1:726504/1‑91 (MQ=255)
       cTGACCTAGTCAACTTTGCGGGGGAAATTCTATGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATAAATTGATa                                                                                   <  1:3543466/91‑1 (MQ=255)
       cTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGAtt                                                                                   <  1:3814056/91‑2 (MQ=38)
       cTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATa                                                                                   >  1:1110863/1‑91 (MQ=255)
       cTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATa                                                                                   >  1:662832/1‑91 (MQ=255)
       cTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATa                                                                                   <  1:3709189/91‑1 (MQ=255)
        tGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATaa                                                                                  >  1:2134785/1‑91 (MQ=255)
         gACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGAttag                                                                                 <  1:3025930/91‑4 (MQ=38)
          aCCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGATTAgt                                                                                <  1:958421/91‑1 (MQ=38)
          aCCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGATTAgt                                                                                <  1:1022266/91‑1 (MQ=38)
          aCCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAgt                                                                                <  1:3823589/91‑1 (MQ=255)
          aCCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAgt                                                                                >  1:1741394/1‑91 (MQ=255)
           ccTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGATTAgtg                                                                               <  1:4028479/91‑1 (MQ=37)
           ccTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGATTAgtg                                                                               <  1:534164/91‑1 (MQ=37)
           ccTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTTATCAATTTATTGATTAgtg                                                                               <  1:893416/91‑1 (MQ=37)
              aGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGttt                                                                            <  1:3104317/91‑1 (MQ=255)
                      ttGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAAtt                                                                    <  1:988092/91‑1 (MQ=255)
                        gCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTc                                                                  <  1:166665/91‑1 (MQ=255)
                          gggggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTa                                                                <  1:45204/91‑1 (MQ=255)
                          gggggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTa                                                                <  1:45353/91‑1 (MQ=255)
                           ggggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTAc                                                               <  1:1036146/91‑1 (MQ=255)
                           ggggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTAc                                                               <  1:3871171/91‑1 (MQ=255)
                           ggggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTAc                                                               >  1:2514471/1‑91 (MQ=255)
                           ggggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTAc                                                               >  1:2592598/1‑91 (MQ=255)
                            gggAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACt                                                              <  1:2735549/91‑1 (MQ=255)
                               aaaTTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGt                                                           <  1:3475757/91‑1 (MQ=255)
                               aaaTTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGt                                                           <  1:3475747/91‑1 (MQ=255)
                               aaaTTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGt                                                           <  1:3130949/91‑1 (MQ=255)
                                  ttCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGtttt                                                        <  1:2871671/91‑1 (MQ=255)
                                   tCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGttttt                                                       <  1:522372/91‑1 (MQ=255)
                                   tCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGttttt                                                       <  1:1090951/91‑1 (MQ=255)
                                    cTAAGCAACCTGGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGtttttt                                                      <  1:145491/91‑1 (MQ=255)
                                     tAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTa                                                     <  1:1330015/91‑1 (MQ=255)
                                     tAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTa                                                     >  1:2279326/1‑91 (MQ=255)
                                       aGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGa                                                   >  1:640646/1‑91 (MQ=255)
                                       aGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGa                                                   <  1:2262373/91‑1 (MQ=255)
                                        gCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGaa                                                  <  1:3630219/91‑1 (MQ=255)
                                        gCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGaa                                                  <  1:3629231/91‑1 (MQ=255)
                                        gCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGaa                                                  <  1:1274922/91‑1 (MQ=255)
                                           aCCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAAttt                                               >  1:411454/1‑91 (MQ=255)
                                            ccTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTa                                              <  1:1836586/91‑1 (MQ=255)
                                             cTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTaa                                             >  1:1273661/1‑91 (MQ=255)
                                             cTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTaa                                             >  1:1136486/1‑91 (MQ=255)
                                             cTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTaa                                             >  1:2205043/1‑91 (MQ=255)
                                                gATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAAt                                          >  1:4064917/1‑91 (MQ=255)
                                                gATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAAt                                          >  1:3099793/1‑91 (MQ=255)
                                                gATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAAt                                          >  1:360674/1‑91 (MQ=255)
                                                  tAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATgg                                        >  1:2677024/1‑91 (MQ=255)
                                                  tAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATgg                                        >  1:697425/1‑91 (MQ=255)
                                                  tAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATgg                                        <  1:3633087/91‑1 (MQ=255)
                                                     ggTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGaa                                     <  1:3557026/91‑1 (MQ=255)
                                                      gTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGaaa                                    >  1:2813671/1‑91 (MQ=255)
                                                      gTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGaaa                                    >  1:2813659/1‑91 (MQ=255)
                                                        tCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAat                                  >  1:3942414/1‑91 (MQ=255)
                                                         ccAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAata                                 >  1:2432658/1‑91 (MQ=255)
                                                          cAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAatat                                >  1:1908275/1‑91 (MQ=255)
                                                            gAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATAtag                              >  1:2633437/1‑91 (MQ=255)
                                                              aTTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATAtagta                            >  1:373218/1‑91 (MQ=255)
                                                              aTTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATAtagta                            >  1:2762734/1‑91 (MQ=255)
                                                               tttCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATAtagtag                           >  1:843587/1‑91 (MQ=255)
                                                               tttCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATAtagtag                           <  1:2608940/91‑1 (MQ=255)
                                                               tttCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATAtagtag                           <  1:2609354/91‑1 (MQ=255)
                                                                 tctcCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTc                         <  1:406906/91‑1 (MQ=255)
                                                                  ctcCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCt                        >  1:1186871/1‑91 (MQ=255)
                                                                   tcCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTa                       >  1:2336939/1‑91 (MQ=255)
                                                                    cccTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTAt                      >  1:2613411/1‑91 (MQ=255)
                                                                      cTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGt                    >  1:4043584/1‑91 (MQ=255)
                                                                       ttaGAAACACTTATCAATTTATTGATTAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGta                   >  1:2750642/3‑91 (MQ=25)
                                                                       tAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAATAGTCTATGta                   >  1:1611082/1‑91 (MQ=255)
                                                                         aGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGtata                 <  1:2006305/91‑1 (MQ=255)
                                                                         aGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGtata                 <  1:3159652/91‑1 (MQ=255)
                                                                          gaagcacTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAg                <  1:3271895/87‑1 (MQ=255)
                                                                            aaCACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGc              <  1:4025889/91‑1 (MQ=255)
                                                                            aaCACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGc              <  1:4026075/91‑1 (MQ=255)
                                                                            aaCACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGc              <  1:1528297/91‑1 (MQ=255)
                                                                             acacTAATCAATAAATTGATAAGTGATTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCa             <  1:3863556/91‑1 (MQ=255)
                                                                              cacTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAt            >  1:266989/1‑91 (MQ=255)
                                                                              cacTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAt            >  1:2192713/1‑91 (MQ=255)
                                                                              cacTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAt            >  1:2905551/1‑91 (MQ=255)
                                                                                cTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAttt          <  1:1652729/91‑1 (MQ=255)
                                                                                cTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAttt          <  1:400633/91‑1 (MQ=255)
                                                                                cTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAttt          <  1:942398/91‑1 (MQ=255)
                                                                                cTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAttt          >  1:1406743/1‑91 (MQ=255)
                                                                                cTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAttt          >  1:2908402/1‑91 (MQ=255)
                                                                                  taTCAATTTATTGATTAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCAttttt        >  1:3307782/2‑91 (MQ=39)
                                                                                   aTCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTa       >  1:1242994/1‑91 (MQ=255)
                                                                                    tCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTaa      >  1:3087908/1‑91 (MQ=255)
                                                                                    tCAATAAATTGATAAGTATTTCTAAAATTTCTTCTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTaa      <  1:3553016/91‑1 (MQ=255)
                                                                                      aaTAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAag    >  1:1590659/1‑91 (MQ=255)
                                                                                      aaTAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAag    <  1:604793/91‑1 (MQ=255)
                                                                                        tAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAagga  >  1:381907/1‑91 (MQ=255)
                                                                                        tAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAagga  >  1:382340/1‑91 (MQ=255)
                                                                                         |                                                                                         
GAATTGTCTGACCTAGTCAACTTTGCGGGGGAAATTCTAAGCAACCTAGATAAGGTTCCAGAATTTCTCCCTAAGAAACACTAATCAATAAATTGATAAGTGTTTCTAAAATTTCTACTTGTTTTTTAGAATTTAAAATGGGAAAATATAGTAGTCTATGTATAGGCATTTTTAAAGGA  >  USA300TCH1516_ALE/1245248‑1245426

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: