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,5080CA94.5% 286.0 / 7.6 92S338Y (TCT→TAT) ftsZCell division protein FtsZ
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base A (37/49);  total (39/53)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGT  >  USA300TCH1516_ALE/1196372‑1196647
                                                                                                                                        |                                                                                                                                           
aTTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGc                                                                                                                                         >  4:15763/1‑141 (MQ=255)
aTTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGc                                                                                                                                         >  2:1003738/1‑141 (MQ=255)
  tttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCaa                                                                                                                                       <  1:211134/141‑1 (MQ=255)
  tttCGGGACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCaa                                                                                                                                       <  2:580157/141‑1 (MQ=255)
    tCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAAAATCACATGGACGTAATAGAGGTAGCACTGGATTCGAAACAAGCGTAAATACTTATAGCAATg                                                                                                                                     >  1:208415/1‑141 (MQ=255)
     cGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                                                                                    >  3:183194/1‑141 (MQ=255)
     cGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                                                                                    >  2:15767/1‑141 (MQ=255)
      ggTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCa                                                                                                                                   >  2:314511/1‑141 (MQ=255)
      ggTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCa                                                                                                                                   <  1:638421/141‑1 (MQ=255)
        tACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAAc                                                                                                                                 <  1:211602/141‑1 (MQ=255)
         aCAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACt                                                                                                                                <  3:159989/141‑1 (MQ=255)
            gTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCt                                                                                                                             <  1:653952/141‑1 (MQ=255)
               attaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTaaa                                                                                                                          <  1:31132/141‑1 (MQ=255)
                      cTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAAt                                                                                                                   >  1:339469/1‑141 (MQ=255)
                        gAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACCTCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCa                                                                                                                 >  2:615303/1‑141 (MQ=255)
                        gAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCa                                                                                                                 >  2:130525/1‑141 (MQ=255)
                        gAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCa                                                                                                                 >  1:471229/1‑141 (MQ=255)
                        gAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCa                                                                                                                 >  4:54385/1‑141 (MQ=255)
                             aCAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                                                                            <  1:600873/141‑1 (MQ=255)
                             aCAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                                                                            <  2:215375/141‑1 (MQ=255)
                             aCAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTAATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                                                                            <  2:373177/141‑1 (MQ=255)
                               aaGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACtt                                                                                                          >  1:767031/1‑141 (MQ=255)
                                aGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTc                                                                                                         <  3:215204/141‑1 (MQ=255)
                                aGATGAGATTGTTGTAACAGTTATTGCAACTGGTGTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTc                                                                                                         <  1:954603/141‑1 (MQ=255)
                                   tGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCaaa                                                                                                      <  1:381966/141‑1 (MQ=255)
                                    gagaTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAt                                                                                                     >  1:519304/1‑141 (MQ=255)
                                    gagaTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAt                                                                                                     <  1:320483/141‑1 (MQ=255)
                                     agaTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAtt                                                                                                    <  2:705328/141‑1 (MQ=255)
                                        ttgttgTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtcat                                                                                                 >  1:410104/1‑141 (MQ=255)
                                         tgttgtAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtcatc                                                                                                <  2:223454/141‑1 (MQ=255)
                                          gttgtAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtcatca                                                                                               <  3:15763/141‑1 (MQ=255)
                                           ttgtAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCaa                                                                                              <  2:666252/141‑1 (MQ=255)
                                            tgtAACAGTTATTGCAACCGGTTTTGATGACAAACCAACGGCACATGGTCGTCAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCaaa                                                                                             <  3:41183/141‑1 (MQ=255)
                                              tAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATg                                                                                           >  2:339046/1‑141 (MQ=255)
                                                aCAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGca                                                                                         >  2:228216/1‑141 (MQ=255)
                                                    ttATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAg                                                                                     <  4:203684/141‑1 (MQ=255)
                                                    ttATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAg                                                                                     <  1:313179/141‑1 (MQ=255)
                                                       ttGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtc                                                                                                   >  1:838110/1‑124 (MQ=255)
                                                       ttGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtc                                                                                                   <  2:838110/124‑1 (MQ=255)
                                                        tGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                                 >  1:185438/1‑141 (MQ=255)
                                                            aCTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGAt                                                                             <  1:246939/141‑1 (MQ=255)
                                                                gTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtg                                                                         <  3:419005/141‑1 (MQ=255)
                                                                gTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtg                                                                         <  1:121599/141‑1 (MQ=255)
                                                                gTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtg                                                                         <  3:241081/141‑1 (MQ=255)
                                                                gTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtg                                                                         <  1:350056/141‑1 (MQ=255)
                                                                   ttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTaa                                                                      <  2:301223/141‑1 (MQ=255)
                                                                    tgatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAg                                                                     >  1:86029/1‑141 (MQ=255)
                                                                    tgatgaCAAACAAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACATCTAAAGATGAATCATTCACTTCAAAATCATCAAATGCACAAGCAACTGATAGTGTAAg                                                                     <  4:41840/141‑1 (MQ=255)
                                                                     gatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGt                                                                    <  1:570408/141‑1 (MQ=255)
                                                                     gatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGt                                                                    <  4:346093/141‑1 (MQ=255)
                                                                     gatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGt                                                                    <  2:689059/141‑1 (MQ=255)
                                                                          cAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAg                                                               >  1:792469/1‑141 (MQ=255)
                                                                              ccAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAaca                                                           <  1:130525/141‑1 (MQ=255)
                                                                                 aCATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACAt                                                        >  4:1757/1‑141 (MQ=255)
                                                                                 aCATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACAt                                                        >  4:117380/1‑141 (MQ=255)
                                                                                    tCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACa                                                     >  2:652430/1‑141 (MQ=255)
                                                                                    tCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACa                                                     >  2:492/1‑141 (MQ=255)
                                                                                    tCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACa                                                     >  2:579111/1‑141 (MQ=255)
                                                                                      acaTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAAc                                                   >  1:141479/1‑141 (MQ=255)
                                                                                       aaTGGTCGTAAATCTGGGAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACt                                                  <  2:682803/140‑1 (MQ=255)
                                                                                         tGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTaa                                                >  3:281171/1‑141 (MQ=255)
                                                                                          ggTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTaaa                                               <  1:41056/141‑1 (MQ=255)
                                                                                              gTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAaagaag                                           >  2:976424/1‑141 (MQ=255)
                                                                                               tAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAagaaga                                          <  2:519304/141‑1 (MQ=255)
                                                                                               tAAATCTGGTAGCACTGGAGTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGAGGAAGCATTCACTTCAAACTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAagaaga                                          <  4:302090/141‑1 (MQ=255)
                                                                                                  aTCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAgatga                                       <  2:290772/141‑1 (MQ=255)
                                                                                                   tCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAgatgat                                      <  1:1003879/141‑1 (MQ=255)
                                                                                                    cTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACTACTAAAGAAGATGata                                     >  4:222071/1‑141 (MQ=255)
                                                                                                    cTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGata                                     >  3:109628/1‑141 (MQ=255)
                                                                                                     tGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGatat                                    <  3:4398/141‑1 (MQ=255)
                                                                                                      ggTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATAtt                                   >  4:356213/1‑141 (MQ=255)
                                                                                                       gTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGa                                                                                                                        <  4:181333/55‑1 (MQ=255)
                                                                                                       gTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGa                                                                                                                        >  3:181333/1‑55 (MQ=255)
                                                                                                        tAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTcc                                 <  3:50279/141‑1 (MQ=255)
                                                                                                         aGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCt                                <  4:175084/141‑1 (MQ=255)
                                                                                                         aGCACTGGATTCGGAACAAGCCGAAATACTTGGATCAACGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAACTGCACAAGCAACTGTTAGTGTAATTGAAAGAACACATACAACTAAAGAAAAGGATATTCCt                                >  4:145394/1‑141 (MQ=255)
                                                                                                           cACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAGCTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAg                              <  1:495630/141‑1 (MQ=255)
                                                                                                           cACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAg                              <  2:21846/141‑1 (MQ=255)
                                                                                                             cTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCt                            <  2:767031/141‑1 (MQ=255)
                                                                                                                   tCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTa                      >  4:56493/1‑141 (MQ=255)
                                                                                                                      gAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGaa                   >  1:108829/1‑141 (MQ=255)
                                                                                                                       aaCAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGaaa                  <  1:623517/141‑1 (MQ=255)
                                                                                                                         cAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATa                <  1:832453/141‑1 (MQ=255)
                                                                                                                         cAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATa                <  1:466737/141‑1 (MQ=255)
                                                                                                                          aaGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATag               >  1:383445/1‑141 (MQ=255)
                                                                                                                           aGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATaga              >  1:421496/1‑141 (MQ=255)
                                                                                                                           aGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGAAACTGATAGTGTAAGTGAAAGAACAAATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATaga              >  4:167423/1‑141 (MQ=255)
                                                                                                                           aGCGTAAATACTGATAGCAATGCAACTTCGAAAGATGAATCATGCACGTCTAATTCAGCAAATGCACAAGCAACTGATAGTGGAAGTGAAAGAACACATACAACCAAAGAAGATGATATTCACAGCTTCATTAGAAATaga              >  4:84883/1‑141 (MQ=255)
                                                                                                                                aaaTACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGagaaga         <  4:82320/141‑1 (MQ=255)
                                                                                                                                  aTACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGaagaaa       <  4:254003/141‑1 (MQ=255)
                                                                                                                                   tACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACGCATACAACTAAAGAAGATGATATTCCAAGCTTCATTAGAAATAGAGaagaaag      <  2:637613/141‑1 (MQ=255)
                                                                                                                                       tCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGt  <  4:361273/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                           
ATTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAGACGT  >  USA300TCH1516_ALE/1196372‑1196647

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