Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,779,807 A→G 15.5% intergenic (‑113/+29) USA300TCH1516_01664 ← / ← rplT hypothetical protein/50S ribosomal protein L20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,779,8070AG15.5% 329.2 / 48.3 129intergenic (‑113/+29)USA300TCH1516_01664/rplThypothetical protein/50S ribosomal protein L20
Reads supporting (aligned to +/- strand):  ref base A (60/49);  new base G (13/7);  total (73/56)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 131 total aligned reads displayed.

TGTTCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTG  >  USA300TCH1516_ALE/1779721‑1779897
                                                                                      |                                                                                          
tgttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAAt                                                                                        >  1:2027158/1‑91 (MQ=255)
  ttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGc                                                                                      >  1:2835969/1‑91 (MQ=255)
  ttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGc                                                                                      <  1:1983988/91‑1 (MQ=255)
  ttCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGc                                                                                      <  1:1983744/91‑1 (MQ=255)
   tctcTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGcc                                                                                     <  1:944578/91‑1 (MQ=255)
   tctcTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGcc                                                                                     <  1:1359257/91‑1 (MQ=255)
       ttAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCAtct                                                                                 <  1:1027760/91‑1 (MQ=38)
        tAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCAtctc                                                                                <  1:386121/91‑1 (MQ=38)
        tAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAtctc                                                                                <  1:832496/91‑1 (MQ=255)
        tAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAtctc                                                                                >  1:340045/1‑91 (MQ=255)
        tAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCAtccc                                                                                >  1:339501/1‑89 (MQ=255)
         aaaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATctct                                                                               >  1:2300713/1‑91 (MQ=255)
         aaaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATctct                                                                               >  1:2300730/1‑91 (MQ=255)
         aaaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATctct                                                                               >  1:244526/1‑91 (MQ=255)
         aaaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATctct                                                                               <  1:2031576/91‑1 (MQ=255)
          aaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCtt                                                                              <  1:3720384/91‑1 (MQ=38)
          aaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCtt                                                                              <  1:3252313/91‑1 (MQ=38)
          aaTGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCtt                                                                              <  1:3252674/91‑1 (MQ=38)
             gCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCttttt                                                                           <  1:2930998/91‑1 (MQ=38)
             gCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCttttt                                                                           <  1:1362257/91‑1 (MQ=255)
              cATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCtttttt                                                                          <  1:1509434/91‑1 (MQ=255)
              cATTTTAGCAAATATGTCATATTTTAGCGATTAAATGAGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCtttttt                                                                          >  1:3512239/1‑91 (MQ=255)
               aTTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTta                                                                         >  1:2327423/1‑91 (MQ=255)
               aTTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTta                                                                         >  1:1465890/1‑91 (MQ=255)
                ttttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTtat                                                                        >  1:1333055/1‑91 (MQ=255)
                ttttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTtat                                                                        >  1:4002986/1‑91 (MQ=255)
                 tttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGGCATTGTCATCTCTTTTTtatt                                                                       <  1:2796268/91‑2 (MQ=37)
                 tttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTtata                                                                       <  1:3100860/91‑1 (MQ=255)
                 tttAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTtata                                                                       >  1:526746/1‑91 (MQ=255)
                    aGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgt                                                                    >  1:1634870/1‑91 (MQ=255)
                     gCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAgtg                                                                   <  1:1464141/91‑1 (MQ=255)
                       aaaTATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGtt                                                                 <  1:2832999/91‑1 (MQ=255)
                       aaaTATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGtt                                                                 >  1:869254/1‑91 (MQ=255)
                        aaTATGTCATAGTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGtta                                                                <  1:3765917/91‑1 (MQ=255)
                         atatGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttat                                                               >  1:2629452/1‑91 (MQ=255)
                          tatGTCATATTTTAGCGATTAAATGTTGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGttatt                                                              <  1:1828317/91‑1 (MQ=255)
                               cATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTCTTAtttt                                                         >  1:3912051/1‑91 (MQ=255)
                               cATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTAtttt                                                         <  1:1997279/91‑1 (MQ=255)
                                 tatTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTCCTAAGTGTTATTATTTTaa                                                       >  1:1488654/1‑91 (MQ=255)
                                    tttAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGc                                                    >  1:620143/1‑91 (MQ=255)
                                      tAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCAt                                                  <  1:1161759/91‑1 (MQ=255)
                                       aGCGATTTAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATc                                                 <  1:2074013/91‑1 (MQ=255)
                                       aGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATc                                                 >  1:2448908/1‑91 (MQ=255)
                                       aGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATc                                                 >  1:54287/1‑91 (MQ=255)
                                       aGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATc                                                 >  1:874254/1‑91 (MQ=255)
                                        gCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATct                                                <  1:1039725/91‑1 (MQ=255)
                                         cGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctt                                               >  1:550023/1‑91 (MQ=255)
                                         cGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctt                                               >  1:1999439/1‑91 (MQ=255)
                                          gATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATcttt                                              <  1:858992/91‑1 (MQ=255)
                                             tAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagc                                           <  1:1633119/91‑1 (MQ=255)
                                                 tgtgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagcttta                                       <  1:3818799/91‑1 (MQ=255)
                                                 tgtgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagcttta                                       >  1:3594633/1‑91 (MQ=255)
                                                 tgtgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagcttta                                       >  1:3101757/1‑91 (MQ=255)
                                                 tgtgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagcttta                                       <  1:2975461/91‑1 (MQ=255)
                                                  gtgGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATctttagctttag                                      <  1:2219140/91‑1 (MQ=255)
                                                   tgGCATTTTGCTATTTAAAAAATAAAAAAGAGACGACAATGCCATCACTTTTTTATAAGCGTTATTATTTTAAAGCATCTTTAGCTTTAGt                                     >  1:3314658/1‑91 (MQ=255)
                                                     gCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTa                                   >  1:2212605/1‑91 (MQ=255)
                                                        ttttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTa                                <  1:150377/91‑1 (MQ=255)
                                                          ttGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAAt                              >  1:4015848/1‑91 (MQ=255)
                                                           tGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAAtt                             >  1:3925865/1‑91 (MQ=255)
                                                            gCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTg                            >  1:1155515/1‑91 (MQ=255)
                                                            gCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTg                            >  1:1136966/1‑91 (MQ=255)
                                                            gCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTg                            >  1:1536410/1‑91 (MQ=255)
                                                             cTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGa                           >  1:2669458/1‑91 (MQ=255)
                                                             cTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGa                           >  1:2604610/1‑91 (MQ=255)
                                                              tATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAg                          >  1:1196752/1‑91 (MQ=255)
                                                                tttAAAAAATAAAAAAGAGATGACAATGCAATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTCAGCa                        >  1:3447386/1‑91 (MQ=255)
                                                                  taaaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCaaa                      <  1:3941544/91‑1 (MQ=255)
                                                                   aaaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAAt                     >  1:1698550/1‑91 (MQ=255)
                                                                    aaaaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATg                    >  1:877965/1‑91 (MQ=255)
                                                                      aaataaaaaaGAGATGACAATGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCt                  <  1:2867064/91‑1 (MQ=255)
                                                                      aaataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCt                  >  1:2044259/1‑91 (MQ=255)
                                                                       tataaaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtt                 <  1:83107/90‑1 (MQ=255)
                                                                         taaaaaaGAGATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt               >  1:2568938/1‑91 (MQ=37)
                                                                         taaaaaaGAGATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt               >  1:2568111/1‑91 (MQ=37)
                                                                         taaaaaaGAGATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt               >  1:2568493/1‑91 (MQ=37)
                                                                         taaaaaaGAGATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCtttt               >  1:2568180/1‑91 (MQ=37)
                                                                          aaaaaaGAGATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCttttt              >  1:3499046/1‑91 (MQ=37)
                                                                           aaaaaGAGATGACAATGCCATCTCTTTTTTATACGTGATATTATTTAAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTc             >  1:3303951/1‑91 (MQ=255)
                                                                           aaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTc             >  1:682025/1‑91 (MQ=255)
                                                                           aaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTc             >  1:2776357/1‑91 (MQ=255)
                                                                           aaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTc             >  1:1229095/1‑91 (MQ=255)
                                                                           aaaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGATATTATTTCAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTc             >  1:3303468/1‑91 (MQ=255)
                                                                            aaaaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCg            <  1:595445/91‑1 (MQ=255)
                                                                              aaGAGATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTc          >  1:924392/1‑91 (MQ=38)
                                                                              aaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTc          <  1:366193/91‑1 (MQ=255)
                                                                              aaGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTc          <  1:3208980/91‑1 (MQ=255)
                                                                               agagATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCa         >  1:488843/1‑91 (MQ=38)
                                                                               agagATGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCa         >  1:2041074/1‑91 (MQ=38)
                                                                                gagaTGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAg        >  1:2732372/1‑91 (MQ=38)
                                                                                gagaTGGCATTGTCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAg        >  1:2732398/1‑91 (MQ=38)
                                                                                    tGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAt    <  1:3163968/91‑1 (MQ=255)
                                                                                    tGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAt    <  1:1570348/91‑1 (MQ=255)
                                                                                    tGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAt    >  1:1848927/1‑91 (MQ=255)
                                                                                    tGACAATGCCATCTCTTGTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAt    <  1:339534/91‑1 (MQ=255)
                                                                                     gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt   <  1:3843129/91‑1 (MQ=255)
                                                                                     gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt   <  1:938091/91‑1 (MQ=255)
                                                                                     gACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAAtt   >  1:2433643/1‑91 (MQ=255)
                                                                                      aCAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTg  >  1:3613251/1‑91 (MQ=255)
                                                                                      aCAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTg  >  1:293195/1‑91 (MQ=255)
                                                                                      |                                                                                          
TGTTCTCTTAAATGCATTTTAGCAAATATGTCATATTTTAGCGATTAAATGTGGCATTTTGCTATTTAAAAAATAAAAAAGAGATGACAATGCCATCTCTTTTTTATAAGTGTTATTATTTTAAAGCATCTTTAGCTTTAGTTACTAATTGAGCAAATGCTTTTTCGTCAGAAATTG  >  USA300TCH1516_ALE/1779721‑1779897

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: