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.2% 466.8 / ‑4.2 131S338Y (TCT→TAT) ftsZCell division protein FtsZ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (68/62);  total (69/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 133 total aligned reads displayed.

ACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATA  >  USA300TCH1516_ALE/1196420‑1196594
                                                                                        |                                                                                      
aCAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTtata                                                                                      >  2:2299168/1‑91 (MQ=255)
 cAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAg                                                                                     >  2:2216501/1‑91 (MQ=255)
 cAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAg                                                                                     <  2:658995/91‑1 (MQ=255)
  aGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGc                                                                                    >  1:2826079/1‑91 (MQ=255)
    ttATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCaa                                                                                  >  2:696138/1‑91 (MQ=255)
      aTTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAAta                                                                                >  1:986970/1‑90 (MQ=255)
      aTTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATg                                                                                >  1:1987169/1‑91 (MQ=255)
      aTTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATg                                                                                <  1:1235451/91‑1 (MQ=255)
       ttGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                               <  2:727012/91‑1 (MQ=255)
       ttGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGc                                                                               >  1:2065357/1‑91 (MQ=255)
        tGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCa                                                                              >  1:2062161/1‑91 (MQ=255)
        tGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCa                                                                              <  2:183541/91‑1 (MQ=255)
         gCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCaa                                                                             >  1:1324039/1‑91 (MQ=255)
          cAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAAc                                                                            <  1:749486/91‑1 (MQ=255)
            aCTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACtt                                                                          >  1:533743/1‑91 (MQ=255)
             cTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTc                                                                         >  1:937372/1‑91 (MQ=255)
             cTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTc                                                                         >  2:822968/1‑91 (MQ=255)
               ggTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTa                                                                       <  2:1613916/91‑1 (MQ=255)
                 ttttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTaaa                                                                     >  2:2837400/1‑91 (MQ=255)
                 ttttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTaaa                                                                     >  2:2837224/1‑91 (MQ=255)
                  tttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAg                                                                    >  2:1784769/1‑91 (MQ=255)
                   ttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGa                                                                   <  1:302213/91‑1 (MQ=255)
                   ttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGa                                                                   >  1:1653986/1‑91 (MQ=255)
                   ttGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGa                                                                   >  1:1403652/1‑91 (MQ=255)
                    tgatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACGAGCGTAAATACTGATAGCAATGCATCTTCTAAAGAt                                                                  >  2:2368368/1‑91 (MQ=255)
                     gatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATg                                                                 <  2:2741600/91‑1 (MQ=255)
                     gatgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATg                                                                 <  1:2832801/91‑1 (MQ=255)
                     aatgaCAAACCATCATCACATTGTCGTAAATCTGGTAGCACTGGATTCAGAACAAGCGTAAATTCTTATAGCAATGCAACTTCTAAAGATg                                                                 <  1:2012490/90‑1 (MQ=255)
                       tgaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGaa                                                               >  2:1246983/1‑91 (MQ=255)
                        gaCAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAAt                                                              >  1:166024/1‑91 (MQ=255)
                          cAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCa                                                            >  2:288317/1‑91 (MQ=255)
                            aaCCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCAtt                                                          <  1:2844449/91‑1 (MQ=255)
                              ccAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCa                                                        >  1:2419014/1‑91 (MQ=255)
                              ccAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCa                                                        >  1:3008695/1‑91 (MQ=255)
                               cAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                       <  2:1480168/91‑1 (MQ=255)
                               cAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                       >  2:1183821/1‑91 (MQ=255)
                               cAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                       >  2:312458/1‑91 (MQ=255)
                               cAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                       <  1:496294/91‑1 (MQ=255)
                               cAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                       <  1:1147862/91‑1 (MQ=255)
                               cAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCAc                                                       <  1:2304286/91‑1 (MQ=255)
                                aaCATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCCCAAATACTTATAGCAATGCAACTGCTAAAGCTCAACcactcact                                                      >  2:186777/1‑91 (MQ=255)
                                 aCATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACtt                                                     >  2:2252870/1‑91 (MQ=255)
                                 aCATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACtt                                                     >  1:1313267/1‑91 (MQ=255)
                                 aCATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACtt                                                     >  1:1313349/1‑91 (MQ=255)
                                 aCATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACtt                                                     >  2:1362941/1‑91 (MQ=255)
                                 aCATCACATGGACGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACtt                                                     <  2:2363064/91‑1 (MQ=255)
                                 aCATCACATGGACGTAAATCTGGGAGCACTGGATTCGGAACAAGCGTAAATACTTATAGGAATGCAACTTCTAAAGATGAATCATTCACtt                                                     <  2:2362940/91‑1 (MQ=255)
                                 aCATCACAAGGCCGTAAATCTGGTCGCATTGGATTCGGTACAAGCGTAATTACTTATAGCAATGCAACTTCTAAAGATGAATCATTCCCtt                                                     >  2:2418275/1‑91 (MQ=255)
                                  cATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTc                                                    <  2:761376/91‑1 (MQ=255)
                                  cATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTc                                                    <  2:1122116/91‑1 (MQ=255)
                                  cATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTc                                                    >  2:1967600/1‑91 (MQ=255)
                                    tCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCaa                                                  >  2:2929942/1‑91 (MQ=255)
                                    tCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCaa                                                  >  1:1859869/1‑91 (MQ=255)
                                     cacaTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCaaa                                                 <  1:2768443/91‑1 (MQ=255)
                                      acaTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAt                                                >  1:1656897/1‑91 (MQ=255)
                                       caTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAtt                                               >  2:2705076/1‑91 (MQ=255)
                                       caTGGTCGTAAATCTGGTACCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAtt                                               >  2:2707207/1‑91 (MQ=255)
                                        aTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAAtgc                                              <  2:2099699/91‑3 (MQ=255)
                                        aTGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtc                                              >  2:2822115/1‑91 (MQ=255)
                                         tGGTCGAAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtca                                             <  2:237521/91‑1 (MQ=255)
                                          ggTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtcat                                            >  1:1532926/1‑91 (MQ=255)
                                          ggTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtcat                                            <  1:2060556/91‑1 (MQ=255)
                                          ggTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATtcat                                            >  1:1462866/1‑91 (MQ=255)
                                          ggTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTACAAATtcat                                            >  1:1463223/1‑91 (MQ=255)
                                             cGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCaa                                         <  2:1832249/91‑1 (MQ=255)
                                              gTAAATCTGGTAGCACTGGATTCGGATCAAGCGTAAATACTTATAGCAATGCAACGTCTAAAGATGAATCATTCACTTCAAATTCATCaaa                                        <  2:964197/91‑1 (MQ=255)
                                                aaaTCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCCAATg                                      >  1:2573926/1‑91 (MQ=255)
                                                aaaTCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATg                                      >  1:1332459/1‑91 (MQ=255)
                                                aaaTCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATg                                      <  2:1060029/91‑1 (MQ=255)
                                                     tGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACaa                                 >  1:453247/1‑91 (MQ=255)
                                                     tGGTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACaa                                 >  1:2245254/1‑91 (MQ=255)
                                                      ggTAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAg                                >  1:1247237/1‑91 (MQ=255)
                                                        tAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCa                              <  1:393489/91‑1 (MQ=255)
                                                        tAGCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCa                              >  1:831937/1‑91 (MQ=255)
                                                          gCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                            <  2:72080/91‑1 (MQ=255)
                                                          gCACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAAc                            <  2:482061/91‑1 (MQ=255)
                                                           cACTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACt                           >  1:2781350/1‑91 (MQ=255)
                                                            aCTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTg                          <  1:3011969/91‑1 (MQ=255)
                                                            aCTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTg                          <  1:486807/91‑1 (MQ=255)
                                                             cTGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGa                         <  2:1081039/91‑1 (MQ=255)
                                                              tGGATTCGGAATAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGAt                        <  1:2809655/91‑1 (MQ=255)
                                                              tGGATTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGAt                        <  2:316769/91‑1 (MQ=255)
                                                                 aTTCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgt                     <  1:1313518/91‑1 (MQ=255)
                                                                   tCGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAgtgt                   >  1:527651/1‑91 (MQ=255)
                                                                    cGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTa                  <  1:1925538/91‑1 (MQ=255)
                                                                    cGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTa                  >  2:2505404/1‑91 (MQ=255)
                                                                    cGGAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTa                  >  1:2015376/1‑91 (MQ=255)
                                                                     ggAACAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTaa                 <  1:2424387/91‑1 (MQ=255)
                                                                      gAACTAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAg                <  2:1547908/91‑1 (MQ=255)
                                                                       aaCAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGt               <  1:1102107/91‑1 (MQ=255)
                                                                        aCAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTg              <  2:1003987/91‑1 (MQ=255)
                                                                         cAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGa             <  1:682353/91‑1 (MQ=255)
                                                                         cAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGa             >  2:168041/1‑91 (MQ=255)
                                                                         cAAGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGa             >  1:911934/1‑91 (MQ=255)
                                                                          aaGCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGaa            <  1:2115109/91‑1 (MQ=255)
                                                                            gCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAg          <  1:1837032/91‑1 (MQ=255)
                                                                            gCGTAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAg          >  2:930270/1‑91 (MQ=255)
                                                                               tAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCACATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAac       <  2:1382346/91‑1 (MQ=255)
                                                                               tAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAac       >  2:2993963/1‑91 (MQ=255)
                                                                               tAAATACTTATAGCAATGCAACTTCTAAAGATGAATCATCCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAac       >  2:934063/1‑91 (MQ=255)
                                                                                        |                                                                                      
ACAGTTATTGCAACTGGTTTTGATGACAAACCAACATCACATGGTCGTAAATCTGGTAGCACTGGATTCGGAACAAGCGTAAATACTTCTAGCAATGCAACTTCTAAAGATGAATCATTCACTTCAAATTCATCAAATGCACAAGCAACTGATAGTGTAAGTGAAAGAACACATA  >  USA300TCH1516_ALE/1196420‑1196594

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: