Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,702,674 G→C 6.1% F116L (TTC→TTG tylM1 ← dTDP‑3‑amino‑3,6‑dideoxy‑alpha‑D‑glucopyranose N,N‑dimethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,702,6740GC 6.1% 691.8 / 28.3 197F116L (TTC→TTGtylM1dTDP‑3‑amino‑3,6‑dideoxy‑alpha‑D‑glucopyranose N,N‑dimethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (107/78);  new base C (2/10);  total (109/88)
Fisher's exact test for biased strand distribution p-value = 6.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 216 total aligned reads displayed.

TTGATTATTAAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCAAA  >  USA300TCH1516_ALE/1702584‑1702763
                                                                                          |                                                                                         
ttGATTATTAAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATtgatg                                                                                           >  1:387689/1‑91 (MQ=255)
 tGATTATTAAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATtgatga                                                                                          <  2:2704767/91‑1 (MQ=255)
  gATTATTAAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGGCAGATGATGATAAACATTGATGaa                                                                                         <  1:1558598/91‑1 (MQ=255)
      attaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGtt                                                                                     >  2:965633/1‑91 (MQ=255)
      attaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGtt                                                                                     >  2:3044336/1‑91 (MQ=255)
      attaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATGtt                                                                                     <  1:3129960/91‑1 (MQ=255)
      attaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATGtt                                                                                     <  1:2692078/91‑1 (MQ=255)
       ttaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGttt                                                                                    <  1:1399832/91‑1 (MQ=255)
       ttaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATGttt                                                                                    <  1:1237950/91‑1 (MQ=255)
       ttaAACAAAGTCATCATTTTATAAACAGTATATACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGttt                                                                                    <  2:1884300/91‑1 (MQ=255)
        taAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTc                                                                                   >  2:1307876/1‑91 (MQ=255)
        taAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTc                                                                                   <  2:1561515/91‑1 (MQ=255)
        taAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTc                                                                                   <  2:3333287/91‑1 (MQ=255)
            cAAAGTCATCATTTTATAAGCAGCATGTACATCAAAAATAAATACACCTGAATCAGTCAGATCATGATAAAGATTGGTGAATGTTTCAATc                                                                               <  2:1159338/91‑1 (MQ=255)
             aaaGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATcg                                                                              >  2:3852490/1‑90 (MQ=255)
              aaGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCAc                                                                             <  1:1676335/91‑1 (MQ=255)
               aGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACt                                                                            >  1:3216711/1‑91 (MQ=255)
               aGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACt                                                                            >  2:2194509/1‑91 (MQ=255)
               aGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACt                                                                            >  2:3285397/1‑91 (MQ=255)
                 tcatcaTTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGc                                                                          >  2:2547605/1‑91 (MQ=255)
                  catcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCa                                                                         >  1:3542081/1‑91 (MQ=255)
                   atcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGtt                                                                                     >  2:2248613/1‑78 (MQ=255)
                   atcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAg                                                                        >  2:1953010/1‑91 (MQ=255)
                   atcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAg                                                                        <  2:425126/91‑1 (MQ=255)
                    tcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGt                                                                       >  2:2870919/1‑91 (MQ=255)
                    tcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGt                                                                       >  2:1097226/1‑91 (MQ=255)
                    tcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGt                                                                       <  1:1563578/91‑1 (MQ=255)
                     catTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATTATGATAAACATTGATGAATGTTTCAATCACTGCAGtt                                                                      <  1:3276691/91‑1 (MQ=255)
                      atTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGttt                                                                     >  2:3883204/1‑91 (MQ=255)
                       ttttATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTc                                                                    >  1:3107849/1‑91 (MQ=255)
                        tttATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCa                                                                   >  2:726743/1‑91 (MQ=255)
                            taAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCtt                                                               >  2:1799451/1‑91 (MQ=255)
                            taAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCtt                                                               >  1:1768629/1‑91 (MQ=255)
                            taAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCtt                                                               >  2:2647724/1‑91 (MQ=255)
                            taAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTATCATCtt                                                               >  2:3441893/1‑91 (MQ=255)
                               aCAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCa                                                            >  1:2531386/1‑91 (MQ=255)
                               aCAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCa                                                            >  1:2247023/1‑91 (MQ=255)
                                cAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCaa                                                           >  2:1801059/1‑91 (MQ=255)
                                cAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCaa                                                           >  2:923752/1‑91 (MQ=255)
                                 aGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCaaa                                                          <  1:898114/91‑1 (MQ=255)
                                  gTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAaat                                                         >  1:2961390/1‑91 (MQ=255)
                                   tATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAaata                                                        >  1:1444227/1‑91 (MQ=255)
                                     tGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAaataat                                                      <  1:1125359/91‑1 (MQ=255)
                                       tACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAAttt                                                    <  1:3794118/91‑1 (MQ=255)
                                        aCATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTa                                                   >  2:2321519/1‑91 (MQ=255)
                                        aCATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTa                                                   <  2:580714/91‑1 (MQ=255)
                                        aCATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTa                                                   >  2:3019770/1‑91 (MQ=255)
                                         cATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTag                                                  <  2:947409/91‑1 (MQ=255)
                                         cATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTag                                                  >  1:3078227/1‑91 (MQ=255)
                                          aTCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTaga                                                 <  1:1506382/91‑1 (MQ=255)
                                               aaataaatACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTc                                            >  2:2130178/1‑91 (MQ=255)
                                               aaataaatACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTc                                            >  1:1825538/1‑91 (MQ=255)
                                                aataaataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTca                                           >  2:1820180/1‑91 (MQ=255)
                                                 ataaataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCa                                                            <  1:300320/73‑1 (MQ=255)
                                                 ataaataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTcac                                          <  1:196575/91‑1 (MQ=255)
                                                 ataaataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTcac                                          >  2:1193267/1‑91 (MQ=255)
                                                  taaataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTcaca                                         >  2:3091173/1‑91 (MQ=255)
                                                  taaataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTcaca                                         <  1:1799451/91‑1 (MQ=255)
                                                   aaataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaa                                        <  2:1782236/91‑1 (MQ=255)
                                                    aataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaa                                       >  1:128184/1‑91 (MQ=255)
                                                    aataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaa                                       >  2:14396/1‑91 (MQ=255)
                                                    aataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaa                                       >  1:455228/1‑91 (MQ=255)
                                                    aataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaa                                       >  1:2772474/1‑91 (MQ=255)
                                                     ataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaaa                                      >  2:2813832/1‑91 (MQ=255)
                                                     ataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaaa                                      >  1:87393/1‑91 (MQ=255)
                                                     ataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaaa                                      <  2:3322128/91‑1 (MQ=255)
                                                     ataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaaa                                      <  2:1396912/91‑1 (MQ=255)
                                                       acacCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaaaaa                                    >  2:1020158/1‑91 (MQ=255)
                                                       acacCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaaa                                      >  1:2103223/1‑89 (MQ=255)
                                                        cacCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAAt                                   >  1:561484/1‑91 (MQ=255)
                                                          ccTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAg                                 >  1:2902939/1‑91 (MQ=255)
                                                          ccTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAg                                 <  1:1307876/91‑1 (MQ=255)
                                                            tGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGtg                               <  2:1638252/91‑1 (MQ=255)
                                                            tGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAATCAAAAAAAATAGtg                               <  2:3590830/91‑1 (MQ=255)
                                                               aTCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGtgatg                            <  1:1412768/91‑1 (MQ=255)
                                                                 cagtcagATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAAACACAAAAAATAGTgatgat                          <  2:2034999/91‑1 (MQ=255)
                                                                  agtcagATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGata                         >  1:712265/1‑91 (MQ=255)
                                                                  agtcagATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGata                         >  1:3066721/1‑91 (MQ=255)
                                                                  agtcagATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCACCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGata                         >  2:3553651/1‑91 (MQ=255)
                                                                  agtcagATGATGATAAACATTGATGAATGTTCCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGata                         >  2:2953745/1‑91 (MQ=255)
                                                                     cagATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCa                      <  2:792927/91‑1 (MQ=255)
                                                                       gatgatgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCaaa                    <  1:2592186/91‑1 (MQ=255)
                                                                           atgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgt                <  2:1662255/91‑1 (MQ=255)
                                                                           atgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgt                <  1:219780/91‑1 (MQ=255)
                                                                            tgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgtt               >  2:3275285/1‑91 (MQ=255)
                                                                            tgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgtt               <  2:3726369/91‑1 (MQ=255)
                                                                             gatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgttg              >  1:281528/1‑91 (MQ=255)
                                                                             gatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgttg              >  1:239276/1‑91 (MQ=255)
                                                                             gatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgttg              >  2:3114349/1‑91 (MQ=255)
                                                                              atAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgttgt             >  1:2291985/1‑91 (MQ=255)
                                                                                aaaCATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAAGAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgttgttg           <  2:3485755/91‑1 (MQ=255)
                                                                                 aaCATTGATGAATGTTTCAATCACTGCAGTTTCCTCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGc          >  1:263416/1‑91 (MQ=255)
                                                                                 aaCATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGc          >  1:24658/1‑91 (MQ=255)
                                                                                 aaCATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGc          >  2:1498146/1‑91 (MQ=255)
                                                                                  aCATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCa         >  1:3376348/1‑91 (MQ=255)
                                                                                  aCATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAATAGTGATGATATCAAATTGTTGTTGCaa        <  2:1252363/91‑1 (MQ=255)
                                                                                      tgatgaCTGTTTCAATCACTGCAGTTTCATCTTGCAAATACTTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAACTc     >  2:489638/1‑91 (MQ=255)
                                                                                       gatgaATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCa    <  2:3376348/91‑1 (MQ=255)
                                                                                       gatgaATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCa    <  1:2130191/91‑1 (MQ=255)
                                                                                       gatgaATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCa    <  2:3138013/91‑1 (MQ=255)
                                                                                          |                                                                                         
TTGATTATTAAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCAAA  >  USA300TCH1516_ALE/1702584‑1702763

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