Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,829,661 T→A 5.2% T487S (ACA→TCA)  ezrA ← Septation ring formation regulator EzrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,829,6610TA 5.2% 941.9 / 27.7 273T487S (ACA→TCA) ezrASeptation ring formation regulator EzrA
Reads supporting (aligned to +/- strand):  ref base T (123/135);  new base A (7/7);  total (130/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 279 total aligned reads displayed.

TATAGTCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCACGAAC  >  USA300TCH1516_ALE/1829571‑1829751
                                                                                          |                                                                                          
tataGTCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGt                                                                                            >  2:2830303/1‑91 (MQ=255)
tataGTCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGt                                                                                            >  1:554921/1‑91 (MQ=255)
 ataGTCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTa                                                                                           >  1:318877/1‑91 (MQ=255)
  taGTCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTAt                                                                                          >  1:1703187/1‑91 (MQ=255)
  taGTCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTAt                                                                                          <  1:3818391/91‑1 (MQ=255)
  taGTCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTAt                                                                                          >  1:1443721/1‑91 (MQ=255)
   aGTCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTAtt                                                                                         <  2:1145357/91‑1 (MQ=255)
    gTCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTc                                                                                        <  2:2858805/91‑1 (MQ=255)
     tCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCa                                                                                       <  2:2423946/91‑1 (MQ=255)
        ttACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGAATACAttt                                                                                    <  1:822990/91‑1 (MQ=255)
         tACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGAATACATTTg                                                                                   <  1:2337598/91‑1 (MQ=255)
         tACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATATTCAAATGAATACATTTg                                                                                   <  1:94729/91‑1 (MQ=255)
          aCGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGa                                                                                  >  2:3412902/1‑91 (MQ=255)
          aCGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGa                                                                                  >  1:494449/1‑91 (MQ=255)
          aCGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGa                                                                                  >  2:345352/1‑91 (MQ=255)
          aCGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGa                                                                                  >  1:2622966/1‑91 (MQ=255)
          aCGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGa                                                                                  >  2:2487145/1‑91 (MQ=255)
          aCGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGa                                                                                  >  2:1040648/1‑91 (MQ=255)
           cGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGaa                                                                                 <  1:869087/91‑1 (MQ=255)
           cGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGaa                                                                                 <  2:1243600/91‑1 (MQ=255)
            gATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAAt                                                                                <  2:585181/91‑1 (MQ=255)
            gATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAAt                                                                                <  1:264788/91‑1 (MQ=255)
            gATATCTATCTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAAt                                                                                <  1:3711361/91‑1 (MQ=255)
             atatCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATc                                                                               <  1:69901/91‑1 (MQ=255)
              tatCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATca                                                                              <  1:3435039/91‑1 (MQ=255)
                tCTATTTCCATATTGAATTAATTTCTCTGCATACACAGCATTAACAAGAACATCATTTGCTTCATCTCAAAATGTATTCATTTGAATcaca                                                                            >  2:3133366/1‑91 (MQ=255)
                tCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATcaca                                                                            >  2:1503392/1‑91 (MQ=255)
                tCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATcaca                                                                            >  1:112002/1‑91 (MQ=255)
                  tATTTCGATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAAt                                                                          <  2:552997/91‑1 (MQ=255)
                  tATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAAt                                                                          >  1:1343325/1‑91 (MQ=255)
                    tttCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAAttt                                                                        <  1:3761558/91‑1 (MQ=255)
                       ccATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAg                                                                     >  1:2604601/1‑91 (MQ=255)
                       ccATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAg                                                                     <  2:662226/91‑1 (MQ=255)
                          tatTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGaca                                                                  >  1:1096418/1‑91 (MQ=255)
                          tatTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGaca                                                                  >  1:1993379/1‑91 (MQ=255)
                          tatTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGaca                                                                  <  1:279440/91‑1 (MQ=255)
                           atTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGacac                                                                 >  1:2020097/1‑91 (MQ=255)
                                attaattTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTAt                                                            >  2:59780/1‑91 (MQ=255)
                                   aattTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTACTCACCTGAATCACAATTTTAGACACTCTATCtt                                                         >  2:136122/1‑91 (MQ=255)
                                     tttCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTCATCtttt                                                       >  1:1946791/1‑91 (MQ=255)
                                      ttCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTa                                                      >  2:1741253/1‑91 (MQ=255)
                                       tctcTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTaa                                                     <  2:1187096/91‑1 (MQ=255)
                                        ctctGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAAc                                                    <  2:1843800/91‑1 (MQ=255)
                                        ctctGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAAc                                                    >  1:367783/1‑91 (MQ=255)
                                        ctctGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAAc                                                    >  1:2452681/1‑91 (MQ=255)
                                          ctGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTg                                                  >  2:1124793/1‑91 (MQ=255)
                                            gCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGtt                                                >  2:3538079/1‑91 (MQ=255)
                                             cATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGttt                                               >  1:3146004/1‑91 (MQ=255)
                                               tAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTaa                                             <  2:1765840/91‑1 (MQ=255)
                                               tAAACAACATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTaa                                             >  1:3020504/1‑91 (MQ=255)
                                                  aCAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACgt                                          >  2:2176348/1‑91 (MQ=255)
                                                     gCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTa                                       <  1:2487145/91‑1 (MQ=255)
                                                     gCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTa                                       <  2:3034811/91‑1 (MQ=255)
                                                      cATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTAt                                      <  1:613835/91‑1 (MQ=255)
                                                      cATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTAt                                      >  2:2008300/1‑91 (MQ=255)
                                                      cATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTAt                                      >  2:197902/1‑91 (MQ=255)
                                                      cATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTAt                                      >  2:2657981/1‑91 (MQ=255)
                                                      cATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTAt                                      <  2:2936815/91‑1 (MQ=255)
                                                      cATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTAt                                      <  1:1069545/91‑1 (MQ=255)
                                                       aTTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTAtt                                     <  2:3621042/91‑1 (MQ=255)
                                                       aTTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTAtt                                     <  2:494449/91‑1 (MQ=255)
                                                         tAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTgg                                   <  1:3548648/91‑1 (MQ=255)
                                                           aCAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGAc                                 >  2:3013284/1‑91 (MQ=255)
                                                           aCAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGAc                                 >  2:3721045/1‑91 (MQ=255)
                                                            cAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACg                                <  1:3339710/91‑1 (MQ=255)
                                                             aaGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGGCACTTTATCTTTTAACTGTTTAACGTGTATTGGACGt                               <  1:3289324/91‑1 (MQ=255)
                                                             aaGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGt                               <  1:1582712/91‑1 (MQ=255)
                                                               gAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTc                             <  1:1265713/91‑1 (MQ=255)
                                                                aaCATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCa                            <  2:21880/91‑1 (MQ=255)
                                                                  catcatTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACt                          <  2:1343325/91‑1 (MQ=255)
                                                                    tcatTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTaa                        <  1:2374575/91‑1 (MQ=255)
                                                                     catTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTaaa                       >  2:2396754/1‑91 (MQ=255)
                                                                     catTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTaaa                       >  1:576538/1‑91 (MQ=255)
                                                                       tttGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAAtt                     <  1:1384875/91‑1 (MQ=255)
                                                                        ttGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTg                    >  1:761113/1‑91 (MQ=255)
                                                                        ttGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTg                    >  2:2541484/1‑91 (MQ=255)
                                                                        ttGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGATCACTAAATTg                    >  2:3337760/1‑91 (MQ=255)
                                                                         tGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGt                   >  2:1713514/1‑91 (MQ=255)
                                                                          gCTTCATCTTCGAATGTATTCAGCTGAATCACAATTTTAGACACTTTATGTTTTAACTGGTTAACGTGTATTGGACGTTCACTAAATTGtt                  <  2:2955948/91‑1 (MQ=255)
                                                                            ttCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCg                <  2:1935992/91‑1 (MQ=255)
                                                                             tCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGt               <  1:1953389/91‑1 (MQ=255)
                                                                             tCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGt               <  2:2548021/91‑1 (MQ=255)
                                                                                tCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTaa            >  2:3080864/1‑91 (MQ=255)
                                                                                tCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTaa            <  2:797436/91‑1 (MQ=255)
                                                                                gattCAAATGAATACATTTGAATCACACTTTTAGACACCTTATCTTTTAACTGTTTAACGTGTATTGGACGCTCACTACATTGTTCGTTaa            >  2:3265019/3‑91 (MQ=255)
                                                                                 cTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAAc           <  2:387497/91‑1 (MQ=255)
                                                                                   tCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACAt         <  1:1564995/91‑1 (MQ=255)
                                                                                   tCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACAt         <  2:1542879/91‑1 (MQ=255)
                                                                                   tCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACAt         >  1:1128399/1‑91 (MQ=255)
                                                                                    cAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATc        >  1:1240183/1‑91 (MQ=255)
                                                                                     aaaTGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCa       >  1:3647356/1‑91 (MQ=255)
                                                                                      aaTGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCAc      <  1:2098019/91‑1 (MQ=255)
                                                                                      aaTGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCAc      <  2:2609024/91‑1 (MQ=255)
                                                                                       aTGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCACg     <  2:983813/91‑1 (MQ=255)
                                                                                       aTGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCACg     <  1:2947659/91‑1 (MQ=255)
                                                                                       aTGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCACg     >  1:2834096/1‑91 (MQ=255)
                                                                                        tGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCACGa    <  1:2745868/91‑1 (MQ=255)
                                                                                        tGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCACGa    >  2:3389118/1‑91 (MQ=255)
                                                                                         gTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCACGaa   <  2:2331415/91‑1 (MQ=255)
                                                                                         gTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACCTTCACTAAATTGTTCGTTAACATCACGaa   >  1:3115465/1‑91 (MQ=255)
                                                                                          |                                                                                          
TATAGTCCTTACGATATCTATTTCCATATTGAATTAATTTCTCTGCATAAACAGCATTAACAAGAACATCATTTGCTTCATCTTCAAATGTATTCATTTGAATCACAATTTTAGACACTTTATCTTTTAACTGTTTAACGTGTATTGGACGTTCACTAAATTGTTCGTTAACATCACGAAC  >  USA300TCH1516_ALE/1829571‑1829751

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: ‑