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,5080CA98.9% 293.0 / NA 94S338Y (TCT→TAT) ftsZCell division protein FtsZ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (39/54);  total (40/54)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAG  >  USA300TCH1516_ALE/1196373‑1196643
                                                                                                                                       |                                                                                                                                       
ttttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTGGGAACAAGCGTAAATAATTATAGCa                                                                                                                                    >  2:627396/1‑141 (MQ=255)
  ttCGGTACAGTTATTAATCCTGAATTACCAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAAt                                                                                                                                  <  1:742306/141‑1 (MQ=255)
  ttCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAACCGTAAATACTTATAGCAAt                                                                                                                                  <  1:742826/141‑1 (MQ=255)
   tCGGTACAGTTATTAATCCTGACTTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATg                                                                                                                                 >  2:725714/1‑141 (MQ=255)
    cGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACCTCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                                                                                <  2:244494/141‑1 (MQ=255)
    cGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                                                                                <  1:578082/141‑1 (MQ=255)
    cGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                                                                                <  2:943713/141‑1 (MQ=255)
    cGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGATTCGGAACAAGCGTAAATACTTATAGCAATGCa                                                                                                                               >  2:621296/1‑141 (MQ=255)
     ggTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCa                                                                                                                               >  1:136367/1‑141 (MQ=255)
      gTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCaa                                                                                                                              >  1:48877/1‑141 (MQ=255)
      gTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCaa                                                                                                                              >  2:908146/1‑141 (MQ=255)
       tACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAAc                                                                                                                             >  1:248081/1‑141 (MQ=255)
       tACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAAc                                                                                                                             >  1:168956/1‑141 (MQ=255)
       tACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAAc                                                                                                                             >  1:155770/1‑141 (MQ=255)
         cacttattaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACtt                                                                                                                           <  1:155572/138‑1 (MQ=255)
         cAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACtt                                                                                                                           <  2:461979/141‑1 (MQ=255)
         aaGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACtt                                                                                                                           <  2:659603/140‑1 (MQ=255)
               ttaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAAGGCAACTTCTAAAg                                                                                                                     >  1:619022/1‑141 (MQ=255)
               ttaATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGAGTCGGAACAAGCGTAAATACTTATAGCAATGCAGCTTCTAAAg                                                                                                                     >  1:803779/1‑141 (MQ=255)
               ttaATCCTGAATTACAAGATGAGATTGTTGTAAAAGTTATTACAACTGGTTTTGATGACAAACAAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTAAAGCAATGCAACTCCCAAAg                                                                                                                     >  2:409635/1‑141 (MQ=255)
                    ccTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGaa                                                                                                                >  1:166149/1‑141 (MQ=255)
                    ccTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGaa                                                                                                                >  1:706088/1‑141 (MQ=255)
                        aaTTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCAt                                                                                                            <  1:812019/141‑1 (MQ=255)
                        aaTTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCAt                                                                                                            <  1:551829/141‑1 (MQ=255)
                            aCAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                                                                        <  2:173627/141‑1 (MQ=255)
                              aaGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACtt                                                                                                      >  1:211512/1‑141 (MQ=255)
                              aaGATGAGATCGTTTTAACAGTTACTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACtt                                                                                                      <  2:261893/141‑1 (MQ=255)
                                  tgcgaTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCaaa                                                                                                  <  1:206840/138‑1 (MQ=255)
                                  taagaTTGTTGTAACAGTTATTGCTACTGGTTTTGATGACCAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCTTTAATACTTATAGCAATGCAACTTCTACAGATGAATCATTCACTTCaaa                                                                                                  <  1:221561/139‑1 (MQ=255)
                                  ggagaTTGTTGTAACAGTTATTGCAACTGGTTTGGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGAGTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCaaa                                                                                                  <  1:685099/140‑1 (MQ=255)
                                     gaTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGGACAAGCGTAAATCCTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtc                                                                                               >  1:212472/1‑141 (MQ=255)
                                      aTTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtca                                                                                              >  1:83587/1‑141 (MQ=255)
                                        tgttgtAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtcatc                                                                                            <  1:559613/141‑1 (MQ=255)
                                               aCAGTTATTGCAACTGGTTTTGATGACAAACCCCCTGCACAGGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGca                                                                                     >  1:983876/1‑141 (MQ=255)
                                               aCAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGca                                                                                     >  2:355455/1‑141 (MQ=255)
                                                   ttATTGCAACTGGTTTTGATGACAAACAAAAATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATAATTCACTTCACATTCATCAAATGCACAAg                                                                                 >  1:275449/1‑141 (MQ=255)
                                                     aTTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCa                                                                               <  2:569485/141‑1 (MQ=255)
                                                       tGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                             <  1:333472/141‑1 (MQ=255)
                                                       tGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                             <  1:834792/141‑1 (MQ=255)
                                                       tGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                             <  1:683018/141‑1 (MQ=255)
                                                       tGCAACTGGTTTTGATGACAAACCAACATCACATGGCCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAACACTTATAGCAATGCAACTTCTAAAGATGACTCATTCACTTCAACTTCATCAAATTCACCAgacac                                                                             >  2:468201/1‑137 (MQ=255)
                                                       tGCAACTGGTTTTGAGGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                                                                             <  1:365849/141‑1 (MQ=255)
                                                          aaCTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGa                                                                          <  1:247649/141‑1 (MQ=255)
                                                          aaCTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGa                                                                          >  2:606349/1‑141 (MQ=255)
                                                           aCTGGTTTTGATGACAACCCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGAt                                                                         <  2:985036/141‑1 (MQ=255)
                                                           aCTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGAt                                                                         <  1:505625/141‑1 (MQ=255)
                                                           aCTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGAt                                                                         <  2:428854/141‑1 (MQ=255)
                                                           aCTGGTGTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGAt                                                                         <  2:125459/141‑1 (MQ=255)
                                                              ggTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgt                                                                      >  2:592195/1‑141 (MQ=255)
                                                              ggTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTCGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCCACTGCACAAGCAACTGATAgt                                                                      >  2:400612/1‑141 (MQ=255)
                                                              ggTTTTGATGACAAACCAACATCACATGGTCGTAAATCGGGTAGCACAGGATTCGGAACAAGCGGAAATACTTATAGCAATGCAACTTCTAACGATGAATGATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgt                                                                      >  2:289883/1‑141 (MQ=255)
                                                                   tgatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCCCAAGCAACTGATAGTGTAAg                                                                 >  2:182441/1‑141 (MQ=255)
                                                                   tgatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAg                                                                 >  2:288716/1‑141 (MQ=255)
                                                                    gctgccAAACCAACAACACATCGTCGAAAAACCTGTAGCACCGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGt                                                                <  1:972546/136‑1 (MQ=255)
                                                                    gatgaCCAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGt                                                                <  1:93445/141‑1 (MQ=255)
                                                                    gatgaCAAACCAACATCATATGGTCGTAAAGCTGGTAGCAATGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGt                                                                <  1:942312/141‑1 (MQ=255)
                                                                    gatgaCAAACCAACAGCCGAGGGTCGCAAATCCGTTAGCACTGTATTCGGAACAAGCGTACATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATATTGTAAGt                                                                <  2:371258/141‑1 (MQ=255)
                                                                    catgaCAAACCAACATCACATGGTCGTAAATCTAGCACCACCGAATTCTGAACAACCGTAAATACTTATAACAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGt                                                                <  2:647610/140‑1 (MQ=255)
                                                                       gaCAAACCAACATCACATGGTCGTAAATCTGGTTGCACGGGATTCGGAGCAAGCGTAAATACATATAGCAATGCAACTTCTAAAGATGAATCATTACCTTCAAATTCATCAAATGCACAAGCAACTCATGGTGTAAGTGaa                                                             >  2:336291/1‑141 (MQ=255)
                                                                       gaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGaa                                                             >  2:112509/1‑141 (MQ=255)
                                                                           aaCCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGaa                                                         <  2:118079/141‑1 (MQ=255)
                                                                             ccAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAaca                                                       >  2:394527/1‑141 (MQ=255)
                                                                                    cacaTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACaa                                                <  1:857660/141‑1 (MQ=255)
                                                                                      caTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACt                                              <  2:211512/141‑1 (MQ=255)
                                                                                       aTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTa                                             <  2:83587/141‑1 (MQ=255)
                                                                                              tcaaTCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAagaaga                                      <  2:821685/139‑1 (MQ=255)
                                                                                              tAAATCTGGTAGCACTGGATTCGGAAAAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATGCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAagaaga                                      <  2:201662/141‑1 (MQ=255)
                                                                                              gaaaTTTGGTAGAACTGGATTCGGAACAATCGTAAATACTTATAGCAATGCAACTTCTAAAAATGAATCATTCACTTCAAATTCATCAAATTCACAAGCAACTGAGAGTGTAAGTGAAAAAACATATACAACTAaagaaaa                                      <  2:282147/140‑3 (MQ=255)
                                                                                                 aTCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAgatga                                   <  2:922162/141‑1 (MQ=255)
                                                                                                 aTCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAgatga                                   <  1:592195/141‑1 (MQ=255)
                                                                                                    tGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGatat                                >  2:755431/1‑141 (MQ=255)
                                                                                                      gTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGCACACATACAAACAAAAAAGATGAAATTc                              >  1:851705/1‑141 (MQ=255)
                                                                                                      gTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTc                              >  1:808650/1‑141 (MQ=255)
                                                                                                      gTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTc                              <  1:792466/141‑1 (MQ=255)
                                                                                                      gTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTc                              >  2:224555/1‑141 (MQ=255)
                                                                                                      gTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGAAACTGATAGTGTAAGTGAAAGAACACATAAAAATAAAGAAGATGATATTc                              >  2:42438/1‑141 (MQ=255)
                                                                                                      gTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGAAATACACCTTCTCAAGATGAGTCATTCACTTCAAATCAATAAAATGCACAAGCAACTTATATTTTAAGTTAAAGAACACATACAACTAAAAAAGATGATCTTc                              >  2:577145/1‑141 (MQ=255)
                                                                                                      gTAGCACTGGATTCGGAACAAGAGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTc                              <  1:337599/141‑1 (MQ=255)
                                                                                                      gTAGCACTGGATTCGGAACAAGACTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTc                              <  1:910055/141‑1 (MQ=255)
                                                                                                       tAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTcc                             <  2:848409/141‑1 (MQ=255)
                                                                                                       tAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTcc                             <  1:49622/141‑1 (MQ=255)
                                                                                                       tAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTcc                             <  1:908146/141‑1 (MQ=255)
                                                                                                         gCACTGGATTCGGAACAAGCGGAAATACTTATAGCAATGTAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAGGATTATATTCCTa                           <  2:190895/141‑1 (MQ=255)
                                                                                                             tGGCTTCGGAACAAGCATAAATACTTATAGCAATGCAACTTCTAAAGAGGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCtt                       <  1:971671/141‑1 (MQ=255)
                                                                                                             tGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCtt                       >  1:639220/1‑141 (MQ=255)
                                                                                                                  tCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTa                  <  1:334172/141‑1 (MQ=255)
                                                                                                                      aaCAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGaaa              <  1:355455/141‑1 (MQ=255)
                                                                                                                      aaCAAACGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGaaa              >  2:415560/1‑141 (MQ=255)
                                                                                                                           gCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATagag         <  1:353698/141‑1 (MQ=255)
                                                                                                                           gCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATagag         >  1:353689/1‑141 (MQ=255)
                                                                                                                           gCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATagag         >  1:893172/1‑141 (MQ=255)
                                                                                                                                 aTACTTATCGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGaagaaa   <  2:212472/141‑1 (MQ=255)
                                                                                                                                 aTACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGaagaaa   <  2:48877/141‑1 (MQ=255)
                                                                                                                                  aaCTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGaagaaag  <  2:929501/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TTTTCGGTACAGTTATTAATCCTGAATTACAAGATGAGATTGTTGTAACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATACAACTAAAGAAGATGATATTCCTAGCTTCATTAGAAATAGAGAAGAAAG  >  USA300TCH1516_ALE/1196373‑1196643

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