Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 79,619 T→A 11.2% intergenic (+228/‑32) USA300HOU_RS00355 → / → cmoB hypothetical protein/tRNA U34 carboxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE79,6190TA11.2% 641.6 / 59.6 215intergenic (+228/‑32)USA300HOU_RS00355/cmoBhypothetical protein/tRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (90/100);  new base A (13/11);  total (103/112)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

Only 100 of 242 total aligned reads displayed.

GACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGC  >  USA300TCH1516_ALE/79529‑79709
                                                                                          |                                                                                          
gACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAatt                                                                                            <  1:779644/91‑1 (MQ=255)
 aCCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAatta                                                                                           <  1:2451090/91‑1 (MQ=255)
 aCCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAatta                                                                                           <  2:3331087/91‑1 (MQ=255)
 aCCACAGCATAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTAGATAGTGTAGTACTATTCAAAACGTAatta                                                                                           <  2:3373559/91‑1 (MQ=255)
  ccACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattat                                                                                          <  1:1388705/91‑1 (MQ=255)
   cacaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattatt                                                                                         <  2:2273891/91‑1 (MQ=255)
    acaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAAtaatta                                                                                        <  2:669529/91‑1 (MQ=255)
    acaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAAtaatta                                                                                        <  2:574487/91‑1 (MQ=255)
     caGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTCTTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTAc                                                                                       >  2:2993294/1‑91 (MQ=255)
     caGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTAc                                                                                       >  2:3548324/1‑91 (MQ=255)
     caGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTAc                                                                                       <  1:2053032/91‑1 (MQ=255)
       gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACgt                                                                                     <  2:3530230/91‑2 (MQ=255)
       gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTACTTATTACGa                                                                                     >  1:3719761/1‑91 (MQ=255)
       gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                                     <  1:3214195/91‑1 (MQ=255)
        cTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAt                                                                                    <  1:1918724/91‑1 (MQ=255)
         ttAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAtt                                                                                   <  1:1547243/91‑1 (MQ=255)
         ttAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTAATATTCAAAACGTAATTATTACGAtt                                                                                   >  1:3615031/1‑91 (MQ=255)
          tAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAttt                                                                                  <  2:344465/91‑1 (MQ=255)
           aaTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTg                                                                                 >  2:2961614/1‑91 (MQ=255)
             tttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGaa                                                                               <  1:835400/91‑1 (MQ=38)
             tttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGaa                                                                               <  2:1369301/91‑1 (MQ=255)
              ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAt                                                                              <  1:3247488/91‑1 (MQ=38)
              ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAt                                                                              >  2:1900015/1‑91 (MQ=255)
              ttttGGTGTTCATTAATGAAAATAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAt                                                                              <  2:2948486/91‑1 (MQ=255)
               tttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAta                                                                             <  2:2525708/91‑2 (MQ=38)
                 tGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTaa                                                                           >  1:2313726/1‑91 (MQ=255)
                 tGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTaa                                                                           >  2:2906916/1‑91 (MQ=255)
                    tgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCg                                                                        >  2:1266854/1‑91 (MQ=255)
                     gtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCga                                                                       >  1:3349696/1‑91 (MQ=255)
                      ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                                      >  1:3318879/1‑91 (MQ=255)
                      ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                                      <  2:362270/91‑1 (MQ=255)
                       tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATAAAGCgagg                                                                     <  2:507449/91‑1 (MQ=255)
                        cattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagga                                                                    <  2:618589/91‑1 (MQ=255)
                          ttcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGgaga                                                                  >  1:1888088/1‑91 (MQ=255)
                           tcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaa                                                                 <  2:2707728/91‑1 (MQ=255)
                               gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                             >  2:2662574/1‑91 (MQ=255)
                               gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                             >  2:1129561/1‑91 (MQ=255)
                               gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                             >  1:1664991/1‑91 (MQ=255)
                                aCTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaa                                                            <  2:2522005/91‑1 (MQ=255)
                                 cTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaat                                                           <  2:784949/91‑1 (MQ=255)
                                   ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                         <  2:3754426/91‑1 (MQ=255)
                                   ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                         >  2:3568814/1‑91 (MQ=255)
                                     aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                                       <  1:1074913/91‑1 (MQ=255)
                                      aaaaaTCCTTATCGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTa                                                      >  1:2991474/1‑91 (MQ=255)
                                       aaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTaa                                                     <  1:1005179/91‑1 (MQ=255)
                                        aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                                    >  2:416274/1‑91 (MQ=255)
                                        aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                                    >  2:358713/1‑91 (MQ=255)
                                        aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                                    >  2:3276113/1‑91 (MQ=255)
                                        aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                                    >  2:160838/1‑91 (MQ=255)
                                         aaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGa                                                   <  1:2846743/91‑1 (MQ=255)
                                           tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                                 >  2:3175005/1‑91 (MQ=255)
                                           tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                                 <  2:1456635/91‑1 (MQ=255)
                                             cTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTa                                               <  1:797298/91‑1 (MQ=255)
                                             cTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTa                                               <  1:2212933/91‑1 (MQ=255)
                                               tATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATg                                             >  1:1778701/1‑91 (MQ=255)
                                               tATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATg                                             <  1:3705208/91‑1 (MQ=255)
                                                aTTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGa                                            <  1:614087/91‑1 (MQ=255)
                                                aTTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGa                                            >  2:183249/1‑91 (MQ=255)
                                                 ttGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGAc                                           >  1:2440577/1‑91 (MQ=255)
                                                  tGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACg                                          <  2:766525/91‑1 (MQ=255)
                                                  tGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACg                                          <  1:310872/91‑1 (MQ=255)
                                                   gCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGt                                         <  1:1751392/91‑1 (MQ=255)
                                                     aTAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGttt                                       <  1:991031/91‑1 (MQ=255)
                                                       aaaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtgg                                     >  1:2266479/1‑91 (MQ=255)
                                                        aaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggt                                    <  1:2580410/91‑1 (MQ=255)
                                                         aTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggtg                                   >  1:323615/1‑91 (MQ=255)
                                                           gTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGTCTTATGACAATtggtgga                                 <  2:212091/91‑2 (MQ=255)
                                                            tACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCa                                >  1:1686593/1‑91 (MQ=255)
                                                              cATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCaaa                              >  2:2108680/1‑91 (MQ=255)
                                                                 aGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAg                           <  1:1487137/91‑1 (MQ=255)
                                                                   tgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGt                         <  2:3457490/91‑1 (MQ=255)
                                                                      agtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCaa                      >  1:1871401/1‑91 (MQ=255)
                                                                      agtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCaa                      <  2:3088060/91‑1 (MQ=255)
                                                                        taCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGa                    <  2:1548970/91‑1 (MQ=255)
                                                                         aCTATTCAAGACGTAATTATGACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGaa                   <  2:635175/91‑1 (MQ=255)
                                                                         aCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGaa                   >  1:1330787/1‑91 (MQ=255)
                                                                          cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt                  >  2:3131456/1‑91 (MQ=255)
                                                                          cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt                  >  1:2867653/1‑91 (MQ=255)
                                                                          cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt                  <  2:3791157/91‑1 (MQ=255)
                                                                            aTTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCa                <  1:1177998/91‑1 (MQ=255)
                                                                            aTTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCa                <  1:3833802/91‑1 (MQ=255)
                                                                              tCAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga              >  1:943868/1‑91 (MQ=38)
                                                                              tCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga              <  2:1651952/91‑1 (MQ=255)
                                                                                aaatcGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatg            >  1:2872510/5‑91 (MQ=39)
                                                                                aaatcGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatg            >  2:422727/5‑91 (MQ=39)
                                                                                 aaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatga           >  2:2888754/1‑91 (MQ=255)
                                                                                   aCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatg         <  1:1782932/91‑1 (MQ=255)
                                                                                   aCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatg         <  1:2675800/91‑1 (MQ=255)
                                                                                    cGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatga        >  1:2112970/1‑91 (MQ=255)
                                                                                    cGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatga        >  2:692128/1‑91 (MQ=255)
                                                                                     gTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatgat       >  2:146614/1‑91 (MQ=255)
                                                                                      tAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGata      >  1:2090989/1‑91 (MQ=255)
                                                                                       aaTTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat     <  1:2359473/91‑1 (MQ=255)
                                                                                       aaTTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat     >  2:3117885/1‑91 (MQ=255)
                                                                                       aaTTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat     <  1:3337993/91‑1 (MQ=255)
                                                                                       aaTGATTACTATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat     <  1:3488554/91‑1 (MQ=255)
                                                                                         ttattaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATgg   >  2:121775/1‑91 (MQ=255)
                                                                                         ttattaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATgg   >  2:612290/1‑91 (MQ=255)
                                                                                         ttattaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATgg   >  2:1139801/1‑91 (MQ=255)
                                                                                          tattaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGc  >  2:36637/1‑91 (MQ=255)
                                                                                          |                                                                                          
GACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGC  >  USA300TCH1516_ALE/79529‑79709

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: