Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 3,352:1 +T 5.7% intergenic (+36/‑354) dnaN → / → USA300HOU_RS00015 DNA polymerase III subunit beta/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE3,3521.T 5.7% 899.0 / 43.2 210intergenic (+36/‑354)dnaN/USA300HOU_RS00015DNA polymerase III subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (113/85);  new base T (8/4);  total (121/89)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 244 total aligned reads displayed.

AACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAAAAGACGCC  >  USA300TCH1516_ALE/3261‑3442
                                                                                            |                                                                                          
 aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTATAAATAAATATAAATAAAGGATGACGTGAttaattaa                                                                                             >  2:3054731/1‑91 (MQ=255)
   cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                          <  1:2103671/91‑4 (MQ=38)
   cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                          <  2:1892960/91‑4 (MQ=38)
   cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                          <  2:626609/91‑4 (MQ=38)
     aaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACg                                                                                        >  2:3153235/1‑91 (MQ=255)
      aaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGt                                                                                       >  1:1754546/1‑91 (MQ=255)
       aGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAATGGATGACGTGATTAATTAAAACGTc                                                                                      >  1:3603830/1‑91 (MQ=255)
        ggTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCa                                                                                     >  2:2125397/1‑91 (MQ=255)
        ggTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCa                                                                                     >  2:2241199/1‑91 (MQ=255)
          tGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATc                                                                                   >  1:3825898/1‑91 (MQ=255)
            acgacTCGGTAACGCAATAAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                                 <  2:1173415/91‑1 (MQ=255)
             cgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                                <  1:874030/91‑1 (MQ=255)
             cgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                                <  2:3666082/91‑1 (MQ=255)
             cgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAAAGTCATCCtt                                                                                >  2:3559312/1‑91 (MQ=255)
              gacTCGGTAACGCAATTAATTTTACCAATCAGACCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCttt                                                                               <  1:1179703/91‑1 (MQ=255)
              gacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCttt                                                                               >  2:1524747/1‑91 (MQ=255)
               acTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAGTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTa                                                                              >  2:3002883/1‑91 (MQ=255)
               acTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTa                                                                              <  1:34048/91‑1 (MQ=255)
                  cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                           >  2:234126/1‑91 (MQ=255)
                  cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                           >  2:780857/1‑91 (MQ=255)
                   ggTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTTTAATTAATCACGTCATCCTTTAttta                                                                          <  2:1114343/91‑2 (MQ=21)
                    gTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTTTAATTAATCACGTCATCCTTTAtttat                                                                         <  2:1557554/91‑3 (MQ=21)
                     tAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttttt                                                                        <  1:1711303/91‑1 (MQ=255)
                      aaCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTTTAATTAATCACGTCATCCTTTAtttatat                                                                       <  1:650938/91‑5 (MQ=18)
                       aCGCAATTACTTCTGCCACTCAGATCTTTCGAACACTAAATATACATGAAGGATGACGTGATCAATTAAAGCGTCATCCTTTATTTTTTgg                                                                      <  2:1670221/91‑1 (MQ=255)
                       aCGCAATTAATTTTACGAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                                      <  2:3707865/91‑1 (MQ=255)
                       aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                                      >  2:3725768/1‑91 (MQ=255)
                       aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                                      <  2:1306967/91‑1 (MQ=255)
                       aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                                      >  2:1422614/1‑91 (MQ=255)
                       aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                                      >  1:2658878/1‑91 (MQ=255)
                        cGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGc                                                                     <  1:1219687/91‑1 (MQ=255)
                          cAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaa                                                                   >  2:1568647/1‑91 (MQ=255)
                          cAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaa                                                                   >  1:1130296/1‑91 (MQ=255)
                            attaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaa                                                                 >  1:111750/1‑91 (MQ=255)
                            attaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaa                                                                 <  1:413603/91‑1 (MQ=255)
                             ttaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaaa                                                                >  1:1114657/1‑91 (MQ=255)
                             ttaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaaa                                                                <  1:2199768/91‑1 (MQ=255)
                               aattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAATTTTGCCaaaaaata                                                             <  2:873107/91‑4 (MQ=255)
                                  tttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAAtt                                                           >  1:2885187/1‑91 (MQ=255)
                                   ttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTc                                                          >  2:3864246/1‑91 (MQ=255)
                                     aCCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                                >  2:2822517/1‑67 (MQ=255)
                                      ccAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAg                                                       >  2:1117063/1‑91 (MQ=255)
                                      ccAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAg                                                       <  2:921760/91‑1 (MQ=255)
                                       cAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAgg                                                      <  2:2170461/91‑1 (MQ=255)
                                       cAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAgg                                                      <  1:2554501/91‑1 (MQ=255)
                                           cAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCg                                                  <  1:2615323/91‑1 (MQ=255)
                                              aaCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                                        <  1:1859494/82‑1 (MQ=255)
                                              aaCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTAt                                               >  1:55307/1‑91 (MQ=255)
                                               aCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATg                                              <  1:1245052/91‑1 (MQ=255)
                                                   aCTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaa                                          <  1:2481248/91‑1 (MQ=255)
                                                    cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                         <  1:1568593/91‑1 (MQ=255)
                                                    cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                         >  1:315520/1‑91 (MQ=255)
                                                    cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                         <  2:2398306/91‑1 (MQ=255)
                                                      aaaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaata                                       <  2:325331/91‑1 (MQ=255)
                                                       aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                                      >  1:2369707/1‑91 (MQ=255)
                                                       aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                                      <  1:1524747/91‑1 (MQ=255)
                                                        aaataaatATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataaa                                     >  2:3774136/1‑91 (MQ=255)
                                                        aaataaatATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATATTTTGGCAAAAATAATTCTAGAAGCGAAAGTAaaataaa                                     >  2:3376461/1‑91 (MQ=255)
                                                          ataaatataaatAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAtt                                   >  1:2058033/1‑91 (MQ=255)
                                                          ataaatataaatAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAtt                                   >  1:1876445/1‑91 (MQ=255)
                                                           taaatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAttt                                  <  2:1834285/91‑1 (MQ=255)
                                                            aaatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTg                                 <  1:2948453/91‑1 (MQ=255)
                                                             aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                                <  1:2423620/91‑1 (MQ=255)
                                                             aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                                <  1:344949/91‑1 (MQ=255)
                                                              atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                           <  2:2864121/79‑1 (MQ=255)
                                                              atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                           >  1:2864121/1‑79 (MQ=255)
                                                              atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgc                               <  2:1249885/91‑1 (MQ=255)
                                                               tataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgca                              <  2:3799819/91‑1 (MQ=255)
                                                                ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                             >  1:2012408/1‑91 (MQ=255)
                                                                ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                             >  2:3636107/1‑91 (MQ=255)
                                                                ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                             <  2:2219053/91‑1 (MQ=255)
                                                                ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                             >  1:1973675/1‑91 (MQ=255)
                                                                   aataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAt                          >  2:2704723/1‑91 (MQ=255)
                                                                    ataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAtt                         <  1:3733924/91‑1 (MQ=255)
                                                                      aaaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAtttt                       >  1:2263683/1‑91 (MQ=21)
                                                                       aaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa                      <  2:1712068/91‑1 (MQ=255)
                                                                       aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa                      >  1:387791/1‑91 (MQ=25)
                                                                       aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa                      >  1:1354653/1‑91 (MQ=25)
                                                                       aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa                      >  1:921903/1‑91 (MQ=25)
                                                                       aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa                      >  1:2817271/1‑91 (MQ=25)
                                                                        aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa                     >  2:272411/1‑91 (MQ=255)
                                                                        aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa                     <  1:2222238/91‑1 (MQ=255)
                                                                        aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa                     <  2:909137/91‑1 (MQ=255)
                                                                             gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc                >  2:1162925/1‑91 (MQ=255)
                                                                             gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc                >  1:2720441/1‑91 (MQ=255)
                                                                             gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc                >  2:1163174/1‑91 (MQ=255)
                                                                               cGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCaa              >  2:457899/1‑91 (MQ=255)
                                                                                 gttttaattaatCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATACATTTGGCAGCATTTTAAACAGCaaa             >  2:2465737/4‑91 (MQ=38)
                                                                                  gATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATa           >  1:3067461/1‑91 (MQ=255)
                                                                                  gATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATa           <  1:3802731/91‑1 (MQ=255)
                                                                                   attaattaAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATaa          >  2:21822/1‑91 (MQ=255)
                                                                                   attaattaAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATaa          >  2:880797/1‑91 (MQ=255)
                                                                                   attaattaAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATaa          >  2:22061/1‑91 (MQ=255)
                                                                                   attaattaAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATaa          >  1:100595/1‑91 (MQ=255)
                                                                                    ttaattaaAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATaaa         >  2:1746034/1‑91 (MQ=255)
                                                                                       attaaAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAAAAGa      >  2:1690399/1‑91 (MQ=255)
                                                                                       attaaAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAAAAGa      >  2:522938/1‑91 (MQ=255)
                                                                                       attaaAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAAAAGa      >  1:768318/1‑91 (MQ=255)
                                                                                         taaAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAAAAGACg    >  2:1607338/1‑91 (MQ=255)
                                                                                          aaaaCGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAAAAGACGc   >  1:23381/1‑91 (MQ=255)
                                                                                            |                                                                                          
AACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAAAAGACGCC  >  USA300TCH1516_ALE/3261‑3442

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