Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,702,674 G→C 10.2% 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,6740GC10.2% 387.0 / 27.2 117F116L (TTC→TTGtylM1dTDP‑3‑amino‑3,6‑dideoxy‑alpha‑D‑glucopyranose N,N‑dimethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (52/53);  new base C (3/9);  total (55/62)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

TTGATTATTAAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCAAA  >  USA300TCH1516_ALE/1702584‑1702763
                                                                                          |                                                                                         
ttGATTATTAAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATtgatg                                                                                           <  1:649709/91‑1 (MQ=255)
  gATTATTAAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGaa                                                                                         <  1:1312742/91‑1 (MQ=255)
  gATTATTAAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGaa                                                                                         <  1:1772995/91‑1 (MQ=255)
    ttattaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATg                                                                                       <  1:1038485/91‑1 (MQ=255)
    ttattaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATg                                                                                       >  1:1210257/1‑91 (MQ=255)
    ttattaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATg                                                                                       <  1:1111805/91‑1 (MQ=255)
       ttaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATGttt                                                                                    <  1:1121962/91‑1 (MQ=255)
       ttaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATGttt                                                                                    <  1:1040184/91‑1 (MQ=255)
       ttaAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATGttt                                                                                    <  1:1121953/91‑1 (MQ=255)
         aaaCAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATGtttat                                                                                  <  1:4103842/91‑3 (MQ=255)
         aaaCAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATGtttat                                                                                  <  1:4103862/91‑3 (MQ=255)
         aaaCAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATGtttat                                                                                  <  1:4103863/91‑3 (MQ=255)
           aCAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAAt                                                                                >  1:3477090/1‑91 (MQ=255)
           aCAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTCATCAATGtttatca                                                                                <  1:3939171/91‑5 (MQ=255)
            cAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATc                                                                               <  1:2790748/91‑1 (MQ=255)
             aaaGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCa                                                                              >  1:3660561/1‑91 (MQ=255)
             aaaGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCa                                                                              >  1:2179620/1‑91 (MQ=255)
             aaaGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCa                                                                              <  1:3469537/91‑1 (MQ=255)
              aaGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCAc                                                                             >  1:3425184/1‑91 (MQ=255)
              aaGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCAc                                                                             <  1:2660619/91‑1 (MQ=255)
              aaGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCAc                                                                             <  1:2233006/91‑1 (MQ=255)
              aaGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCAc                                                                             >  1:2114562/1‑91 (MQ=255)
              aaGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCAc                                                                             <  1:1630391/91‑1 (MQ=255)
               aGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACt                                                                            >  1:390491/1‑91 (MQ=255)
               aGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACt                                                                            >  1:392597/1‑91 (MQ=255)
               aGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACt                                                                            >  1:518172/1‑91 (MQ=255)
               aGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACt                                                                            >  1:516468/1‑91 (MQ=255)
                gTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTg                                                                           <  1:2507396/91‑1 (MQ=255)
                 tcatcaTTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGc                                                                          <  1:1611012/91‑1 (MQ=255)
                   atcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAg                                                                        >  1:2126020/1‑91 (MQ=255)
                   atcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAg                                                                        >  1:826154/1‑91 (MQ=255)
                   atcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAg                                                                        <  1:1163824/91‑1 (MQ=255)
                   atcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAg                                                                        >  1:1934269/1‑91 (MQ=255)
                   atcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAg                                                                        <  1:2916449/91‑1 (MQ=255)
                    tcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGt                                                                       <  1:2141457/91‑1 (MQ=255)
                    tcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGt                                                                       <  1:2548414/91‑1 (MQ=255)
                    tcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGt                                                                       >  1:3856983/1‑91 (MQ=255)
                    tcatTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACAGTGATGAATGTTTCAATCACTGCAGt                                                                       >  1:3855418/1‑91 (MQ=255)
                     catTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGtt                                                                      <  1:1193600/91‑1 (MQ=255)
                     catTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGtt                                                                      <  1:3829763/91‑1 (MQ=255)
                      atTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGttt                                                                     <  1:2810677/91‑1 (MQ=255)
                      atTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGttt                                                                     <  1:2810428/91‑1 (MQ=255)
                      atTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGttt                                                                     <  1:290375/91‑1 (MQ=255)
                              aaCAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGc                                                             <  1:814319/91‑1 (MQ=255)
                               aCAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCa                                                            >  1:882069/1‑91 (MQ=255)
                                 aGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCaaa                                                          <  1:3049135/91‑1 (MQ=255)
                                 aGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCaaa                                                          <  1:1874554/91‑1 (MQ=255)
                                  gTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAaat                                                         >  1:564798/1‑91 (MQ=255)
                                      gTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAAtt                                                     >  1:1125750/1‑91 (MQ=255)
                                       tACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAAttt                                                    >  1:2672732/1‑91 (MQ=255)
                                       tACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAAttt                                                    <  1:1637021/91‑1 (MQ=255)
                                        aCATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTa                                                   >  1:3242824/1‑91 (MQ=255)
                                           tCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTagag                                                <  1:3645335/91‑1 (MQ=255)
                                            cAAAAATAAATACACCTGAATCAGTCAGATGATGATACACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTagaga                                               >  1:3504488/1‑91 (MQ=255)
                                            cAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTagaga                                               <  1:1275598/91‑1 (MQ=255)
                                            cAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTagaga                                               >  1:3506943/1‑91 (MQ=255)
                                            cAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTagaga                                               >  1:3503945/1‑91 (MQ=255)
                                              aaaaTAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCCGTTTCATCTTGCAAATAATTCAGAGAGt                                             >  1:1411530/1‑91 (MQ=255)
                                              aaaaTAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGt                                             <  1:342983/91‑1 (MQ=255)
                                               aaataaatACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTc                                            >  1:156452/1‑91 (MQ=255)
                                               aaataaatACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTc                                            >  1:156571/1‑91 (MQ=255)
                                                 ataaataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTcac                                          <  1:2559935/91‑1 (MQ=255)
                                                    aataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaa                                       >  1:3931687/1‑91 (MQ=255)
                                                    aataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaa                                       >  1:2841053/1‑91 (MQ=255)
                                                    aataCACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaa                                       >  1:2841854/1‑91 (MQ=255)
                                                    aataCACCTGAATCAGTCAGATGATGATAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaaa                                      <  1:1881592/91‑1 (MQ=255)
                                                       acacCTGAATCAGTCAGATGATGATAAACATTGATGATTGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACaaaaaa                                    >  1:3497371/1‑91 (MQ=255)
                                                          ccTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAg                                 <  1:2230729/91‑1 (MQ=255)
                                                           cTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGt                                >  1:251368/1‑91 (MQ=255)
                                                           cTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGt                                <  1:3706505/91‑1 (MQ=255)
                                                                tcagtcagATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGtgatga                           <  1:2976996/91‑1 (MQ=255)
                                                                tcagtcagATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGtgatga                           <  1:3696262/91‑1 (MQ=255)
                                                                   gtcagATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGatat                        <  1:2424362/91‑1 (MQ=255)
                                                                   gtcagATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGatat                        <  1:2774919/91‑1 (MQ=255)
                                                                      agATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATATTTTAGAGAGTCACAAAAAATAGTGATGATATCaa                     >  1:361800/1‑91 (MQ=255)
                                                                      agATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCaa                     >  1:317067/1‑91 (MQ=255)
                                                                       gatgatgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCaaa                    >  1:2566320/1‑91 (MQ=255)
                                                                       gatgatgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCaaa                    >  1:2566436/1‑91 (MQ=255)
                                                                       gatgatgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCaaa                    >  1:2052405/1‑91 (MQ=255)
                                                                       gatgatgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCaaa                    <  1:165199/91‑1 (MQ=255)
                                                                         tgatgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAtt                  >  1:679092/1‑91 (MQ=255)
                                                                         tgatgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAtt                  <  1:954672/91‑1 (MQ=255)
                                                                            tgatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgtt               <  1:1569027/91‑1 (MQ=255)
                                                                            attgaaaCATTCATCAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgtt               >  1:3029304/5‑91 (MQ=255)
                                                                            attgaaaCATTCATCAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgtt               >  1:3024125/5‑91 (MQ=255)
                                                                             gatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgttg              >  1:3233763/1‑91 (MQ=255)
                                                                             gatAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgttg              >  1:505876/1‑91 (MQ=255)
                                                                               tAAACATTGATGAATGTTTCAATCACTGCTGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgttgtt            <  1:4111450/91‑1 (MQ=255)
                                                                               tAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgttgtt            >  1:2379390/1‑91 (MQ=255)
                                                                                aaaCATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAAttgttgttg           <  1:1248935/91‑1 (MQ=255)
                                                                                 aaCATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGc          <  1:3305401/91‑1 (MQ=255)
                                                                                  aCATTCATCAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCa         >  1:1872096/1‑91 (MQ=255)
                                                                                   cATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCaa        <  1:2842785/91‑1 (MQ=255)
                                                                                   cATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCaa        >  1:832991/1‑91 (MQ=255)
                                                                                    aTTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCaaa       <  1:3796161/91‑1 (MQ=255)
                                                                                     ttGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAAt      <  1:2118724/91‑1 (MQ=255)
                                                                                        atgaATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCaa   <  1:2662175/91‑1 (MQ=255)
                                                                                        atgaATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCaa   >  1:2467549/1‑91 (MQ=255)
                                                                                         tgaATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCaaa  <  1:82503/91‑1 (MQ=255)
                                                                                         tgaATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCaaa  >  1:2620607/1‑91 (MQ=255)
                                                                                          |                                                                                         
TTGATTATTAAACAAAGTCATCATTTTATAAACAGTATGTACATCAAAAATAAATACACCTGAATCAGTCAGATGATGATAAACATTGATGAATGTTTCAATCACTGCAGTTTCATCTTGCAAATAATTTAGAGAGTCACAAAAAATAGTGATGATATCAAATTGTTGTTGCAAATCAAA  >  USA300TCH1516_ALE/1702584‑1702763

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: