Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 3,343 Δ1 bp 7.6% intergenic (+27/‑363) 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,3430G. 7.6% 495.8 / 33.6 144intergenic (+27/‑363)dnaN/USA300HOU_RS00015DNA polymerase III subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (73/60);  new base . (9/2);  total (82/62)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 158 total aligned reads displayed.

TATTCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAA  >  USA300TCH1516_ALE/3253‑3431
                                                                                          |                                                                                         
   tCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGatt                                                                                        <  1:765009/91‑1 (MQ=255)
   tCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAAGAAAGGATGACGTGatt                                                                                        <  1:4045334/91‑1 (MQ=255)
        aaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAACTATAAATAAAGGATGACGTTttaattaa                                                                                  <  1:2564556/91‑1 (MQ=38)
        aaCCAAAAGGTGACGACTCGGTAACACAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaa                                                                                  <  1:3876603/91‑1 (MQ=39)
         aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGAttaattaa                                                                                  <  1:3887272/91‑1 (MQ=255)
         aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaat                                                                                 <  1:718921/91‑2 (MQ=39)
          ccAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaa                                                                                >  1:3472283/1‑91 (MQ=255)
          ccAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaa                                                                                >  1:2948540/1‑91 (MQ=255)
          ccAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaa                                                                                >  1:4090141/1‑91 (MQ=255)
           cAATAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                               <  1:3339518/91‑4 (MQ=38)
           cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaaa                                                                               <  1:1110244/91‑1 (MQ=255)
           cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                               <  1:1787021/91‑4 (MQ=38)
           cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                               <  1:2858243/91‑4 (MQ=38)
           cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                               <  1:2859233/91‑4 (MQ=38)
           cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                               <  1:4015243/91‑4 (MQ=38)
           cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                               <  1:3565188/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)
              aaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGt                                                                            >  1:3915531/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)
                    acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                      >  1:2614354/1‑91 (MQ=255)
                    acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                      <  1:2429279/91‑1 (MQ=255)
                    acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                      <  1:2427505/91‑1 (MQ=255)
                    acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                      >  1:2614677/1‑91 (MQ=255)
                    acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                      >  1:2829013/1‑91 (MQ=255)
                    acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                      >  1:1860720/1‑91 (MQ=255)
                    acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                      >  1:1860699/1‑91 (MQ=255)
                      gacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCttt                                                                    >  1:2305451/1‑91 (MQ=255)
                      gacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCttt                                                                    <  1:1952003/91‑1 (MQ=255)
                          cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                <  1:2109473/91‑1 (MQ=255)
                          cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                >  1:605212/1‑91 (MQ=255)
                          cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                >  1:550069/1‑91 (MQ=255)
                          cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                >  1:549793/1‑91 (MQ=255)
                          cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTTTAATTAATCACGTCATCCTTTAttt                                                                >  1:1137036/1‑91 (MQ=21)
                           ggTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttt                                                               >  1:1602184/1‑91 (MQ=255)
                            gTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttttt                                                              <  1:1572821/91‑1 (MQ=255)
                            gTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttttt                                                              >  1:4115934/1‑91 (MQ=255)
                             tAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttttt                                                             >  1:3529455/1‑91 (MQ=255)
                             tAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttttt                                                             >  1:3794739/1‑91 (MQ=255)
                              aaCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTTTAATTAATCACGTCATCCTTTAtttatat                                                            <  1:449985/91‑5 (MQ=18)
                               aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                           >  1:1730536/1‑91 (MQ=255)
                               aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                           >  1:2415074/1‑91 (MQ=255)
                                cGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGc                                                          >  1:708768/1‑91 (MQ=255)
                                 gCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCa                                                         >  1:2941065/1‑91 (MQ=255)
                                 gCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCa                                                         <  1:2098853/91‑1 (MQ=255)
                                 gCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCa                                                         <  1:3098405/91‑1 (MQ=255)
                                   aattaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaa                                                       <  1:2130949/91‑1 (MQ=255)
                                   aattaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaa                                                       <  1:2130931/91‑1 (MQ=255)
                                    attaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaa                                                      >  1:275151/1‑91 (MQ=255)
                                       aattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaata                                                   <  1:1635856/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:3013335/1‑91 (MQ=255)
                                                          tACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                <  1:1763678/91‑1 (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:3956931/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)
                                                                                          |                                                                                         
TATTCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAA  >  USA300TCH1516_ALE/3253‑3431

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: