Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 79,612 A→T 8.4% intergenic (+221/‑39) USA300HOU_RS00355 → / → cmoB hypothetical protein/tRNA U34 carboxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE79,6120AT 8.4% 441.1 / 24.2 131intergenic (+221/‑39)USA300HOU_RS00355/cmoBhypothetical protein/tRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (52/68);  new base T (5/6);  total (57/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 140 total aligned reads displayed.

GTCTCGCGACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGAT  >  USA300TCH1516_ALE/79522‑79701
                                                                                          |                                                                                         
 tctcGCGACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAAc                                                                                          >  1:2301961/1‑91 (MQ=255)
     gcgACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTaa                                                                                      <  1:2059546/91‑1 (MQ=255)
       gACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAatt                                                                                    <  1:2782414/91‑1 (MQ=255)
        aCCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAataa                                                                                   <  1:133760/91‑3 (MQ=255)
        aCCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAatta                                                                                   <  1:3030167/91‑1 (MQ=255)
         ccACAGCTTAATTTTTGGTGTTCATTCATGGCTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattat                                                                                  <  1:1325281/91‑1 (MQ=255)
          cacaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattatt                                                                                 <  1:1564172/91‑1 (MQ=255)
          cacaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattatt                                                                                 <  1:1247675/91‑1 (MQ=255)
          cacaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattatt                                                                                 <  1:1288496/91‑1 (MQ=255)
          cacaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAattatt                                                                                 <  1:3296762/91‑1 (MQ=255)
           acaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAAttatta                                                                                >  1:3744646/1‑91 (MQ=255)
           acaGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAAttatta                                                                                >  1:3640174/1‑91 (MQ=255)
             aGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACg                                                                              <  1:2512902/91‑1 (MQ=255)
              gCTTAATTTTTGGTGTTCATTCCTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACgt                                                                             <  1:3672383/91‑2 (MQ=255)
              gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACgt                                                                             <  1:821319/91‑2 (MQ=255)
              gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACgt                                                                             <  1:3672315/91‑2 (MQ=255)
              gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                             >  1:3662658/1‑91 (MQ=255)
              gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                             >  1:2865986/1‑91 (MQ=255)
              gCTTAATTTTTGCTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                             >  1:3208476/1‑91 (MQ=255)
               cTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAt                                                                            <  1:2585559/91‑1 (MQ=255)
                   aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGa                                                                        >  1:2986162/1‑91 (MQ=255)
                   aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGa                                                                        >  1:3422554/1‑91 (MQ=255)
                   aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGa                                                                        >  1:817842/1‑91 (MQ=255)
                   aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGa                                                                        <  1:1096923/91‑1 (MQ=255)
                    tttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGaa                                                                       <  1:1372523/91‑1 (MQ=255)
                    tttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGaa                                                                       <  1:3773405/91‑1 (MQ=255)
                     ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAt                                                                      <  1:739500/91‑1 (MQ=38)
                     ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAt                                                                      <  1:2157259/91‑1 (MQ=255)
                       ttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAtag                                                                    <  1:717849/91‑3 (MQ=37)
                        tGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTaa                                                                   <  1:1834763/91‑1 (MQ=255)
                          gtgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGc                                                                 <  1:928613/91‑1 (MQ=255)
                          gtgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGc                                                                 <  1:929492/91‑1 (MQ=255)
                           tgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCg                                                                <  1:1017342/91‑1 (MQ=255)
                           tgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCg                                                                >  1:1733367/1‑91 (MQ=255)
                            gtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCga                                                               >  1:343625/1‑91 (MQ=255)
                             ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                              >  1:2224117/1‑91 (MQ=255)
                              tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagg                                                             <  1:1947033/91‑1 (MQ=255)
                              tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagg                                                             >  1:3874525/1‑91 (MQ=255)
                              tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagg                                                             <  1:1524226/91‑1 (MQ=255)
                              tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagg                                                             <  1:872869/91‑1 (MQ=255)
                                attcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgaggag                                                           <  1:2947518/91‑1 (MQ=255)
                                attcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgaggag                                                           <  1:2361620/91‑1 (MQ=255)
                                 ttcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGgaga                                                          >  1:2096469/1‑91 (MQ=255)
                                 ttcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGgaga                                                          >  1:2379939/1‑91 (MQ=255)
                                     tGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatga                                                      <  1:2952974/91‑1 (MQ=255)
                                      gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                     >  1:3941401/1‑91 (MQ=255)
                                      gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                     >  1:1074656/1‑91 (MQ=255)
                                      gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                     <  1:650855/91‑1 (MQ=255)
                                      gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                     <  1:3945075/91‑1 (MQ=255)
                                      gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                     >  1:1939048/1‑91 (MQ=255)
                                        cTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaat                                                   >  1:2738242/1‑91 (MQ=255)
                                         tttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatg                                                  >  1:2306472/1‑91 (MQ=255)
                                          ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                 <  1:69851/91‑1 (MQ=255)
                                          ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                 <  1:1017839/91‑1 (MQ=255)
                                          ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                 <  1:1483299/91‑1 (MQ=255)
                                          ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                 <  1:2124387/91‑1 (MQ=255)
                                          ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                 <  1:2124464/91‑1 (MQ=255)
                                            aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                               <  1:382668/91‑1 (MQ=255)
                                            aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAACGAAATGACt                                               >  1:3720470/1‑91 (MQ=255)
                                             aaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTa                                              >  1:1942450/1‑91 (MQ=255)
                                               aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                            >  1:3149392/1‑91 (MQ=255)
                                                  tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                         >  1:1346543/1‑91 (MQ=255)
                                                  tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                         >  1:505083/1‑91 (MQ=255)
                                                  tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                         >  1:4091975/1‑91 (MQ=255)
                                                  tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                         <  1:4083511/91‑1 (MQ=255)
                                                        ttGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGAc                                   <  1:1290270/91‑1 (MQ=255)
                                                        ttGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGAc                                   <  1:3236532/91‑1 (MQ=255)
                                                         tGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACg                                  <  1:1134531/91‑1 (MQ=255)
                                                           cATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGtt                                >  1:1062181/1‑91 (MQ=255)
                                                             tAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtg                              <  1:1225452/91‑1 (MQ=255)
                                                              aaaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtgg                             >  1:3746356/1‑91 (MQ=255)
                                                               aaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggt                            <  1:1303796/91‑1 (MQ=255)
                                                               aaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggt                            <  1:2692871/91‑1 (MQ=255)
                                                               aaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAAACTTATGACGTTtggt                            <  1:259938/91‑1 (MQ=255)
                                                                   tACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCa                        <  1:172919/91‑1 (MQ=255)
                                                                    aCATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCaa                       <  1:319387/91‑1 (MQ=255)
                                                                         gtgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAgg                  <  1:3080168/91‑1 (MQ=255)
                                                                         gtgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAgg                  <  1:69524/91‑1 (MQ=255)
                                                                         gtgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAgg                  <  1:1862287/91‑1 (MQ=255)
                                                                          tgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGt                 >  1:1720107/1‑91 (MQ=255)
                                                                           gtagtaCTATTCATAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                <  1:2500634/91‑1 (MQ=255)
                                                                           gtagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                >  1:2254280/1‑91 (MQ=255)
                                                                           gtagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                >  1:1440520/1‑91 (MQ=255)
                                                                           gtagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                <  1:613265/91‑1 (MQ=255)
                                                                           gtagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                <  1:2500211/91‑1 (MQ=255)
                                                                             agtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCaa              <  1:386136/91‑1 (MQ=255)
                                                                              gtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAg             >  1:3568421/1‑91 (MQ=255)
                                                                                aCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGaa           <  1:3916320/91‑1 (MQ=255)
                                                                                 cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt          >  1:3757110/1‑91 (MQ=255)
                                                                                 cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt          >  1:3756864/1‑91 (MQ=255)
                                                                                  tATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATc         >  1:3265093/1‑91 (MQ=255)
                                                                                     tCAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga      >  1:506817/1‑91 (MQ=38)
                                                                                     tCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGCCTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga      <  1:3811039/91‑1 (MQ=255)
                                                                                      cAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat     >  1:517279/1‑91 (MQ=255)
                                                                                      cAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat     >  1:2771319/1‑91 (MQ=39)
                                                                                      cAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat     <  1:3637718/91‑1 (MQ=255)
                                                                                        aaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatga   >  1:2248792/1‑91 (MQ=255)
                                                                                        aaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatga   <  1:2939930/91‑1 (MQ=255)
                                                                                         aaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgat  <  1:73481/91‑1 (MQ=255)
                                                                                         aaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgat  <  1:1268415/91‑1 (MQ=255)
                                                                                          |                                                                                         
GTCTCGCGACCACAGCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGAT  >  USA300TCH1516_ALE/79522‑79701

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