Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 79,626 A→T 10.6% intergenic (+235/‑25) 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,6260AT10.6% 700.3 / 54.6 225intergenic (+235/‑25)USA300HOU_RS00355/cmoBhypothetical protein/tRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (99/102);  new base T (16/8);  total (115/110)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 230 total aligned reads displayed.

GCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCACGAGAC  >  USA300TCH1516_ALE/79536‑79716
                                                                                          |                                                                                          
gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                                            <  1:3214195/91‑1 (MQ=255)
gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                                            >  2:3488554/1‑91 (MQ=255)
 cTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAt                                                                                           <  1:1918724/91‑1 (MQ=255)
  ttAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAtt                                                                                          <  1:1547243/91‑1 (MQ=255)
    aaTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTg                                                                                        >  2:2524918/1‑91 (MQ=39)
    aaTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTAAGAATTg                                                                                        >  2:2961679/1‑91 (MQ=255)
     aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGa                                                                                       <  2:3437537/91‑1 (MQ=38)
     aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGa                                                                                       >  1:654140/1‑91 (MQ=255)
      tttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGaa                                                                                      <  1:835402/91‑1 (MQ=38)
       ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAt                                                                                     <  1:2524918/91‑1 (MQ=38)
       ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAt                                                                                     <  1:3247334/91‑1 (MQ=38)
       ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAt                                                                                     <  1:3248763/91‑1 (MQ=38)
       ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAt                                                                                     >  2:1900015/1‑91 (MQ=255)
       ttttGGTGTTCATTAATGAAAATAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAt                                                                                     <  2:2948486/91‑1 (MQ=255)
         ttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTa                                                                                   >  1:440955/1‑91 (MQ=255)
         tgggTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTa                                                                                   >  1:175717/3‑91 (MQ=255)
          tGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTaa                                                                                  >  1:2314095/1‑91 (MQ=255)
          tGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTaa                                                                                  >  2:2906916/1‑91 (MQ=255)
             tgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCg                                                                               >  2:1266854/1‑91 (MQ=255)
              gtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCga                                                                              >  1:3349696/1‑91 (MQ=255)
               ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                                             >  1:3318879/1‑91 (MQ=255)
               ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                                             <  1:2140362/91‑1 (MQ=255)
               ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                                             <  2:362270/91‑1 (MQ=255)
                tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATAAAGCgagg                                                                            <  2:507449/91‑1 (MQ=255)
                 cattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagga                                                                           <  2:1399690/91‑1 (MQ=255)
                 cattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagga                                                                           <  2:618589/91‑1 (MQ=255)
                  attcctGACTTTAAAAAATCCTTATTGCATAAATTTACCTAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgaggag                                                                          >  2:3872096/1‑91 (MQ=255)
                   ttcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGgaga                                                                         <  2:1920644/91‑1 (MQ=255)
                   ttcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGgaga                                                                         >  1:1888088/1‑91 (MQ=255)
                    tcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaa                                                                        <  1:2093893/91‑1 (MQ=255)
                    tcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaa                                                                        <  2:2707728/91‑1 (MQ=255)
                     catGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaat                                                                       <  1:1427853/91‑1 (MQ=255)
                       tGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatga                                                                     >  1:2643991/1‑91 (MQ=255)
                        gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                                    >  2:2662574/1‑91 (MQ=255)
                        gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                                    >  2:1129608/1‑91 (MQ=255)
                        gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                                    >  1:1664991/1‑91 (MQ=255)
                          cTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaat                                                                  <  2:784949/91‑1 (MQ=255)
                            ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGTAGaatgaaatga                                                                <  1:2974757/91‑1 (MQ=255)
                            ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgcaatga                                                                >  2:3568775/1‑91 (MQ=255)
                            ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                                <  2:3754426/91‑1 (MQ=255)
                            ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                                <  1:3756757/91‑1 (MQ=255)
                             tAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGAc                                                               <  2:1147594/91‑1 (MQ=255)
                              aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                                              <  1:3140758/91‑1 (MQ=255)
                              aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                                              <  1:1074913/91‑1 (MQ=255)
                                aaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTaa                                                            <  1:1005179/91‑1 (MQ=255)
                                 aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                                           >  2:160838/1‑91 (MQ=255)
                                 aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                                           >  2:3276113/1‑91 (MQ=255)
                                 aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                                           >  2:358713/1‑91 (MQ=255)
                                  aaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGa                                                          <  1:2846743/91‑1 (MQ=255)
                                    tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                                        >  2:3175005/1‑91 (MQ=255)
                                    tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                                        <  2:1456635/91‑1 (MQ=255)
                                     ccTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACtt                                                       >  2:2144798/1‑91 (MQ=255)
                                      cTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTa                                                      <  1:797298/91‑1 (MQ=255)
                                        tATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATg                                                    <  1:3705238/91‑1 (MQ=255)
                                        tATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATg                                                    >  1:1778701/1‑91 (MQ=255)
                                         aTTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGa                                                   >  2:183249/1‑91 (MQ=255)
                                         aTTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGa                                                   <  1:3543802/91‑1 (MQ=255)
                                          ttGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGAc                                                  >  1:2440017/1‑91 (MQ=255)
                                           tGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACg                                                 <  2:766525/91‑1 (MQ=255)
                                              aTAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGttt                                              <  1:991031/91‑1 (MQ=255)
                                               tAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtg                                             >  2:419879/1‑91 (MQ=255)
                                                     tACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCa                                       >  1:1686593/1‑91 (MQ=255)
                                                      aCATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCaa                                      >  1:1994523/1‑91 (MQ=255)
                                                       cATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCaaa                                     >  1:3692787/1‑91 (MQ=255)
                                                           gtgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAgg                                 >  1:3088060/1‑91 (MQ=255)
                                                            tgtAGTAATATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGt                                <  2:3457478/91‑1 (MQ=255)
                                                               agtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCaa                             <  2:3088060/91‑1 (MQ=255)
                                                                 taCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGa                           <  2:1548970/91‑1 (MQ=255)
                                                                  aCTATTCAAGACGTAATTATGACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGaa                          <  2:635175/91‑1 (MQ=255)
                                                                  aCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGaa                          >  1:1330787/1‑91 (MQ=255)
                                                                   cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt                         >  1:2866141/1‑91 (MQ=255)
                                                                   cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt                         >  1:2868591/1‑91 (MQ=255)
                                                                     aTTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCa                       <  1:1177998/91‑1 (MQ=255)
                                                                       tCAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga                     >  2:3130023/1‑91 (MQ=38)
                                                                       tCAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga                     >  1:943868/1‑91 (MQ=38)
                                                                       tCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga                     >  1:2307832/1‑91 (MQ=255)
                                                                       tCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga                     <  2:1651952/91‑1 (MQ=255)
                                                                       tCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga                     >  2:649894/1‑91 (MQ=255)
                                                                        cAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat                    >  2:606377/1‑91 (MQ=38)
                                                                        cAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat                    >  2:3284788/1‑91 (MQ=38)
                                                                         aaatCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAaatg                   >  1:2257364/1‑91 (MQ=38)
                                                                         aaaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatg                   >  2:295146/1‑91 (MQ=255)
                                                                         aaaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatg                   <  2:2800319/91‑1 (MQ=255)
                                                                            aCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatg                <  1:1782932/91‑1 (MQ=255)
                                                                             cGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatga               >  2:692179/1‑91 (MQ=255)
                                                                              gTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatgat              >  2:146614/1‑91 (MQ=255)
                                                                               taataaTTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGata             >  2:439597/1‑91 (MQ=255)
                                                                                aaTTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat            <  1:3338101/91‑1 (MQ=255)
                                                                                aaTTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat            <  1:2359473/91‑1 (MQ=255)
                                                                                aaTTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat            >  2:3117885/1‑91 (MQ=255)
                                                                                  ttattaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATgg          >  1:1842257/1‑91 (MQ=255)
                                                                                  ttattaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATgg          >  2:2417584/1‑91 (MQ=255)
                                                                                  ttattaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATgg          >  2:121775/1‑91 (MQ=255)
                                                                                  taattaCGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATgg          >  1:2568442/3‑91 (MQ=255)
                                                                                   tattaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGc         >  2:36637/1‑91 (MQ=255)
                                                                                    atcaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGAAGAAAAAAAAAAAAAAGGTCAAGAATCAGATGATGATATGGCa        <  2:210054/88‑1 (MQ=255)
                                                                                     ttaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCAc       <  1:780700/91‑1 (MQ=255)
                                                                                     ttaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCAc       >  1:1431329/1‑91 (MQ=255)
                                                                                        cGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCACgag    <  1:1212324/91‑1 (MQ=255)
                                                                                          aTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCACGAGAc  >  2:3584662/1‑91 (MQ=255)
                                                                                          |                                                                                          
GCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCACGAGAC  >  USA300TCH1516_ALE/79536‑79716

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