Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 779,959 C→T 7.6% intergenic (+15/+57) csbB → / ← saeS Putative glycosyltransferase CsbB/Histidine protein kinase SaeS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE779,9590CT 7.6% 414.0 / 16.7 131intergenic (+15/+57)csbB/saeSPutative glycosyltransferase CsbB/Histidine protein kinase SaeS
Reads supporting (aligned to +/- strand):  ref base C (62/59);  new base T (4/6);  total (66/65)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 135 total aligned reads displayed.

ATATATTATTCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTA  >  USA300TCH1516_ALE/779869‑780047
                                                                                          |                                                                                        
atatatTATTCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGcc                                                                                          <  1:937240/91‑1 (MQ=255)
atatatTATTCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGcc                                                                                          <  1:3578286/91‑1 (MQ=255)
atatatTATTCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGcc                                                                                          >  1:603552/1‑91 (MQ=255)
atatatTATTCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGcc                                                                                          >  1:603584/1‑91 (MQ=255)
  atatTATTCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCt                                                                                        <  1:3167672/91‑1 (MQ=255)
   tattatTCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTc                                                                                       <  1:28210/91‑1 (MQ=255)
    attattCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCa                                                                                      <  1:4079346/91‑1 (MQ=255)
    attattCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCa                                                                                      <  1:2123515/91‑1 (MQ=255)
      tattCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCAtt                                                                                    <  1:53099/91‑1 (MQ=255)
      tattCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCAtt                                                                                    >  1:1149604/1‑91 (MQ=255)
      tattCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCCCAtt                                                                                    <  1:2903126/91‑1 (MQ=255)
      tattCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCCCAtt                                                                                    <  1:2903141/91‑1 (MQ=255)
       attCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTa                                                                                   >  1:564955/1‑91 (MQ=255)
       attCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTa                                                                                   >  1:564538/1‑91 (MQ=255)
            aGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATggg                                                                              >  1:38128/1‑91 (MQ=255)
             gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCTATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                             <  1:1538125/91‑1 (MQ=255)
             gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                             >  1:3871866/1‑91 (MQ=255)
             gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                             <  1:2015244/91‑1 (MQ=255)
             gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                             >  1:89381/1‑91 (MQ=255)
             gCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGa                                                                             <  1:2788662/91‑1 (MQ=255)
              cTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAg                                                                            <  1:2939254/91‑1 (MQ=255)
               tACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGc                                                                           <  1:2489437/91‑1 (MQ=255)
               tACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGc                                                                           <  1:366500/91‑1 (MQ=255)
               tACAAATCTCTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGc                                                                           <  1:3967846/91‑1 (MQ=255)
                aCAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGct                                                                          <  1:1858547/91‑1 (MQ=255)
                aCAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGct                                                                          >  1:526692/1‑91 (MQ=255)
                 cAAATCTTTCAAGTCTTGAACATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGctt                                                                         <  1:2522786/91‑1 (MQ=255)
                  aaaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGcttc                                                                        <  1:1111664/91‑1 (MQ=37)
                  aaaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGcttc                                                                        <  1:175468/91‑1 (MQ=37)
                   aaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttct                                                                       <  1:3832468/91‑1 (MQ=255)
                   aaTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttct                                                                       >  1:2363487/1‑91 (MQ=255)
                    aTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttctt                                                                      >  1:3625904/1‑91 (MQ=255)
                    aTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttctt                                                                      >  1:2910165/1‑91 (MQ=255)
                    aTCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGcttctt                                                                      >  1:2910690/1‑91 (MQ=255)
                     tCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGCTTCttt                                                                     <  1:1291955/91‑1 (MQ=37)
                      cTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCTCCCATTAATGAGGGCTTCtttt                                                                    <  1:3983295/91‑1 (MQ=37)
                          cAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGt                                                                <  1:2377976/91‑1 (MQ=255)
                            aGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCt                                                              >  1:1399328/1‑91 (MQ=255)
                            aGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCt                                                              >  1:1601521/1‑91 (MQ=255)
                            aGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCt                                                              >  1:1952273/1‑91 (MQ=255)
                             gTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCtt                                                             <  1:2515885/91‑1 (MQ=255)
                             gTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAAGAAACAAAAAAAGAAGCCCTCATTAAGGGGAGCTTCTTTTTAGTCtt                                                             <  1:3296817/91‑1 (MQ=255)
                                ttGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGc                                                          >  1:373344/1‑91 (MQ=255)
                                 tGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCa                                                         <  1:3566517/91‑1 (MQ=255)
                                 tGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCa                                                         <  1:1982861/91‑1 (MQ=255)
                                 tGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCa                                                         >  1:221720/1‑91 (MQ=255)
                                  gAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCat                                                        <  1:1853026/91‑1 (MQ=255)
                                   aaaaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGCAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatt                                                       >  1:75437/1‑91 (MQ=255)
                                   aaaaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatt                                                       <  1:3432025/91‑1 (MQ=255)
                                     aaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttt                                                     <  1:1670276/91‑1 (MQ=255)
                                     aaTGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttt                                                     >  1:1990662/1‑91 (MQ=255)
                                      atgatgAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCatttta                                                    >  1:3549546/1‑91 (MQ=255)
                                       tgatgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttat                                                   >  1:1552388/1‑91 (MQ=255)
                                       tgatgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttat                                                   <  1:583820/91‑1 (MQ=255)
                                        gatgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttatt                                                  >  1:362975/1‑91 (MQ=255)
                                         atgaGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCattttattt                                                 >  1:2841204/1‑91 (MQ=255)
                                           gagaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCAttttatttta                                               >  1:1137533/1‑91 (MQ=255)
                                           gagaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCAttttatttta                                               <  1:2491130/91‑1 (MQ=255)
                                           gagaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCAttttatttta                                               >  1:2040908/1‑91 (MQ=255)
                                            agaAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATtttattttat                                              >  1:658925/1‑91 (MQ=255)
                                                 gATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaata                                         >  1:527508/1‑91 (MQ=255)
                                                  aTACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTataaataa                                        <  1:843677/91‑1 (MQ=255)
                                                   tACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTAtaaataaa                                       <  1:2480988/91‑1 (MQ=255)
                                                    aCCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATaaataaat                                      <  1:1986578/91‑1 (MQ=255)
                                                    aCCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATaaataaat                                      <  1:584665/91‑1 (MQ=255)
                                                    aCCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATaaataaat                                      <  1:3932510/91‑1 (MQ=255)
                                                      ccATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCg                                    <  1:1492247/91‑1 (MQ=255)
                                                      ccATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCg                                    <  1:3121715/91‑1 (MQ=255)
                                                      ccATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCg                                    <  1:1581721/91‑1 (MQ=255)
                                                       cATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCgg                                   >  1:585906/1‑91 (MQ=255)
                                                             gTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATg                             <  1:856498/91‑1 (MQ=255)
                                                                tATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACg                          >  1:1641550/1‑91 (MQ=255)
                                                                  ttCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTa                        >  1:364147/1‑91 (MQ=255)
                                                                      aaataaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGt                    <  1:2912304/91‑1 (MQ=255)
                                                                       aataaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTc                   >  1:2060998/1‑91 (MQ=255)
                                                                       aataaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTc                   >  1:283296/1‑91 (MQ=255)
                                                                         taaaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTa                 >  1:363414/1‑91 (MQ=255)
                                                                           aaCAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAt               >  1:2175011/1‑91 (MQ=255)
                                                                             cAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAttt             >  1:2530888/1‑91 (MQ=255)
                                                                             cAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAttt             >  1:1963111/1‑91 (MQ=255)
                                                                             cAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAttt             >  1:2774028/1‑91 (MQ=255)
                                                                             cAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGCCTCTCCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAAttt             >  1:475733/1‑91 (MQ=255)
                                                                              aaaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtg            >  1:2927139/1‑91 (MQ=255)
                                                                               aaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgt           <  1:526886/91‑1 (MQ=255)
                                                                               aaaaaaGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTtgt           <  1:2526866/91‑1 (MQ=255)
                                                                                 aaaaGAAGCTCCCATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTgtgt         >  1:45642/1‑91 (MQ=37)
                                                                                 aaaaGAAGCTCCCATTAATGAGGGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTgtgt         >  1:2667524/1‑91 (MQ=37)
                                                                                    agaagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAAt      >  1:929841/1‑91 (MQ=255)
                                                                                     gaagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTACGACGTAATGTCTAATTTGTGTAATg     >  1:3847150/1‑91 (MQ=255)
                                                                                      aagCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGt    <  1:2331382/91‑1 (MQ=255)
                                                                                       agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt   >  1:2870800/1‑91 (MQ=255)
                                                                                       agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt   >  1:2925999/1‑91 (MQ=255)
                                                                                       agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt   >  1:702688/1‑91 (MQ=255)
                                                                                       agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt   >  1:685368/1‑91 (MQ=255)
                                                                                       agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt   >  1:2595843/1‑91 (MQ=255)
                                                                                       agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt   >  1:3705598/1‑91 (MQ=255)
                                                                                       agCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGtt   >  1:1478149/1‑91 (MQ=255)
                                                                                        gCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTa  <  1:1334562/91‑1 (MQ=255)
                                                                                        gCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTa  <  1:1335369/91‑1 (MQ=255)
                                                                                        gCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTa  >  1:3154343/1‑91 (MQ=255)
                                                                                          |                                                                                        
ATATATTATTCAAGCTACAAATCTTTCAAGTCTTGAAAATGATGAGAAGGATACCCATAAAGTTTATTCTAAATAAACAAAAAAAGAAGCCCTCATTAATGGGAGCTTCTTTTTAGTCTTTGCATTTTATTTTATAAATAAATCGGATTATGACGTAATGTCTAATTTGTGTAATGTTA  >  USA300TCH1516_ALE/779869‑780047

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