Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 3,344 A→T 15.6% 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,3440AT15.6% 444.5 / 64.9 153intergenic (+28/‑362)dnaN/USA300HOU_RS00015DNA polymerase III subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (72/57);  new base T (9/15);  total (81/72)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 158 total aligned reads displayed.

TCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAA  >  USA300TCH1516_ALE/3256‑3434
                                                                                        |                                                                                           
tCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGatt                                                                                           <  1:765009/91‑1 (MQ=255)
tCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAAGAAAGGATGACGTGatt                                                                                           <  1:4045334/91‑1 (MQ=255)
     aaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAACTATAAATAAAGGATGACGTTttaattaa                                                                                     <  1:2564556/91‑1 (MQ=38)
      aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCATTCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGAttaattaa                                                                                     <  1:3893561/91‑1 (MQ=255)
      aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGAttaattaa                                                                                     <  1:3887272/91‑1 (MQ=255)
      aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaat                                                                                    <  1:718921/91‑2 (MQ=39)
       ccAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaa                                                                                   >  1:2948540/1‑91 (MQ=255)
       ccAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaa                                                                                   >  1:3472283/1‑91 (MQ=255)
       ccAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaa                                                                                   >  1:4090141/1‑91 (MQ=255)
        cAATAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  1:3339518/91‑4 (MQ=38)
        cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaaa                                                                                  <  1:3778989/91‑1 (MQ=255)
        cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaaa                                                                                  <  1:1110244/91‑1 (MQ=255)
        cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  1:2859233/91‑4 (MQ=38)
        cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  1:2858243/91‑4 (MQ=38)
        cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  1:1787021/91‑4 (MQ=38)
        cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  1:3562972/91‑4 (MQ=38)
        cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  1:4015243/91‑4 (MQ=38)
        cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  1:3818300/91‑4 (MQ=38)
         aaaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAAc                                                                                 <  1:2163419/91‑1 (MQ=255)
          aaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAATTAAAGGATGACGTGATTAATTAAAACg                                                                                <  1:810026/91‑1 (MQ=255)
          aaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACg                                                                                >  1:1437897/1‑91 (MQ=255)
             ggTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCa                                                                             >  1:1384952/1‑91 (MQ=255)
             ggTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCa                                                                             >  1:1384951/1‑91 (MQ=255)
                gacgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAGCGTCATcc                                                                          <  1:1489727/91‑1 (MQ=255)
                 acgacTCGGTAACGCAATTAATTTTGCCAATCAGAACTTGCTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         <  1:3850384/91‑1 (MQ=255)
                 acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         >  1:2829013/1‑91 (MQ=255)
                 acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         >  1:2614677/1‑91 (MQ=255)
                 acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         >  1:2614354/1‑91 (MQ=255)
                 acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         >  1:1860699/1‑91 (MQ=255)
                 acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         >  1:1860720/1‑91 (MQ=255)
                 acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         <  1:2429279/91‑1 (MQ=255)
                 acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         <  1:2427505/91‑1 (MQ=255)
                   gacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCttt                                                                       >  1:2305451/1‑91 (MQ=255)
                   gacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCttt                                                                       <  1:1952003/91‑1 (MQ=255)
                       cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                   >  1:550069/1‑91 (MQ=255)
                       cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                   >  1:549793/1‑91 (MQ=255)
                       cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                   >  1:605212/1‑91 (MQ=255)
                       cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                   <  1:2109473/91‑1 (MQ=255)
                       cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTTTAATTAATCACGTCATCCTTTAttt                                                                   >  1:1137036/1‑91 (MQ=21)
                        ggTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttt                                                                  >  1:1602184/1‑91 (MQ=255)
                         gTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttttt                                                                 >  1:4115934/1‑91 (MQ=255)
                         gTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttttt                                                                 <  1:1572821/91‑1 (MQ=255)
                           aaCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTTTAATTAATCACGTCATCCTTTAtttatat                                                               <  1:449985/91‑5 (MQ=18)
                            aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                              >  1:2415074/1‑91 (MQ=255)
                            aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                              >  1:1730536/1‑91 (MQ=255)
                             cGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGc                                                             >  1:708768/1‑91 (MQ=255)
                              gCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCa                                                            <  1:3098405/91‑1 (MQ=255)
                              gCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCa                                                            <  1:2098853/91‑1 (MQ=255)
                              gCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCa                                                            >  1:2941065/1‑91 (MQ=255)
                                aattaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaa                                                          <  1:2130931/91‑1 (MQ=255)
                                aattaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaa                                                          <  1:2130949/91‑1 (MQ=255)
                                 attaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaa                                                         >  1:275151/1‑91 (MQ=255)
                                    aattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaata                                                      <  1:1635856/91‑1 (MQ=255)
                                    aattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaata                                                      <  1:3917221/91‑1 (MQ=255)
                                      ttttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaataat                                                    <  1:3989017/91‑1 (MQ=255)
                                        ttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTc                                                  <  1:2800055/91‑1 (MQ=255)
                                          aCCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                                >  1:221879/1‑91 (MQ=255)
                                             aaTCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGt                                             >  1:1316094/1‑91 (MQ=255)
                                              aTCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTg                                            <  1:1750470/91‑1 (MQ=255)
                                               tCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGc                                           >  1:1721488/1‑91 (MQ=255)
                                                cAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCg                                          <  1:2456794/91‑1 (MQ=255)
                                                cAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGATGCg                                          >  1:3458779/1‑91 (MQ=255)
                                                 aGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGt                                         >  1:2868891/1‑91 (MQ=255)
                                                 aGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGt                                         >  1:2253993/1‑91 (MQ=255)
                                                       tACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                   <  1:1763678/91‑1 (MQ=255)
                                                       tACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                   >  1:3013335/1‑91 (MQ=255)
                                                        aCTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaa                                  <  1:1263316/91‑1 (MQ=255)
                                                         cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                 >  1:1568503/1‑91 (MQ=255)
                                                         cTAAAAATAAATATAAATAAAGGATGACGTGATGAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                 <  1:3392231/91‑1 (MQ=255)
                                                          tAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaat                                >  1:2377861/1‑91 (MQ=255)
                                                           aaaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCGTCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaata                               <  1:3293953/91‑1 (MQ=255)
                                                            aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                              <  1:80753/91‑1 (MQ=255)
                                                            aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                              >  1:920985/1‑91 (MQ=255)
                                                                taaatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAttt                          <  1:3510915/91‑1 (MQ=255)
                                                                  aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                        <  1:1680079/91‑1 (MQ=255)
                                                                  aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                        >  1:1690092/1‑91 (MQ=255)
                                                                   atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgc                       >  1:3419621/1‑91 (MQ=255)
                                                                   atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgc                       >  1:609094/1‑91 (MQ=255)
                                                                    tataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgca                      <  1:3078348/91‑1 (MQ=255)
                                                                     atacataaATGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     <  1:3733680/91‑1 (MQ=255)
                                                                     ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     >  1:104724/1‑91 (MQ=255)
                                                                        aataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAt                  >  1:451061/1‑91 (MQ=255)
                                                                          taaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAttt                <  1:2477736/91‑1 (MQ=255)
                                                                           aaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAtttt               >  1:1364210/1‑91 (MQ=255)
                                                                           aaaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAtttt               >  1:3254077/1‑91 (MQ=21)
                                                                           aaaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAtttt               >  1:5849/1‑91 (MQ=21)
                                                                            aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa              >  1:606578/1‑91 (MQ=25)
                                                                            aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa              >  1:2922726/1‑91 (MQ=25)
                                                                            aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa              >  1:642244/1‑91 (MQ=25)
                                                                            aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa              >  1:3976106/1‑91 (MQ=25)
                                                                             aGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaa             >  1:2295322/1‑91 (MQ=255)
                                                                              ggATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaaa            >  1:2165205/1‑91 (MQ=255)
                                                                              ggATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaaa            >  1:1058779/1‑91 (MQ=255)
                                                                                 tGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAg         >  1:1534476/1‑91 (MQ=255)
                                                                                 tGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAg         <  1:885518/91‑1 (MQ=255)
                                                                                 tGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAg         <  1:885609/91‑1 (MQ=255)
                                                                                   aCGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCa       <  1:2234330/91‑1 (MQ=255)
                                                                                   aCGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCa       <  1:3508703/91‑1 (MQ=255)
                                                                                     gTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCaaa     >  1:2294862/1‑91 (MQ=255)
                                                                                        attaattaAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATaa  >  1:3963397/1‑91 (MQ=255)
                                                                                        |                                                                                           
TCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAA  >  USA300TCH1516_ALE/3256‑3434

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