Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 3,344 A→T 9.3% 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,3440AT 9.3% 830.6 / 56.0 256intergenic (+28/‑362)dnaN/USA300HOU_RS00015DNA polymerase III subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (137/94);  new base T (7/17);  total (145/111)
Fisher's exact test for biased strand distribution p-value = 8.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 278 total aligned reads displayed.

ATTCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAA  >  USA300TCH1516_ALE/3254‑3434
                                                                                          |                                                                                           
    tAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaa                                                                                         >  2:3606707/1‑91 (MQ=255)
     aaaaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaat                                                                                        >  1:1136518/1‑91 (MQ=255)
      aaaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaatt                                                                                       >  1:153937/1‑91 (MQ=255)
       aaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaatta                                                                                      <  1:553383/91‑1 (MQ=255)
       aaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaatta                                                                                      <  1:378016/91‑1 (MQ=255)
       aaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaatta                                                                                      <  1:2659035/91‑1 (MQ=255)
       aaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaa                                                                                     <  1:149435/91‑1 (MQ=255)
       aaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaa                                                                                     <  1:3549503/91‑1 (MQ=255)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGAttaattaa                                                                                     <  2:2416139/91‑1 (MQ=255)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGAttaattaa                                                                                     >  2:1126955/1‑91 (MQ=255)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGAttaattaa                                                                                     <  1:1740965/91‑1 (MQ=255)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaat                                                                                    <  2:984691/91‑2 (MQ=39)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaat                                                                                    <  2:454369/91‑2 (MQ=39)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaat                                                                                    <  2:1674786/91‑2 (MQ=39)
        aCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaat                                                                                    <  1:1826799/91‑2 (MQ=39)
         ccAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaa                                                                                   >  1:3699940/1‑91 (MQ=255)
          cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaaa                                                                                  >  1:1220612/1‑91 (MQ=255)
          cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaaa                                                                                  <  2:3439357/91‑1 (MQ=255)
          cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTttaattaatca                                                                                  <  1:3523534/91‑4 (MQ=38)
             aaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGt                                                                               >  1:1309222/1‑91 (MQ=255)
             aaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGt                                                                               >  1:3116877/1‑91 (MQ=255)
             aaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGt                                                                               >  2:3718448/1‑91 (MQ=255)
              aGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAGAAATAAATATAAATAAAGGATGTCGTGATTAATTAAAACGTc                                                                              >  2:1153501/1‑91 (MQ=255)
              aGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTc                                                                              >  1:3289891/1‑91 (MQ=255)
                gTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCAt                                                                            >  1:3883236/1‑91 (MQ=255)
                gTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCAt                                                                            >  2:3151762/1‑91 (MQ=255)
                 tGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTTTTAATTAATCACGTCATc                                                                           <  2:2036366/91‑1 (MQ=25)
                  gacgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATcc                                                                          <  2:679720/91‑1 (MQ=255)
                  gacgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATcc                                                                          >  2:2191873/1‑91 (MQ=255)
                  gacgacTCGGTAACGCAATGACTTTTGCCAATCAGAACGTACTAACAACAACTATAAACAAAGCACTACGGGCTTCATTAAACCGTCATcc                                                                          >  2:1591754/1‑91 (MQ=255)
                   acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         <  1:781704/91‑1 (MQ=255)
                   acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         <  2:3076888/91‑1 (MQ=255)
                   acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         >  1:2167107/1‑91 (MQ=255)
                   acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         <  2:1076923/91‑1 (MQ=255)
                    cgacTCGGTAACGCAATTAATTTTACCAATCAGGACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                        <  2:1235466/91‑1 (MQ=255)
                    cgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAGCGTCATCCtt                                                                        <  1:2356937/91‑1 (MQ=255)
                    cgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                        >  2:3098097/1‑91 (MQ=255)
                    cgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                        <  2:2820539/91‑1 (MQ=255)
                     gacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCttt                                                                       >  2:3585152/1‑91 (MQ=255)
                      acTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTa                                                                      >  1:2271461/1‑91 (MQ=255)
                      acTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTa                                                                      <  1:2853645/91‑1 (MQ=255)
                       cTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAt                                                                     <  1:3533133/91‑1 (MQ=255)
                          ggTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttt                                                                  >  1:297188/1‑91 (MQ=255)
                           gTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttttt                                                                 <  2:1365824/91‑1 (MQ=255)
                              aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                              >  2:2766551/1‑91 (MQ=255)
                                gCAATTGAGTTTACCAGTCAGGACTTACTAAAGATAAATATAGATAGGGGATGACGTGATTAATTAAACCGTCATCCTGTATTTTTTGGCa                                                            <  2:2903754/91‑1 (MQ=255)
                                gCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCa                                                            <  1:2874410/91‑1 (MQ=255)
                                   attaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaa                                                         >  1:1859561/1‑91 (MQ=255)
                                    ttaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaaa                                                        >  2:2934146/1‑91 (MQ=255)
                                        ttttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaataat                                                    <  1:2507904/91‑1 (MQ=255)
                                          ttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTc                                                  >  2:3386287/1‑91 (MQ=255)
                                          ttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTc                                                  <  1:202653/91‑1 (MQ=255)
                                            aCCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                                <  2:1930326/91‑1 (MQ=255)
                                            aCCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                                <  1:2691296/91‑1 (MQ=255)
                                             ccAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAg                                               >  1:3425058/1‑91 (MQ=255)
                                                aTCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCGTCCTTTATTTTTTGGCAAAAATAATTCTAGGTg                                            >  1:3240146/1‑91 (MQ=255)
                                                aTCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGGGATTAATTAAAACATCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTg                                            <  2:2468227/91‑1 (MQ=255)
                                                  cAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCg                                          <  1:1194244/91‑1 (MQ=255)
                                                    gAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTa                                        >  1:2156520/1‑91 (MQ=255)
                                                    gAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTa                                        >  1:136820/1‑91 (MQ=255)
                                                      aCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATg                                      >  2:1522282/1‑91 (MQ=255)
                                                       cTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTTTTTTTTATTATTCTAGGTGCGTATGt                                     >  2:1832417/1‑91 (MQ=37)
                                                       cTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGt                                     <  2:1402029/91‑1 (MQ=255)
                                                        ttACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTa                                    >  2:3348459/1‑91 (MQ=255)
                                                          aCTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaa                                  >  1:2611301/1‑91 (MQ=255)
                                                           cTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaaaa                                 >  2:60770/1‑91 (MQ=255)
                                                             aaaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaata                               >  2:3817260/1‑91 (MQ=255)
                                                              aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                              >  1:1636034/1‑91 (MQ=255)
                                                                aataaataTAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataaat                            >  1:3066986/1‑91 (MQ=255)
                                                                 ataaatataaatAAAGGATGACGCGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTGGGTGCGTATGTAAAATAAAtt                           >  2:3490084/1‑91 (MQ=255)
                                                                  taaatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAttt                          <  1:3612946/91‑1 (MQ=255)
                                                                   aaatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTg                         >  2:1901747/1‑91 (MQ=255)
                                                                     atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgc                       >  1:888890/1‑91 (MQ=255)
                                                                     atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgc                       >  1:2759843/1‑91 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGtcag                     >  2:3944818/1‑91 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     <  1:2627847/91‑1 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     >  2:255131/1‑91 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     <  1:994465/91‑1 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     >  1:522178/1‑91 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     <  2:10433/91‑1 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     >  1:3461357/1‑91 (MQ=255)
                                                                            taaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAttt                >  2:601438/1‑91 (MQ=255)
                                                                            taaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAttt                <  1:3638686/91‑1 (MQ=255)
                                                                              aaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa              >  2:1119263/1‑91 (MQ=255)
                                                                              aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa              >  2:3149735/1‑91 (MQ=25)
                                                                                ggATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCATAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaaa            >  2:2495541/1‑91 (MQ=255)
                                                                                ggATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTaaa            >  2:391591/1‑91 (MQ=255)
                                                                                ggATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTata            >  2:1552451/1‑89 (MQ=25)
                                                                                 gATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATATTAAAc           >  2:2098170/1‑91 (MQ=255)
                                                                                   tGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAg         >  1:1936277/1‑91 (MQ=255)
                                                                                   tGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAg         >  1:3913575/1‑91 (MQ=255)
                                                                                   tGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAg         >  1:2474474/1‑91 (MQ=255)
                                                                                   ggtcGTGATTAATTAAAACGTCATCCTTTAATTTTTAGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTTGCAGCATTTTAAACAg         <  2:903189/88‑1 (MQ=255)
                                                                                    gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc        >  1:1052704/1‑91 (MQ=255)
                                                                                      cGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCaa      >  1:491627/1‑91 (MQ=255)
                                                                                       gTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCaaa     >  1:826062/1‑91 (MQ=255)
                                                                                       gTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCaaa     >  2:3243781/1‑91 (MQ=255)
                                                                                       gTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCaaa     >  1:3837244/1‑91 (MQ=255)
                                                                                       gTGATTAATTAAAACGTCATCCTTTATTTGTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCaaa     <  2:1289207/91‑1 (MQ=255)
                                                                                          attaattaAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATaa  >  1:161210/1‑91 (MQ=255)
                                                                                          |                                                                                           
ATTCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAA  >  USA300TCH1516_ALE/3254‑3434

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