Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 590,688 G→A 5.4% G302S (GGT→AGT)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE590,6880GA 5.4% 196.8 / 3.0 75G302S (GGT→AGT) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (36/34);  new base A (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-01

TTATACCGTCGTGTAATTAATCGAAATAATCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAGTTATTGCAGTAGGT  >  USA300TCH1516_ALE/590553‑590825
                                                                                                                                       |                                                                                                                                         
ttATAACGTCGTGTAAGTTATGGAAATAATCGATTGAATCGTTTATTAGATTTAAGTGCTCCTGGTATCATCGTTCAAAACGAAATTCGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGgtcg                                                                                                                                       <  2:548791/140‑1 (MQ=255)
      cgtcgtGTAATTAATCGAAATAATCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAATCGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGt                                                                                                                                 <  2:299663/140‑1 (MQ=255)
      cgtcgtGTAATTAATCGAAATAATCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGt                                                                                                                                 <  1:500594/140‑1 (MQ=255)
      cgtcgtGTAATTAATCGAAATAATCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGt                                                                                                                                 <  1:219313/140‑1 (MQ=255)
        tcttgtAATTAATTTATATTACCGTTTGAATCGTTTATTAGATTTAGGTGCAGCTGGTATCATCGTTCACATCGAAAAACGTATGTTACAAGCAGCCGTTGCCGCCTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTa                                                                                                                               <  2:393002/137‑1 (MQ=255)
           tgtAATTAATCGAAATAATCGTTTGAAACGTTCATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTAGTACAGTTACTg                                                                                                                            >  2:174432/1‑140 (MQ=255)
            gtAATTAATCGAAATAATCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTgg                                                                                                                           >  1:625281/1‑140 (MQ=255)
                 taatCGAAATAATCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTAATAATGGTCGTCGTGGTCGTCCAGTTACTGGCACAg                                                                                                                      >  1:14612/1‑140 (MQ=255)
                    tctaaaTAATCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTa                                                                                                                   <  1:39074/137‑1 (MQ=255)
                     cGAAATAGTCGATCGAGACGTTTATTAGATTTAGGTGCACTTGGTATAATCGTTCATAACGAAAAACGTATGTTACAAGAAGCCGTGGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTaa                                                                                                                  <  2:299713/140‑1 (MQ=255)
                      gcaatattCTTTTGAAGCGTTTATTAGATTTAGGTGCACCGGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGTCTCATGTAAc                                                                                                                 <  2:87751/138‑1 (MQ=255)
                      gagagaatCGTTAGGAGGGTTTATGAGAGTTTGGTGCGCGTGGTAACATCGTGAAAAACGAAATAGGTATGTTACAAGAAGCCGGTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAAc                                                                                                                 <  2:665295/135‑1 (MQ=255)
                        aataatCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCg                                                                                                               >  1:71540/1‑140 (MQ=255)
                            atCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCa                                                                                                           >  1:193000/1‑140 (MQ=255)
                             tCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAGCGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTAGTCGTCCAGTTACTGGCCCAGGTAACCGTCCAt                                                                                                          >  2:346395/1‑140 (MQ=255)
                             tCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTAGTCGTCCAGTTACTGGCCCAGGTAACCGTCCAt                                                                                                          >  2:346402/1‑140 (MQ=255)
                              cttttAAATCGTTTATTAGATTTAGGTGCATCTGGTATCATCGTTAAACGCTAAATTCGTTTTTTGCTGGGTGCCGTGGACGCTTTAATTGAGATTGGGCGTCTTTGTCGTCCAGTTACTGGCCCAGGTATTCGTctagt                                                                                                         <  2:114917/138‑5 (MQ=255)
                              cGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCAtt                                                                                                         <  2:554715/140‑1 (MQ=255)
                                 ttGAAACGTTTACTAGATTTATGTTCCTCTGGTATCAGCGTTGAAAACGAAAAACGTATTTTACTAGAAGCCTTTGAGGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTaaa                                                                                                      <  2:366169/140‑1 (MQ=255)
                                     aaagTTTATTAAATTTAGGTGTAACTGGTATCATCGTTCAGAATGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCt                                                                                                  <  2:68568/137‑1 (MQ=255)
                                     aaCGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCt                                                                                                  >  2:511135/1‑140 (MQ=255)
                                        gTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTa                                                                                               >  1:400741/1‑140 (MQ=255)
                                            attaGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTAtctc                                                                                           >  1:439011/1‑140 (MQ=255)
                                               aggtttAGGTGCACCTGTGATCATCCGGCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCata                                                                                        <  2:134807/137‑1 (MQ=255)
                                                  tttAGGTGCACCTGGTATCAACGATCAAAAAAGATAGCGGATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCATCGTGGCAGTCCAGTTACTGGCCCAGTGAACCGTCAATTAAACTCCTTAGCTCATATGt                                                                                     >  2:318597/1‑140 (MQ=255)
                                                    tAGGTGCACCTGGTATCATCGTTCAAAACTAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTa                                                                                   <  2:28140/140‑1 (MQ=255)
                                                    tAGGTGCACCTGGGATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTa                                                                                   <  1:254686/140‑1 (MQ=255)
                                                     aGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTaa                                                                                  >  1:533912/1‑140 (MQ=255)
                                                     aGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTACAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGATaa                                                                                  >  2:362220/1‑140 (MQ=255)
                                                         gCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAgg                                                                              <  1:603814/140‑1 (MQ=255)
                                                              tGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATGAAAATCTTTATCTCATAGTTTAAGAGGTAAAc                                                                         >  2:492843/1‑140 (MQ=255)
                                                               ggTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGCCGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTACCATCTCTATCTCATATCCTAAAAGGtcacca                                                                        >  2:425935/1‑135 (MQ=255)
                                                               ggTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAACTCTTTATCTCATATGTTAAACGGTAAACa                                                                        >  2:197359/1‑140 (MQ=255)
                                                               ggTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCCACTATGTTCAAAGGTAGACa                                                                        >  2:551994/1‑140 (MQ=255)
                                                               ggTATCATCGTTCAAAACCAAAAACGTATGTTACAAGAAGCCGTTGCCGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACa                                                                        <  1:347643/140‑1 (MQ=255)
                                                                 tATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAg                                                                      >  1:583413/1‑140 (MQ=255)
                                                                  atcatcGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAgg                                                                     <  2:71540/140‑1 (MQ=255)
                                                                   tcatcGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGt                                                                    <  1:615668/140‑1 (MQ=255)
                                                                   tcatcGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCATAGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTGATCTCATATGTTAAAACGTAACCAAGGt                                                                    >  2:471143/1‑140 (MQ=255)
                                                                      tcGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATATCATATGTTAAAAGGAAAACAAGGACGt                                                                 >  2:424788/1‑140 (MQ=255)
                                                                              aaCGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTCTCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCa                                                         >  1:175600/1‑140 (MQ=255)
                                                                                 gAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCaaac                                                      >  2:107768/1‑139 (MQ=255)
                                                                                 gAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTAGTACAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATTTCATATATTAAAAGGAAAACAAGATAGTGTCCGGCaaaa                                                      >  2:97676/1‑140 (MQ=255)
                                                                                 gAAAAACGTATGATAAAAGAAGCCGTTGACGCTTTAATTGATAGTGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCaaaa                                                      >  2:319069/1‑140 (MQ=255)
                                                                                  aaaacCGTATGTTACAAGAAGCCGTTGACGCTTTAATTCATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCCTATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAac                                                     <  1:53914/140‑1 (MQ=255)
                                                                                      aCGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTAGTTTCCGTCAAAacttac                                                 >  2:83160/1‑140 (MQ=255)
                                                                                        gTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATCTGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAacttactt                                               >  1:417966/1‑140 (MQ=255)
                                                                                         tATGTTACAAGAAGCCGTTAGCGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAAGCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTg                                              <  2:57825/140‑1 (MQ=255)
                                                                                          aTGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTgg                                             >  2:654434/1‑140 (MQ=255)
                                                                                           tGTTACAAGAAGCCGTTGACGCTTTAATTGATACTGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGt                                            <  1:367803/140‑1 (MQ=255)
                                                                                             ttACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTaa                                          >  1:276467/1‑140 (MQ=255)
                                                                                              tACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTaaa                                         <  1:196251/140‑1 (MQ=255)
                                                                                                   gaagCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGATCGTTTCCGTCAAAACTTACTTGGTAAACgtgt                                    >  2:322971/1‑140 (MQ=255)
                                                                                                     agCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTg                                  <  1:502761/140‑1 (MQ=255)
                                                                                                            gACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCGAAGTAACCGTTCAGTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTc                           >  1:267995/1‑140 (MQ=255)
                                                                                                               gCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAgg                        >  1:591445/1‑140 (MQ=255)
                                                                                                                    aaTTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAAGCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGtt                   <  1:54746/140‑1 (MQ=255)
                                                                                                                     aTTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTc                  <  2:484573/140‑1 (MQ=255)
                                                                                                                      ttGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCa                 <  2:193000/140‑1 (MQ=255)
                                                                                                                      ttGATAATGGTCGTCGTAGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCa                 <  1:346402/140‑1 (MQ=255)
                                                                                                                      ttGATAATGGTCGTCGTAGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCa                 <  1:346395/140‑1 (MQ=255)
                                                                                                                       tGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAACATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAg                <  1:229105/140‑1 (MQ=255)
                                                                                                                       tGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAg                <  1:378106/140‑1 (MQ=255)
                                                                                                                          tAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAACTCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGACAAAACTACCTTGGTAAACGTGTTGACTATTCCGGCCGtcacgtac             >  2:73575/1‑133 (MQ=255)
                                                                                                                            aTGGTTGTCGTGATCGTACAGTTAGTGGCCCAGGTAACCGACCAGTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAGTTAtt           <  1:33939/140‑1 (MQ=255)
                                                                                                                              ggTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAGCTAtttc         >  2:331864/1‑138 (MQ=255)
                                                                                                                                tcgtcgTGGTTGTCTAGTTATTGGCCCAGGTAACAGTCCATTAAAATCTTTATATCATATGGTAAGAGGTAAACAAGGTCGTCTCCGTCAAAACTTACTTGGTAAACGTGTTGACCATTCAGGACGTTCAGTTATTGCAg       <  1:318597/140‑1 (MQ=255)
                                                                                                                                 cgtcgtGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAGTTATTGCAGt      >  1:111505/1‑140 (MQ=255)
                                                                                                                                 cgtcgtGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAGTTATTGCAGt      >  1:592745/1‑140 (MQ=255)
                                                                                                                                 cgtcgtGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAGTTATTCCAGt      >  1:494135/1‑140 (MQ=255)
                                                                                                                                 cgtcgtGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTCGAGTTATAGCAgc      >  1:282885/1‑139 (MQ=255)
                                                                                                                                 cgtcgtGGTCGTCCAGTTACGGGCCGAGGTAACCGTCGATGAAAATCTTGATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGGTCAGTTATTGCatt      >  1:167920/1‑138 (MQ=255)
                                                                                                                                    cgtGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAGTTATTGCAGTAgg   <  1:6886/140‑1 (MQ=255)
                                                                                                                                     ttGGTCGGCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAGTTATTGCAGTAGGt  <  1:55783/139‑1 (MQ=255)
                                                                                                                                     gtGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTAGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAGTTATTGCAGTAGGt  <  1:57039/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
TTATACCGTCGTGTAATTAATCGAAATAATCGTTTGAAACGTTTATTAGATTTAGGTGCACCTGGTATCATCGTTCAAAACGAAAAACGTATGTTACAAGAAGCCGTTGACGCTTTAATTGATAATGGTCGTCGTGGTCGTCCAGTTACTGGCCCAGGTAACCGTCCATTAAAATCTTTATCTCATATGTTAAAAGGTAAACAAGGTCGTTTCCGTCAAAACTTACTTGGTAAACGTGTTGACTATTCAGGACGTTCAGTTATTGCAGTAGGT  >  USA300TCH1516_ALE/590553‑590825

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