Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 3,344 A→T 5.2% 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 5.2% 1377.3 / 48.6 376intergenic (+28/‑362)dnaN/USA300HOU_RS00015DNA polymerase III subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (172/184);  new base T (4/16);  total (176/200)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 405 total aligned reads displayed.

ATTCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAATAA  >  USA300TCH1516_ALE/3254‑3434
                                                                                          |                                                                                           
 ttCTAAAACCAAATGGTGACGACTCGTTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGat                                                                                            <  2:3763838/91‑1 (MQ=255)
 ttCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGat                                                                                            <  2:5497089/91‑1 (MQ=255)
 ttCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGat                                                                                            <  2:2710175/91‑1 (MQ=255)
  tCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGatt                                                                                           >  1:1090064/1‑91 (MQ=255)
  tCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGatt                                                                                           <  1:228078/91‑1 (MQ=255)
  tCTAAAACCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGatt                                                                                           <  2:1133708/91‑1 (MQ=255)
      aaaCCAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGattaatt                                                                                       >  2:3872398/1‑91 (MQ=255)
         ccAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaa                                                                                   >  1:46605/1‑91 (MQ=255)
          cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaaa                                                                                  <  2:1514678/91‑1 (MQ=255)
          cAAAAGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTaaaa                                                                                  <  1:931960/91‑1 (MQ=255)
           aaaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAAc                                                                                 >  1:3278033/1‑91 (MQ=255)
            aaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACg                                                                                >  1:5274996/1‑91 (MQ=255)
            aaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACg                                                                                >  2:2467956/1‑91 (MQ=255)
            aaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACg                                                                                >  1:1952907/1‑91 (MQ=255)
            aaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACg                                                                                >  2:4094874/1‑91 (MQ=255)
            aaaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACg                                                                                <  1:2574740/91‑1 (MQ=255)
             aaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGCACTTACTAAAAATAAATATAAATAAAGGATGACGTGACTAATTAAAACGt                                                                               >  2:3602999/1‑91 (MQ=255)
             aaGGTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGt                                                                               >  2:206179/1‑91 (MQ=255)
                gTGACGACTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCAt                                                                            >  1:3777675/1‑91 (MQ=255)
                 tGACGACTCGGTAACGCAATGAATTTTACCAGTCAGAACTTACTAAAAATAAATATATATAAAGGATGACGTGATTAATTGAAACGTCATc                                                                           >  1:5158870/1‑91 (MQ=255)
                  gacgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATcc                                                                          <  2:238386/91‑1 (MQ=255)
                  gacgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATcc                                                                          <  1:5641616/91‑1 (MQ=255)
                  gacgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATcc                                                                          >  2:5644990/1‑91 (MQ=255)
                   acgacTCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         <  1:3598218/91‑1 (MQ=255)
                   acgacCCGGTAACGCAATTAATTTTACCAATCAGAAATTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCt                                                                         <  2:3469838/91‑1 (MQ=255)
                    cgacTCGGTAACGCAATTAATTTTACCAATCATAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCtt                                                                        <  1:5533791/91‑1 (MQ=255)
                        tCGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtt                                                                    <  2:4244629/91‑1 (MQ=255)
                         cGGTAACGCAATTAATTTTACCAATCAGAACTTACTACAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                   >  1:4313800/1‑91 (MQ=255)
                         cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                   >  2:138775/1‑91 (MQ=255)
                         cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                   >  2:120602/1‑91 (MQ=255)
                         cGGTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttt                                                                   <  2:2217435/91‑1 (MQ=255)
                          ggTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttt                                                                  >  2:3123132/1‑91 (MQ=255)
                          ggTAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttt                                                                  <  1:5068029/91‑1 (MQ=255)
                           gTAACGAAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAttttt                                                                 <  2:299421/91‑1 (MQ=255)
                            tAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttttt                                                                >  2:1802333/1‑91 (MQ=255)
                            tAACGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTAtttttt                                                                >  2:1690685/1‑91 (MQ=255)
                              aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                              >  1:2038647/1‑91 (MQ=255)
                              aCGCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTgg                                                              <  1:723651/91‑1 (MQ=255)
                                gCAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCa                                                            >  1:1760447/1‑91 (MQ=255)
                                 cAATTAATTTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaa                                                           >  1:1151266/1‑91 (MQ=255)
                                  aattaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaa                                                          <  2:3488467/91‑1 (MQ=255)
                                   attaattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCaaaa                                                         >  1:2541739/1‑91 (MQ=255)
                                     taattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaat                                                       <  1:1006048/91‑1 (MQ=255)
                                      aattTTACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaata                                                      >  1:2057208/1‑91 (MQ=255)
                                        ttttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaataat                                                    <  1:4093256/91‑1 (MQ=255)
                                        ttttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaataat                                                    <  2:2662393/91‑1 (MQ=255)
                                        ttttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAaataat                                                    <  1:4024823/91‑1 (MQ=255)
                                         tttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAAtt                                                   >  2:3829462/1‑91 (MQ=255)
                                          ttgccAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTc                                                  <  2:4994188/88‑1 (MQ=255)
                                          ttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTc                                                  <  2:1556242/91‑1 (MQ=255)
                                          ttACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGAAAAAATAATTc                                                  <  2:5203238/91‑1 (MQ=255)
                                           tACCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTCTGGCAAAAATAATTCt                                                 >  2:5373116/1‑91 (MQ=255)
                                            aCCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                                >  1:3368383/1‑91 (MQ=255)
                                            aCCAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTa                                                >  1:1217573/1‑91 (MQ=255)
                                              cAATCAGAACTTACTAAAAATAAATATAAATAAAGGATGACATGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAgg                                              >  2:3802620/1‑91 (MQ=255)
                                               aaTCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGt                                             <  2:3325806/91‑1 (MQ=255)
                                               aaTCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGt                                             >  1:3203996/1‑91 (MQ=255)
                                                aTCAGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTg                                            <  2:210895/91‑1 (MQ=255)
                                                   aGAACTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGt                                         >  1:1541621/1‑91 (MQ=255)
                                                      aCTTACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATg                                      <  1:2159716/91‑1 (MQ=255)
                                                        ttACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTa                                    <  1:1382438/91‑1 (MQ=255)
                                                         tACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                   <  1:5373712/91‑1 (MQ=255)
                                                         tACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                   >  1:4953156/1‑91 (MQ=255)
                                                         tACTAAAAATAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTaa                                   <  1:1335205/91‑1 (MQ=255)
                                                             aaaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaata                               <  1:3039185/91‑1 (MQ=255)
                                                             aaaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaata                               <  2:4903497/91‑1 (MQ=255)
                                                              aaaaTAAATATAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataa                              >  2:1979043/1‑91 (MQ=255)
                                                                aataaataTAAATAAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAaaataaat                            >  2:1437785/1‑91 (MQ=255)
                                                                 ataaatataaatAAAGGATGACGTGATTAATTAAAACGTCATCCTGTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAtt                           <  2:5458932/91‑1 (MQ=255)
                                                                  taaatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAttt                          <  2:482748/91‑1 (MQ=255)
                                                                  taaatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAAttt                          <  2:41770/91‑1 (MQ=255)
                                                                    aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                        <  2:709996/91‑1 (MQ=255)
                                                                    aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                        >  2:281149/1‑91 (MQ=255)
                                                                    aatataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTgg                        <  2:4542045/91‑1 (MQ=255)
                                                                     atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAGTTCTAGGTGCGTATGTAAAATAAATTTGgc                       <  1:4500930/91‑1 (MQ=255)
                                                                     atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgc                       >  1:1449046/1‑91 (MQ=255)
                                                                     atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgc                       >  1:4915629/1‑91 (MQ=255)
                                                                     atataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgc                       <  1:3174986/91‑1 (MQ=255)
                                                                      tataaataAAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgca                      >  1:2720481/1‑91 (MQ=255)
                                                                       ataaataaAGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     >  2:5309/1‑91 (MQ=255)
                                                                       ataaataaAGAATAACGTGATTAATTAAAACGTAATCTTTAATAAAAAAAAAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcag                     <  2:1305720/91‑1 (MQ=255)
                                                                        taaataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcagc                    <  2:2538739/91‑1 (MQ=255)
                                                                        taaataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcagc                    <  2:4613455/91‑1 (MQ=255)
                                                                        taaataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGgcagc                    <  1:2938281/91‑1 (MQ=255)
                                                                          aataaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAt                  <  1:832050/91‑1 (MQ=255)
                                                                            taaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAttt                >  2:4710318/1‑91 (MQ=255)
                                                                            taaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAttt                >  2:9060/1‑91 (MQ=255)
                                                                             aaaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCAtttt               >  1:4238073/1‑91 (MQ=255)
                                                                              aaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa              <  2:1889195/91‑1 (MQ=255)
                                                                              aaGGATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa              <  1:3403897/91‑1 (MQ=255)
                                                                              aaGGATGACGTTTTAATTAATCACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTa              >  1:1611127/1‑91 (MQ=25)
                                                                                ggATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTTTTTTTTTTTATTATAGGAGCGTATGTAAAATAAATTTTGCAGCATTTTaaa            >  2:944325/1‑91 (MQ=255)
                                                                                ggATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTAGGCAGCATTTTaaa            >  1:1260167/1‑91 (MQ=255)
                                                                                ggATGACGTGATTAATTAAAACGGCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCACCATTTTaca            >  2:3984022/1‑89 (MQ=255)
                                                                                 gATGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAAc           <  2:4356648/91‑1 (MQ=255)
                                                                                  aTGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACa          <  1:612008/91‑1 (MQ=255)
                                                                                   tGACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAg         <  1:1829802/91‑1 (MQ=255)
                                                                                    gACGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGc        >  2:3914787/1‑91 (MQ=255)
                                                                                      cGTGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCaa      >  1:280202/1‑91 (MQ=255)
                                                                                        tGATTAATTAAAACGTCATCCTTTATTTTTTGGCAAAAATAATTCTAGGTGCGTATGTAAAATAAATTTGGCAGCATTTTAAACAGCAAAt    <  1:2383399/91‑1 (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: