Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 536,868 T→A 10.0% intergenic (+61/‑250) rplY → / → pth 50S ribosomal protein L25/Peptidyl‑tRNA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE536,8680TA10.0% 726.9 / 52.8 227intergenic (+61/‑250)rplY/pth50S ribosomal protein L25/Peptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base T (105/99);  new base A (11/12);  total (116/111)
Fisher's exact test for biased strand distribution p-value = 8.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 243 total aligned reads displayed.

AAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGAC  >  USA300TCH1516_ALE/536778‑536958
                                                                                          |                                                                                          
  agaagaCGGAGAAAAAACTGAAGAGTAATATTAACCTGTTACATTAAAGAGTTTATACTTTGCTTAACCAGCCCTGTGCTTATTCTAacat                                                                                          >  2:677071/1‑91 (MQ=255)
   gaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAAtata                                                                                         >  2:2119120/1‑91 (MQ=255)
   gaagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAAtata                                                                                         <  1:2679706/91‑1 (MQ=255)
    aagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATaa                                                                                        <  2:734536/91‑1 (MQ=255)
    aagaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATaa                                                                                        <  2:3626628/91‑1 (MQ=255)
     agaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAg                                                                                       >  2:2487889/1‑91 (MQ=255)
     agaCGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAg                                                                                       <  2:2915724/91‑1 (MQ=38)
           agaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGt                                                                                 >  1:160435/1‑91 (MQ=255)
            gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                                <  1:471207/91‑1 (MQ=255)
            gaaAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTg                                                                                <  1:1623441/91‑1 (MQ=255)
             aaaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc                                                                               <  2:3703195/91‑1 (MQ=255)
             aaaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGc                                                                               >  2:955786/1‑91 (MQ=255)
              aaaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCt                                                                              <  1:1891931/91‑1 (MQ=255)
               aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                             >  2:2035301/1‑91 (MQ=255)
               aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCtt                                                                             >  2:2065200/1‑91 (MQ=255)
               aaaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCACAGTGCtt                                                                             <  2:2204140/91‑1 (MQ=21)
                aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                            <  1:2948012/91‑1 (MQ=255)
                aaaCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTGAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttt                                                                            <  2:3186025/91‑1 (MQ=255)
                  aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                          >  2:1649562/1‑91 (MQ=255)
                  aCTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCttttt                                                                          <  2:1550054/91‑1 (MQ=255)
                   cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtg                                                                         <  2:728902/91‑1 (MQ=255)
                   cTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACCATGCTTATATTAAAATAAGCACAGTGCttgtta                                                                         <  2:3751939/91‑5 (MQ=18)
                    tGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgt                                                                        >  2:2018183/1‑91 (MQ=255)
                    tGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATCGTGCTTTTtgt                                                                        >  1:3518111/1‑91 (MQ=255)
                     gAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgtg                                                                       >  2:2144201/1‑91 (MQ=255)
                     gAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTtgtg                                                                       >  2:256031/1‑91 (MQ=255)
                      aaGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                      >  1:932267/1‑91 (MQ=255)
                      aaGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTgtgt                                                                      >  1:2197994/1‑91 (MQ=255)
                       agagTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGtt                                                                     <  2:2244191/91‑1 (MQ=255)
                          gTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatt                                                                  >  1:638671/1‑91 (MQ=255)
                           tAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatta                                                                 <  1:3557345/91‑1 (MQ=255)
                             aTTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATtata                                                               <  1:2018183/91‑1 (MQ=255)
                               tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                             <  1:1085037/91‑1 (MQ=255)
                               tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                             >  2:2594824/1‑91 (MQ=255)
                               tttAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATaaa                                                             >  1:3704909/1‑91 (MQ=255)
                                 tAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGc                                                           <  1:1551907/91‑1 (MQ=255)
                                 tAATCTGTTACATTAAAGTTTTTATACTTTGTGTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGc                                                           <  1:3313872/91‑1 (MQ=255)
                                 tAATCAGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGc                                                           <  2:3576746/91‑1 (MQ=255)
                                  aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGttatt                                                                  <  1:2628867/83‑1 (MQ=255)
                                  aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                          <  1:2850635/91‑1 (MQ=255)
                                  aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                          >  2:1435378/1‑91 (MQ=255)
                                  aaTCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCt                                                          <  1:795136/91‑1 (MQ=255)
                                    tCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCtta                                                        >  2:1590832/1‑91 (MQ=255)
                                      tGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaat                                                      >  1:1917202/1‑91 (MQ=255)
                                       gTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGTttaatt                                                     <  2:3800139/91‑1 (MQ=255)
                                        ttACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaatta                                                    >  1:2515096/1‑91 (MQ=255)
                                        ttACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaatta                                                    >  1:2515773/1‑91 (MQ=255)
                                        ttACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaatta                                                    <  1:1648984/91‑1 (MQ=255)
                                         tACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaattaa                                                   >  2:648975/1‑91 (MQ=255)
                                         tACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCttaattaa                                                   >  2:1223387/1‑91 (MQ=255)
                                           cATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAAc                                                 >  2:1466626/1‑91 (MQ=255)
                                             ttAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACtt                                               <  2:345735/91‑1 (MQ=255)
                                              tAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACttt                                              <  1:460211/91‑1 (MQ=255)
                                              tAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACttt                                              >  1:3141576/1‑91 (MQ=255)
                                                aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                            >  2:2338824/1‑91 (MQ=255)
                                                aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                            <  1:2026453/91‑1 (MQ=255)
                                                aaGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttat                                            >  2:405459/1‑91 (MQ=255)
                                                 aGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatt                                           <  1:2872594/91‑1 (MQ=255)
                                                  gTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTttatta                                          <  1:1050197/91‑1 (MQ=255)
                                                   tttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                                   >  1:1195792/1‑49 (MQ=255)
                                                   tttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATg                                                                                   <  2:1195792/49‑1 (MQ=255)
                                                    ttttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACt                                        >  1:1907928/1‑91 (MQ=255)
                                                     tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                       <  1:955786/91‑1 (MQ=255)
                                                     tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                       <  2:3209058/91‑1 (MQ=255)
                                                     tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                       >  2:2863888/1‑91 (MQ=255)
                                                     tttATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACtt                                       >  2:2863606/1‑91 (MQ=255)
                                                        ataCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGt                                    <  2:1731610/91‑1 (MQ=255)
                                                         taCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTa                                   <  1:2088422/91‑1 (MQ=255)
                                                         taCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTa                                   <  1:3681623/91‑1 (MQ=255)
                                                          aCTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTAc                                  <  2:3096476/91‑1 (MQ=255)
                                                           cTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                                 >  1:311250/1‑91 (MQ=255)
                                                           cTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                                 >  1:1453327/1‑91 (MQ=255)
                                                           cTTTGGTTAACAAGCACTGTGCTTATGTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACt                                 <  1:3766062/91‑1 (MQ=255)
                                                            tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                                >  1:418465/1‑91 (MQ=255)
                                                            tttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTa                                <  2:2671489/91‑1 (MQ=255)
                                                             ttGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTaa                               >  1:3492174/1‑91 (MQ=255)
                                                                 ttAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGtt                           <  2:2646996/91‑1 (MQ=255)
                                                                  tAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGttt                          <  1:2594824/91‑1 (MQ=255)
                                                                    aaaaGCACCATGCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                        >  2:795440/3‑91 (MQ=21)
                                                                    aCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                        <  2:2772206/91‑1 (MQ=255)
                                                                    aCAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaa                        <  1:1432297/91‑1 (MQ=255)
                                                                    aCAAGCACTGTGCTTAGTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATGAAACTTTATTACTTTGTACGAAAGTgtaa                        <  2:3384984/91‑1 (MQ=255)
                                                                     cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                       >  1:2577730/1‑91 (MQ=255)
                                                                     cAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaat                       <  1:660099/91‑1 (MQ=255)
                                                                      aaGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatt                      <  2:1390381/91‑1 (MQ=255)
                                                                       aGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTttaatta                     >  2:2550121/1‑91 (MQ=255)
                                                                         cacagtgCTTGTTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTtaattaat                   >  1:1678178/5‑91 (MQ=255)
                                                                         cACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTACAGTTtaattaat                   <  1:2986159/91‑1 (MQ=255)
                                                                           ccatgCTTATATTAAAATAAGCACAGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAttt                 >  1:1890819/4‑91 (MQ=37)
                                                                            tgtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAAtttt                >  2:243229/1‑91 (MQ=255)
                                                                             gtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTa               >  1:2693309/1‑91 (MQ=255)
                                                                             gtgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTa               >  2:2372175/1‑91 (MQ=255)
                                                                              tgCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAg              >  1:901108/1‑91 (MQ=255)
                                                                                  tATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAg          >  1:2761704/1‑91 (MQ=255)
                                                                                   aTTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGt         <  1:775002/91‑1 (MQ=255)
                                                                                    ttttAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTa        <  2:442591/91‑1 (MQ=255)
                                                                                       tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa     >  2:494564/1‑91 (MQ=255)
                                                                                       tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa     <  1:3806007/91‑1 (MQ=255)
                                                                                       tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa     >  1:3691537/1‑91 (MQ=255)
                                                                                       tAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTaaaa     <  1:2171068/91‑1 (MQ=255)
                                                                                          |                                                                                          
AAAGAAGACGAAGAAAAAACTGAAGAGTAATTTTAATCTGTTACATTAAAGTTTTTATACTTTGTTTAACAAGCACTGTGCTTATTTTAATATAAGCATGGTGCTTTTTGTGTTATTATAAAGCTTAATTAAACTTTATTACTTTGTACTAAAGTTTAATTAATTTTAGTGAGTAAAAGAC  >  USA300TCH1516_ALE/536778‑536958

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: