Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,672,412 G→A 66.0% A194V (GCA→GTA)  sigA ← RNA polymerase sigma factor SigA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,672,4120GA66.0% 62.1 / 60.1 68A194V (GCA→GTA) sigARNA polymerase sigma factor SigA
Reads supporting (aligned to +/- strand):  ref base G (8/15);  new base A (14/30);  total (23/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

ATTTCTTCTGGTGCTGGATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTACGTATCTTTT  >  USA300TCH1516_ALE/1672277‑1672551
                                                                                                                                       |                                                                                                                                           
aTTTCTTCTGGTGCTGGATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGAGTACTTGt                                                                                                                                        <  1:437006/141‑1 (MQ=255)
   tcttctGGTGCTGGATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTa                                                                                                                                     >  1:674957/1‑141 (MQ=255)
     ttctGGTGCTGGATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAAt                                                                                                                                   >  2:631042/1‑141 (MQ=255)
       ctGGTGCTGGATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATcc                                                                                                                                 <  1:20289/141‑1 (MQ=255)
        tGGTGCTGGATCTCGACCTAAGTCCTGTCATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATCTGCACCGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATcca                                                                                                                                <  1:482015/141‑1 (MQ=255)
        tGGTGCTGGATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATcca                                                                                                                                <  1:116665/141‑1 (MQ=255)
         ggTGCTGGATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATccac                                                                                                                               <  1:89401/141‑1 (MQ=255)
          gTGCTGGATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTCTTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATccacc                                                                                                                              <  1:572861/141‑1 (MQ=255)
             cTGGATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATg                                                                                                                           >  1:185310/1‑141 (MQ=255)
               ggATCTCGTCCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGtt                                                                                                                         <  1:643080/141‑1 (MQ=255)
               ggATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGtt                                                                                                                         >  1:8098/1‑141 (MQ=255)
               ggATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACCAACTAATTTCTAAATTCATTATAACATTTGCAACGGGGTCAGTAACATACGTGCTGGGGCAACAGGTGGACGAGTGAATGCGGGTCGAATCCACCATGtt                                                                                                                         >  2:576807/1‑141 (MQ=255)
                  tctcGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGGGCTTGGTAAGACATTGCACGAGCGATTAATTGTCTAGTCCCACATGTTGcg                                                                                                                      >  2:217585/1‑140 (MQ=255)
                      gACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATg                                                                                                                  >  1:415134/1‑141 (MQ=255)
                              tCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAAc                                                                                                          <  2:258687/141‑1 (MQ=255)
                                cTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACtt                                                                                                        <  2:302524/141‑1 (MQ=255)
                                   taataaTTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTaaa                                                                                                     <  1:475677/141‑1 (MQ=255)
                                       aatTGACGTTGAACGCGAAGTAATTGATTAATTGTTTCTACCATCTGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAACCTTAAATCCt                                                                                                 >  1:381908/1‑141 (MQ=255)
                                       aatTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCt                                                                                                 >  1:641566/1‑141 (MQ=255)
                                       aatTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCt                                                                                                 >  1:552620/1‑141 (MQ=255)
                                         ttGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCttt                                                                                               >  2:548705/1‑141 (MQ=255)
                                            aCGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGtt                                                                                            <  1:127450/141‑1 (MQ=255)
                                                tcaaCATGAGTTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAg                                                                                        <  2:172055/139‑1 (MQ=255)
                                                          ttaattTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCCCCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAAttttt                                                                              >  1:107512/1‑141 (MQ=255)
                                                          ttaattTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAAttttt                                                                              <  1:468310/141‑1 (MQ=255)
                                                           taattTGTTTATTGTTTCGACCATAAGTACAGGGATACTAATTGTACGTGCTTGGTTAGCGATTGCACGAGTGATTAATAGGCTTATACAGCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCATATTTTTc                                                                             <  2:329989/141‑1 (MQ=255)
                                                                    aaTTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCttt                                                                    <  2:347073/141‑1 (MQ=255)
                                                                    aaTTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCCTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCttt                                                                    <  2:92540/141‑1 (MQ=255)
                                                                     aTTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTta                                                                   <  2:499969/141‑1 (MQ=255)
                                                                     aTTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTta                                                                   <  2:687680/141‑1 (MQ=255)
                                                                     aTTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTta                                                                   <  2:481608/141‑1 (MQ=255)
                                                                     aTTGTTTCTACCATATGAACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTAATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTta                                                                   <  1:267863/141‑1 (MQ=255)
                                                                      ttGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTtaa                                                                  <  2:30132/141‑1 (MQ=255)
                                                                      ttGTTTCTAACATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGTAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCAGATGTTGAAAACTTAAATCCTTTCTTACAGCCAAACTTTTCAACAGCTTtaa                                                                  <  2:117157/141‑1 (MQ=255)
                                                                      ttGTTTCCACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTtaa                                                                  <  1:405114/141‑1 (MQ=255)
                                                                       tGTTTCTACTATAGGCACAGGTATACGAATCGTACGTGCTTGGGCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTtaat                                                                 <  2:293112/141‑1 (MQ=255)
                                                                       tGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTtaat                                                                 <  1:658028/141‑1 (MQ=255)
                                                                       tGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTtaat                                                                 <  2:534694/141‑1 (MQ=255)
                                                                                 atatGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCatt                                                       >  2:269213/1‑140 (MQ=255)
                                                                                  tatGCACAGGGATACGAATCGTACGTGCATGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCatat                                                      >  1:439974/1‑141 (MQ=255)
                                                                                      cacaGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTAcc                                                  <  1:262925/141‑1 (MQ=255)
                                                                                         aGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACcttc                                               <  1:478916/141‑1 (MQ=255)
                                                                                            gATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTg                                            <  1:39801/141‑1 (MQ=255)
                                                                                             aTACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAAAAAGACCCATATTACCTTCTTgg                                           <  2:493678/141‑1 (MQ=255)
                                                                                             aTACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTCCTTGTCTAATCCACCATGTTCCATCTGTTGAAAACTTAAATCCTTTCTTAAAGTCAAATTTTTCCACCGTTTTAATAACACCCATATTCCCTTCCTgt                                           >  2:433788/1‑140 (MQ=255)
                                                                                              tACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGa                                          <  2:188976/141‑1 (MQ=255)
                                                                                               aCGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGAt                                         <  2:185310/141‑1 (MQ=255)
                                                                                                  aaTCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTGGTGAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTaa                                      >  1:399989/1‑141 (MQ=255)
                                                                                                         cGTGCTTGGTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACGCATATTACCTTCTTGGATTAAATCAAgg                               >  2:630873/1‑141 (MQ=255)
                                                                                                         cGTGCTTGGTCAGCAATTGCACGAGTGATTACTGGTCTAGTCCACGATGTTGCATATCTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTGAACAGCTTTAATAAGATCCAGATTACCTTATTGAATTGAATCAAga                               >  2:326105/1‑140 (MQ=255)
                                                                                                         cGTGCTTGGTCAGCAATTGCACGAGTGATTACCTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAgg                               <  1:351286/141‑1 (MQ=255)
                                                                                                                 gTCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACAt                       <  2:415134/141‑1 (MQ=255)
                                                                                                                  tCAGCAATTGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATa                      >  1:620421/1‑141 (MQ=255)
                                                                                                                   cAGCAATTGCACGAGAGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTGCTTGGATTAAATCAAGGAATAAAATAc                     >  1:351662/1‑141 (MQ=255)
                                                                                                                        aTTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGa                <  1:242911/141‑1 (MQ=255)
                                                                                                                          tGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACAAcc              >  1:566611/1‑141 (MQ=255)
                                                                                                                          tGCACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTCCCTTCTTGGATTAAATCAAGGAATAACATACCACGAcc              >  1:16229/1‑141 (MQ=255)
                                                                                                                            cACGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTa            <  2:469034/141‑1 (MQ=255)
                                                                                                                             aCGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAATTTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTAc           <  2:121967/141‑1 (MQ=255)
                                                                                                                             aCGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTAc           <  2:464024/141‑1 (MQ=255)
                                                                                                                              cgggTGATTACTTGTCTAACCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCACATTTTTCAACAGTGGTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTACg          <  2:593282/138‑1 (MQ=255)
                                                                                                                              cGAGTGATTACTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTACg          <  2:302025/141‑1 (MQ=255)
                                                                                                                              cGAATGATTACTTGTCTAATCCACCTTGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTGTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTACg          <  2:669811/141‑1 (MQ=255)
                                                                                                                                 gTGATTGCTTGACTAATCCACCATGTTGCATATGTTGAAAACTTAAAGCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTACGTAt       <  2:514252/141‑1 (MQ=255)
                                                                                                                                  tGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTACGTATc      <  1:418657/141‑1 (MQ=255)
                                                                                                                                    aTTGCTTGTCTAATCCACCATGTTGCGTATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAACACCCAGATTACCTTCTTGGATTAAATCAAGGAATAACCTACCACGACCTACGTATCtt    <  1:423916/141‑1 (MQ=255)
                                                                                                                                    aTTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTACGTATCtt    <  2:612981/141‑1 (MQ=255)
                                                                                                                                      tGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCAAATAGGGGGCTCGGATTAAATCAAGGAATAACATACCACGACCTACGTATCtttt  >  1:677515/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                           
ATTTCTTCTGGTGCTGGATCTCGACCTAAGTCCTGTAATAATTGACGTTGAACACGAATTAATTTATTAATTGTTTCTACCATATGCACAGGGATACGAATCGTACGTGCTTGGTCAGCAATTGCACGAGTGATTGCTTGTCTAATCCACCATGTTGCATATGTTGAAAACTTAAATCCTTTGTTAAAGTCAAATTTTTCAACAGCTTTAATAAGACCCATATTACCTTCTTGGATTAAATCAAGGAATAACATACCACGACCTACGTATCTTTT  >  USA300TCH1516_ALE/1672277‑1672551

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