Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,196,508 C→A S338Y (TCT→TAT)  ftsZ → Cell division protein FtsZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,196,5080CA99.1% 333.6 / NA 109S338Y (TCT→TAT) ftsZCell division protein FtsZ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (52/56);  total (53/56)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

TTAATATGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGTTC  >  USA300TCH1516_ALE/1196364‑1196649
                                                                                                                                                  |                                                                                                                                             
ttAATATGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAg                                                                                                                                            >  1:284423/1‑148 (MQ=255)
    tatGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAAt                                                                                                                                        >  2:1383752/1‑148 (MQ=255)
     atGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCTTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATAATTATAGCAATg                                                                                                                                       >  2:556214/1‑148 (MQ=255)
     atGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATg                                                                                                                                       >  2:699902/1‑148 (MQ=255)
      tGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                                                                                      >  1:714982/1‑148 (MQ=255)
      tGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                                                                                      >  2:141179/1‑148 (MQ=255)
      tGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                                                                                      >  1:524143/1‑148 (MQ=255)
       gATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCCCTGGATTCCGCCCAAGCGTAAAGACTTATAGCAATGCa                                                                                                                                     >  1:27763/1‑148 (MQ=255)
       gATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAGATACTTATAGGAATGCa                                                                                                                                     >  2:727863/1‑148 (MQ=255)
        aTTTTCGGTACAGTTATTAATCATGAATTACAAGATGAGATTGTTGTAACAGTTATAGCAACTGGTT‑‑TTGATGAGAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCaa                                                                                                                                    >  2:623801/1‑148 (MQ=255)
         ttttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAAc                                                                                                                                   >  1:766611/1‑148 (MQ=255)
         ttttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAAc                                                                                                                                   >  2:1324968/1‑148 (MQ=255)
         ttttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGGAAATATGGTAGCGCTGGATTCGGAACAAGCGTAAATACTTATAGCAATGGAAc                                                                                                                                   >  2:498516/1‑148 (MQ=255)
          tttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACt                                                                                                                                  <  2:1247764/148‑1 (MQ=255)
          tttCGGTACAGTTATTAAGCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACt                                                                                                                                  <  2:877079/148‑1 (MQ=255)
              ggTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTTACGGATCCTGC‑‑CTGGTTTGTTGGTGGCGAACCACCATCACATGGTGGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTa                                                                                                                              <  2:1353820/148‑1 (MQ=255)
              ggTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTa                                                                                                                              >  1:1280337/1‑148 (MQ=255)
               gTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTaa                                                                                                                             <  1:865655/148‑1 (MQ=255)
                    gTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATg                                                                                                                        <  1:776685/148‑1 (MQ=255)
                     ttattaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGa                                                                                                                       <  1:1354477/148‑1 (MQ=255)
                      tattaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGaa                                                                                                                      >  1:278961/1‑148 (MQ=255)
                      tattaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGaa                                                                                                                      >  2:1342646/1‑148 (MQ=255)
                        ttaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAACAATGCAACTTCTAAAGATGAATc                                                                                                                    >  1:1109174/1‑148 (MQ=255)
                        ttaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACGTATAGCAATGCAACTTCTAAAGATGAATc                                                                                                                    >  2:894803/1‑148 (MQ=255)
                          aaTCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCAt                                                                                                                  >  2:986869/1‑148 (MQ=255)
                          aaTCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCCGAACAAGCGCAAATACTTCTAGCAATGCAACTTCTAAAGATGaaccct                                                                                                                  >  1:913921/1‑144 (MQ=255)
                           aTCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGCCAAACCAACATCACATGGTCGTAAATCTGGTGGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCAtt                                                                                                                 >  2:1132934/1‑148 (MQ=255)
                           aTCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCAtt                                                                                                                 >  2:1231665/1‑148 (MQ=255)
                           aTCCTGAATTACAAGATGAGATAGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCAtt                                                                                                                 >  1:754021/1‑148 (MQ=255)
                             ccTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCAttac                                                                                                               >  1:313099/1‑146 (MQ=255)
                             ccTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCa                                                                                                               >  2:1086959/1‑148 (MQ=255)
                             ccTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCa                                                                                                               >  1:1100115/1‑148 (MQ=255)
                             ccTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCa                                                                                                               >  1:1099991/1‑148 (MQ=255)
                                     aCAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAt                                                                                                       <  1:1213551/148‑1 (MQ=255)
                                       aaGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtc                                                                                                     >  1:1403705/1‑148 (MQ=255)
                                        aGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtca                                                                                                    >  2:1160540/1‑148 (MQ=255)
                                        aGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTACTAAAGATGAATCATTCACTTCAAATtca                                                                                                    >  2:1038226/1‑148 (MQ=255)
                                            gagaTTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCaa                                                                                                >  1:1123364/1‑148 (MQ=255)
                                               aTTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCTTTCACTTCAAATTCATCAAATg                                                                                             <  2:955332/148‑1 (MQ=255)
                                               aTTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATg                                                                                             <  2:451870/148‑1 (MQ=255)
                                                   ttgtAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtca                                                                                                    <  1:1298191/137‑1 (MQ=255)
                                                   ttgtAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGAAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtca                                                                                                    >  2:1298191/1‑137 (MQ=255)
                                                        aCAGTTATTGCTACTGGTT‑‑TTGATGACAAACCGACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTTAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCCCTTCAAATTCATCAAATGCACAAGCaa                                                                                    <  2:787287/148‑1 (MQ=255)
                                                        aCAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGCACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACCAGCac                                                                                    >  1:1323962/1‑147 (MQ=255)
                                                        aCAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATATCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCaa                                                                                    <  2:1110289/148‑1 (MQ=255)
                                                         cAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                                   >  2:1182883/1‑148 (MQ=255)
                                                         cAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                                   >  2:961696/1‑148 (MQ=255)
                                                         cAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                                   >  2:341108/1‑148 (MQ=255)
                                                         cAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                                   >  2:324722/1‑148 (MQ=255)
                                                          aGTTAGTGCAACTGGTT‑‑TTGATGCCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACATCAAATTCATCAAATGCACAAGCAACt                                                                                  <  1:212523/148‑1 (MQ=255)
                                                               ttGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGAAACTGATAg                                                                             <  1:1342646/148‑1 (MQ=255)
                                                                 gccaCTGGTT‑‑TTGATTACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGTTTCGGCACAAGCGTAAATTCTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtg                                                                           <  2:500400/145‑1 (MQ=255)
                                                                       ggTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCTCAAGCAACTGATAGTCTAAGTg                                                                     <  2:956524/148‑1 (MQ=255)
                                                                       ggTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTg                                                                     <  2:474553/148‑1 (MQ=255)
                                                                         tt‑‑ttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGaa                                                                   >  1:1368146/1‑148 (MQ=255)
                                                                               gctgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCAGAAGCAACTGATAGTGTAAGTGAAAGaa                                                               <  2:1343309/146‑1 (MQ=255)
                                                                               gatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGaa                                                               <  1:560210/148‑1 (MQ=255)
                                                                                      aaCCAACATCACATGGTCGTAACTCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATAc                                                        >  1:779335/1‑148 (MQ=255)
                                                                                      aaCCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGCTGACTCATTCACTTCAAGTTCATCAAATGCACAAGCAACTGATAGTGTAATTTAAAAAACACATat                                                        >  1:29352/1‑147 (MQ=255)
                                                                                      aaCCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGGAAAAGCAACTCATAGTGTAAGGGAAAGAACACATAc                                                        >  2:384923/1‑148 (MQ=255)
                                                                                      aaCCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAAGACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATAc                                                        >  1:619098/1‑148 (MQ=255)
                                                                                       aCCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACa                                                       <  1:1230993/148‑1 (MQ=255)
                                                                                       aCCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACa                                                       <  1:228370/148‑1 (MQ=255)
                                                                                        ccAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACaa                                                      <  1:653298/148‑1 (MQ=255)
                                                                                        ccAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACaa                                                      <  2:278961/148‑1 (MQ=255)
                                                                                         cAACATCACATGGTCGTAAATCTGGTAGCACTAGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATAAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAAc                                                     <  2:785826/148‑1 (MQ=255)
                                                                                          aaCATGACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACt                                                    <  2:913921/148‑1 (MQ=255)
                                                                                           aCATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTa                                                   <  1:986869/148‑1 (MQ=255)
                                                                                           aCATCAAATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTa                                                   <  1:1131772/148‑1 (MQ=255)
                                                                                             aTCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTaaa                                                 <  2:766611/148‑1 (MQ=255)
                                                                                             aTCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTaaa                                                 <  1:1324968/148‑1 (MQ=255)
                                                                                              tCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAAATGAAACATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAACGAACACATACAACTaaaa                                                >  2:1207106/1‑147 (MQ=255)
                                                                                              tCACATGGTCGTAAAGCTGGTAGCACTGGACTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAaag                                                <  2:77407/148‑1 (MQ=255)
                                                                                                  aTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACCCATACAACTAAagaaga                                            >  1:1205712/1‑148 (MQ=255)
                                                                                                  aTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAagaaga                                            >  1:551301/1‑148 (MQ=255)
                                                                                                  aTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAagaaga                                            >  1:1205747/1‑148 (MQ=255)
                                                                                                   tGGTCGTAAATCTTGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAgat                                           <  1:473924/148‑1 (MQ=255)
                                                                                                    ggTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAGCACATACAACTAAAGAAgatg                                          <  2:1087769/148‑1 (MQ=255)
                                                                                                     gTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAgatga                                         <  2:1175682/148‑1 (MQ=255)
                                                                                                         tAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGCACACGTACAACTAAAGAAGATGATAtt                                     <  1:775254/148‑1 (MQ=255)
                                                                                                            aTCTGGTAGCACTGGATTCGGAACCAGCGCAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCt                                  >  2:1001686/1‑148 (MQ=255)
                                                                                                            aTCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCt                                  >  2:938239/1‑148 (MQ=255)
                                                                                                               tGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGc                               >  1:711776/1‑148 (MQ=255)
                                                                                                                 gTAGCACTGTATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAACGAGGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGGTAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCtt                             <  2:584734/148‑1 (MQ=255)
                                                                                                                 gTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCtt                             <  2:246627/148‑1 (MQ=255)
                                                                                                                   aGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCa                           >  1:1246276/1‑148 (MQ=255)
                                                                                                                        tGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGa                      <  2:1108675/148‑1 (MQ=255)
                                                                                                                        tGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGa                      <  2:902173/148‑1 (MQ=255)
                                                                                                                        tGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGa                      <  2:902213/148‑1 (MQ=255)
                                                                                                                        tGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATCAATCATTCACTTGAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGa                      <  2:963803/148‑1 (MQ=255)
                                                                                                                         ggATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTCCATTAGaa                     <  2:739407/148‑1 (MQ=255)
                                                                                                                          gATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAGGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGaaa                    <  2:754021/148‑1 (MQ=255)
                                                                                                                           aTTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAAt                   <  1:699902/148‑1 (MQ=255)
                                                                                                                           aTTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAAt                   <  1:556214/148‑1 (MQ=255)
                                                                                                                            ttCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATa                  >  1:1261152/1‑148 (MQ=255)
                                                                                                                              cGGAACCAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATaga                <  2:407971/148‑1 (MQ=255)
                                                                                                                                gAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTAAGTTACAGCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAgaga              <  1:324722/148‑1 (MQ=255)
                                                                                                                                     aGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGaagaaa         <  2:284423/148‑1 (MQ=255)
                                                                                                                                        gTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGAc      <  2:278099/148‑1 (MQ=255)
                                                                                                                                            aTACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGTTc  >  1:1360906/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TTAATATGATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTT‑‑TTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGTTC  >  USA300TCH1516_ALE/1196364‑1196649

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