Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 973,213 C→A 72.2% E175* (GAG→TAG)  yitU_1 ← 5‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase YitU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE973,2130CA72.2% 99.7 / 54.6 79E175* (GAG→TAG) yitU_15‑amino‑6‑(5‑phospho‑D‑ribitylamino)uracil phosphatase YitU
Reads supporting (aligned to +/- strand):  ref base C (15/7);  new base A (30/26);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

CTTCATCACCGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATTACCCAT  >  USA300TCH1516_ALE/973078‑973351
                                                                                                                                       |                                                                                                                                          
tgatgtcaCCGAATGCAATAATATTATTTCGGTCAATATTTAAAGATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAACTTTTACAATTTCAATGACATGGAATGGTGCGCCCCAGCGTCGATGCTAAatat                                                                                                                                       >  2:656692/6‑141 (MQ=255)
tgatgtcaCCGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAatat                                                                                                                                       >  1:235211/6‑141 (MQ=255)
     tcaCCGAATGCAATAAGATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCTCCAGCGTCGATGCTAAATATGAGCg                                                                                                                                  >  1:51696/1‑141 (MQ=255)
      caCCGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCgg                                                                                                                                 >  2:1170453/1‑141 (MQ=255)
       cccGAATGCAATAAAATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGc                                                                                                                                <  1:346072/140‑1 (MQ=255)
        ccGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCa                                                                                                                               >  1:703054/1‑141 (MQ=255)
         cGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCAt                                                                                                                              >  1:650770/1‑141 (MQ=255)
         cGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCAt                                                                                                                              <  1:670103/141‑1 (MQ=255)
         cGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGGCTAACGTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCAt                                                                                                                              <  1:220441/141‑1 (MQ=255)
          gAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATa                                                                                                                             >  2:737739/1‑141 (MQ=255)
           aaTGCAATAATATTTTTTCAGTCAATATTTAAAAATTGTCTATCTTATACAATGCTTCCTGCTTTATTAATACCACGTTTTAAAATTTCATTCACAGGTAATGGTGCGCCCCAGCGTCGATTCTAAGTATGATCGGCATaa                                                                                                                            >  2:420779/1‑141 (MQ=255)
           aaTGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCGTTGCCCGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATaa                                                                                                                            >  2:420791/1‑141 (MQ=255)
            aTGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATaaa                                                                                                                           <  2:632707/141‑1 (MQ=255)
            aTGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATaaa                                                                                                                           >  1:553070/1‑141 (MQ=255)
            aTGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATaaa                                                                                                                           >  2:920987/1‑141 (MQ=255)
            aTGCAATAATATTATTTCGGTCAATATTCAACAATTGTCTAACTTGGTCAATGCCGCTTGCTTTCTTACTACCAAGTCTTACAATCTCACTGACAGGGAATCGTGCGCCCCAGCGGCGACACTAAATATGATCAGCATaaa                                                                                                                           >  2:323660/1‑141 (MQ=255)
             tGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATaaaa                                                                                                                          <  1:353569/141‑1 (MQ=255)
              gCAATAATATTGTTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATaaaaa                                                                                                                         <  2:1118582/141‑1 (MQ=255)
                     tattctTTCGCACAATATTCAAAAATTCTCTACCTTGCTCAACGCCGCTTGCTTTATTAATACCCAGTTTTACAATTTCAATGACAGGCAATGGAGCTCCCCAGCGTCGATGCTAAATAAGAAAGAAATAAAAATGAGTat                                                                                                                  <  2:203468/141‑2 (MQ=255)
                      attattTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAg                                                                                                                 >  1:155203/1‑141 (MQ=255)
                       ttattTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGc                                                                                                                >  2:964745/1‑141 (MQ=255)
                         attTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCat                                                                                                              >  2:940920/1‑141 (MQ=255)
                          tttCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCTAGTTTTACAGTTTCAATGACAGGGAATGGTTCGCCCCAGCGTCGATGCTAAATTTGATCGGCATAAAAATGAGTAAGCatt                                                                                                             >  2:13676/1‑140 (MQ=255)
                          tttCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATTACAGGGAATGTTGCGGCTCAGCGTCGATGCTAAATATGATCGGCATAAAATTGAGTAAGGatt                                                                                                             >  1:965452/1‑140 (MQ=255)
                          tttCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGCCTCGCCACCCCTTTCGTCGATGCTAAATATTATCGGCATAAAAATGAAGAAGCata                                                                                                             >  1:1172249/1‑141 (MQ=255)
                          tttCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCata                                                                                                             >  2:1137178/1‑141 (MQ=255)
                              ggTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATAtttt                                                                                                         >  2:829757/1‑141 (MQ=255)
                                 cAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCCAGAAGTACAAATGCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGa                                                                                                      <  1:274531/141‑1 (MQ=255)
                                    tatTTAAAAATTGTCTAGCTTGCTAAATGCCTCCTGCTTTATTCATAAAAAGTGTAACAATTTCAAAGCCAGGGACTGGTCCGCCCCAGCGTCAATGCCCAGTATGATCGGCATAAAAATGATAGAGCTTATGTTTGAttt                                                                                                   >  2:8465/1‑141 (MQ=255)
                                     atTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCGGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTc                                                                                                  <  1:758384/141‑1 (MQ=255)
                                     atTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTc                                                                                                  <  1:96547/141‑1 (MQ=255)
                                        tAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAgg                                                                                               >  1:1021651/1‑141 (MQ=255)
                                         aaaaaTTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGt                                                                                              >  2:52981/1‑141 (MQ=255)
                                              ttGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGGTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTAtttt                                                                                         <  1:489955/141‑1 (MQ=255)
                                               tGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTa                                                                                        >  1:376085/1‑141 (MQ=255)
                                               tGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTa                                                                                        >  2:208569/1‑141 (MQ=255)
                                                  tcaaCTTGCTCAATGCCTCTTGCTTTATTGAGACCAAGTTTTACACTTTCAATGACAGGAAATGGGGCGCCCCAGCGTGGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACtt                                                                                     <  2:484179/139‑1 (MQ=255)
                                                  cTACCTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACtt                                                                                     <  1:76514/141‑1 (MQ=255)
                                                  cTAACTTGCCTAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACtt                                                                                     <  2:828833/141‑1 (MQ=255)
                                                    aaCTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttc                                                                                   <  1:896697/141‑1 (MQ=255)
                                                    aaCTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttc                                                                                   <  1:1126783/141‑1 (MQ=255)
                                                       ttGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTttcttc                                                                                >  2:21127/1‑141 (MQ=255)
                                                           ccAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCt                                                                            <  1:344539/140‑1 (MQ=255)
                                                            cAATGCCTCTTGCTTTATTAGTACCAAGGTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATCTGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCtt                                                                           >  1:775096/1‑141 (MQ=255)
                                                               tGCCTCTTGCTTTATTAATAACAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCaa                                                                        >  2:576553/1‑141 (MQ=255)
                                                                tccTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCCTATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAt                                                                       <  2:437878/140‑1 (MQ=255)
                                                                gCCTCTTGCTTTATTAAGACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAt                                                                       <  1:694838/141‑1 (MQ=255)
                                                                gCCTCTTGCTTCTTTACTACCACGTTTTACTATTTCACTGACAGGGAATGGTGCGCCCCCGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAAt                                                                       <  1:575181/141‑1 (MQ=255)
                                                                  ctctTGCTTTATTAATACCAAGGTGTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTa                                                                     >  1:700701/1‑141 (MQ=255)
                                                                      tGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTaaaaa                                                                 >  2:1114433/1‑140 (MQ=255)
                                                                      tGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAAt                                                                 >  1:501925/1‑141 (MQ=255)
                                                                      tGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGAGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAAt                                                                 >  1:165462/1‑141 (MQ=255)
                                                                         tttATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAACTTGa                                                              >  2:1033885/1‑141 (MQ=255)
                                                                                   ccAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGAt                                                    <  1:940920/141‑1 (MQ=255)
                                                                                    cAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGATGTAGGGGAtt                                                   >  1:1134989/1‑141 (MQ=255)
                                                                                      aGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTATGGGATTCt                                                 >  1:966854/1‑141 (MQ=255)
                                                                                       gTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCtt                                                <  1:787073/141‑1 (MQ=255)
                                                                                          ttACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCa                                             <  2:981663/141‑1 (MQ=255)
                                                                                            ccAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAACATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAg                                           <  2:177244/140‑1 (MQ=255)
                                                                                            aCAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAGTTAAAATTGAGGTAGGGGATTCTTTCAAg                                           >  2:288137/1‑141 (MQ=255)
                                                                                            aCAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAg                                           >  2:937167/1‑141 (MQ=255)
                                                                                            aCAATTTCAATGACAGAGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAg                                           >  1:595241/1‑141 (MQ=255)
                                                                                              aaTTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCATGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGGAGGGAGTTCTTTCAAGTg                                         >  1:79766/1‑141 (MQ=255)
                                                                                                     aTGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGt                                  >  1:1085868/1‑141 (MQ=255)
                                                                                                     aTGACAGGGAATGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGt                                  <  1:737739/141‑1 (MQ=255)
                                                                                                         cAGGGAATGGTGCGCCCCAGCGTCGGTCAGGAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTAAGGTAGGGGATTCTTTCAAGTGGACAGGTAAAt                              >  2:547675/1‑141 (MQ=255)
                                                                                                              aaTGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTGTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCa                         >  2:342123/1‑141 (MQ=255)
                                                                                                                tGGTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGt                       >  2:1074221/1‑141 (MQ=255)
                                                                                                                  gTGCGCCCCAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGttt                     <  1:382405/141‑1 (MQ=255)
                                                                                                                  gTGCGCCCCAGCGTCGATGCTAAATATGATCCGCACAAAAATGAGTAAGCACATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGttt                     <  2:966508/141‑1 (MQ=255)
                                                                                                                     cgcCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGGGGACAAGTAAATTACCAGTTTGaa                  >  1:360894/1‑141 (MQ=255)
                                                                                                                       ccccAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAAtt                <  1:808549/141‑1 (MQ=255)
                                                                                                                         ccAGCGTCGATGCTAAATATGATCGGCATAAAAATGAGCAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCt              <  1:76193/141‑1 (MQ=255)
                                                                                                                         ccAGCGGCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCt              <  2:825225/141‑1 (MQ=255)
                                                                                                                            gCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGGAGGGGATTCTTTCAAGTGGACAAGTAAATGACCAGTTTGAATTCTTgg           >  1:895530/1‑141 (MQ=255)
                                                                                                                             cGTCGATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGa          <  1:21127/141‑1 (MQ=255)
                                                                                                                                 gATGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAGGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATTAc      >  1:299481/1‑141 (MQ=255)
                                                                                                                                  aTGCTAAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATTAcc     <  1:441953/141‑1 (MQ=255)
                                                                                                                                     cTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATTACCCAt  <  2:1021651/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
CTTCATCACCGAATGCAATAATATTATTTCGGTCAATATTTAAAAATTGTCTAACTTGCTCAATGCCTCTTGCTTTATTAATACCAAGTTTTACAATTTCAATGACAGGGAATGGTGCGCCCCAGCGTCGATGCTCAATATGATCGGCATAAAAATGAGTAAGCATATTTTTGATTTCAGGTATTTTACTTTCTTCGGCTTCAATTAAAATTGAGGTAGGGGATTCTTTCAAGTGGACAAGTAAATTACCAGTTTGAATTCTTGGATTACCCAT  >  USA300TCH1516_ALE/973078‑973351

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