Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 3,344 A→T 10.8% intergenic (+28/‑362) 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,3440AT10.8% 694.4 / 58.8 222intergenic (+28/‑362)dnaN/USA300HOU_RS00015DNA polymerase III subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (110/88);  new base T (8/16);  total (118/104)
Fisher's exact test for biased strand distribution p-value = 5.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 243 total aligned reads displayed.

ATTCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAA  >  USA300TCH1516_ALE/3254‑3434
                                                                                          |                                                                                           
 ttCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGat                                                                                            <  1:3585958/91‑1 (MQ=255)
  tCTAAAACCAAAAGGAGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGatt                                                                                           <  2:3639272/91‑1 (MQ=255)
     aaaaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaat                                                                                        <  2:2547316/91‑1 (MQ=255)
      aaaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaatt                                                                                       >  1:1460003/1‑91 (MQ=255)
       aaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaatta                                                                                      <  1:1095804/91‑1 (MQ=255)
       aaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaa                                                                                     <  1:2309609/91‑1 (MQ=255)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTATAAATAAATATAAATAAAGGATGACGTGAttaattaa                                                                                     >  2:3054731/1‑91 (MQ=255)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGAttaattaa                                                                                     <  1:1162347/91‑1 (MQ=255)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaat                                                                                    <  1:1628299/91‑2 (MQ=39)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaat                                                                                    <  1:18405/91‑2 (MQ=39)
          cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  2:3784655/91‑4 (MQ=38)
          cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  2:1892960/91‑4 (MQ=38)
          cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  2:1749136/91‑4 (MQ=38)
            aaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACg                                                                                >  2:3153235/1‑91 (MQ=255)
             aaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGt                                                                               >  1:3856435/1‑91 (MQ=255)
             aaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGt                                                                               >  1:1754546/1‑91 (MQ=255)
               ggTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCa                                                                             >  2:2828309/1‑91 (MQ=255)
               ggTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCa                                                                             >  2:2241199/1‑91 (MQ=255)
               ggTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCa                                                                             >  2:2125397/1‑91 (MQ=255)
                  gacgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATcc                                                                          >  2:1545437/1‑91 (MQ=255)
                    cgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAAAGTCATCCtt                                                                        >  2:3559312/1‑91 (MQ=255)
                     gacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCttt                                                                       >  2:1524747/1‑91 (MQ=255)
                      acTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAGTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTa                                                                      >  2:3002883/1‑91 (MQ=255)
                      acTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTa                                                                      >  2:1023495/1‑91 (MQ=255)
                      acTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTa                                                                      <  1:34048/91‑1 (MQ=255)
                       cTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAt                                                                     >  2:1197800/1‑91 (MQ=255)
                        tCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtt                                                                    >  2:1675568/1‑91 (MQ=255)
                         cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                   >  2:234126/1‑91 (MQ=255)
                         cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                   <  2:548151/91‑1 (MQ=255)
                          ggTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTTTAATTAATCACGTCATCCTTTAttta                                                                  <  2:1114343/91‑2 (MQ=21)
                            aaaCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttttt                                                                <  2:67076/90‑1 (MQ=255)
                              aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                              <  2:1306967/91‑1 (MQ=255)
                              aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                              >  2:1422614/1‑91 (MQ=255)
                              aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                              >  2:3724862/1‑91 (MQ=255)
                               cGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGc                                                             <  2:2233809/91‑1 (MQ=255)
                                 cAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaa                                                           >  2:1568647/1‑91 (MQ=255)
                                  aattaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaa                                                          <  1:2566207/91‑1 (MQ=255)
                                   attaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaa                                                         >  1:111598/1‑91 (MQ=255)
                                   attaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaa                                                         >  2:887060/1‑91 (MQ=255)
                                    ttaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaaa                                                        <  1:2199768/91‑1 (MQ=255)
                                      aattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaata                                                      <  1:2029081/91‑1 (MQ=255)
                                      aattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaata                                                      <  2:857754/91‑1 (MQ=255)
                                       attTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaataa                                                     >  2:2185196/1‑91 (MQ=255)
                                         tttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAAtt                                                   >  1:1220212/1‑91 (MQ=255)
                                         tttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAAtt                                                   >  1:2885187/1‑91 (MQ=255)
                                          ttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTc                                                  >  2:3864246/1‑91 (MQ=255)
                                          ttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTc                                                  <  1:2508735/91‑1 (MQ=255)
                                             ccAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAg                                               <  2:1659919/91‑1 (MQ=255)
                                             ccAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAg                                               >  2:1117063/1‑91 (MQ=255)
                                             ccAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAg                                               <  2:921760/91‑1 (MQ=255)
                                              cAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAgg                                              <  2:2170461/91‑1 (MQ=255)
                                                 tCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGc                                           >  1:1195786/1‑91 (MQ=255)
                                                  cAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCg                                          <  2:2687042/91‑1 (MQ=255)
                                                  cAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCg                                          <  1:2615323/91‑1 (MQ=255)
                                                   aGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGt                                         >  1:2067362/1‑91 (MQ=255)
                                                     aaCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                                >  2:3817704/1‑82 (MQ=255)
                                                     aaCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                                <  1:3817704/82‑1 (MQ=255)
                                                      aCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATg                                      <  1:878075/91‑1 (MQ=255)
                                                         tACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                   >  1:2692300/1‑91 (MQ=255)
                                                           cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                 <  1:1568593/91‑1 (MQ=255)
                                                           cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                 <  2:2398306/91‑1 (MQ=255)
                                                              aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                              >  1:434629/1‑91 (MQ=255)
                                                              aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                              >  1:3186639/1‑91 (MQ=255)
                                                              aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                              >  1:2369707/1‑91 (MQ=255)
                                                               aaataaatATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataaa                             >  1:1335029/1‑91 (MQ=255)
                                                               aaataaatATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATATTTTGGCAAAAATAATTCTAGAAGCGAAAGTAaaataaa                             >  2:3376461/1‑91 (MQ=255)
                                                                aataaataTAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataaat                            >  1:1443669/1‑91 (MQ=255)
                                                                 ataaatataaatAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAtt                           >  1:61621/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:344949/91‑1 (MQ=255)
                                                                    aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                        <  1:2423620/91‑1 (MQ=255)
                                                                     atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                   >  1:2864121/1‑79 (MQ=255)
                                                                     atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                   <  2:2864121/79‑1 (MQ=255)
                                                                     atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgc                       >  1:1981747/1‑91 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     >  1:1973596/1‑91 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     >  2:3636107/1‑91 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATACTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     >  2:169689/1‑91 (MQ=255)
                                                                          aataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAt                  <  1:666641/91‑1 (MQ=255)
                                                                             aaaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAtttt               >  1:2263683/1‑91 (MQ=21)
                                                                              aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa              >  1:2817271/1‑91 (MQ=25)
                                                                               aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa             >  2:272396/1‑91 (MQ=255)
                                                                               aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa             <  2:909137/91‑1 (MQ=255)
                                                                               aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa             <  1:2222116/91‑1 (MQ=255)
                                                                                   tGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAg         <  1:848609/91‑1 (MQ=255)
                                                                                    gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc        >  1:2720441/1‑91 (MQ=255)
                                                                                    gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc        >  2:1162925/1‑91 (MQ=255)
                                                                                    gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc        >  1:1125608/1‑91 (MQ=255)
                                                                                    gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc        >  2:1163174/1‑91 (MQ=255)
                                                                                     aCGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgca                      <  2:3751190/76‑1 (MQ=255)
                                                                                     aCGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgca                      >  1:3751190/1‑76 (MQ=255)
                                                                                      cGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCaa      >  2:457899/1‑91 (MQ=255)
                                                                                        tGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAAt    >  2:3698802/1‑91 (MQ=255)
                                                                                        tGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAAt    <  1:39576/91‑1 (MQ=255)
                                                                                         gATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATa   >  1:3067461/1‑91 (MQ=255)
                                                                                         gATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATa   <  1:3698802/91‑1 (MQ=255)
                                                                                          attaattaAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATaa  >  1:100595/1‑91 (MQ=255)
                                                                                          attaattaAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATaa  >  2:880536/1‑91 (MQ=255)
                                                                                          attaattaAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATaa  >  2:22061/1‑91 (MQ=255)
                                                                                          |                                                                                           
ATTCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAA  >  USA300TCH1516_ALE/3254‑3434

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