Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,239,562 A→G 7.9% intergenic (+29/‑161) USA300HOU_RS06160 → / → recG hypothetical protein/ATP‑dependent DNA helicase RecG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,239,5620AG 7.9% 725.2 / 34.1 235intergenic (+29/‑161)USA300HOU_RS06160/recGhypothetical protein/ATP‑dependent DNA helicase RecG
Reads supporting (aligned to +/- strand):  ref base A (117/99);  new base G (9/10);  total (126/109)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

TAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTC  >  USA300TCH1516_ALE/1239472‑1239651
                                                                                           |                                                                                         
  gAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATaaa                                                                                         >  1:1295800/1‑91 (MQ=255)
    aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                       >  2:3042142/1‑91 (MQ=255)
    aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                       >  2:14538/1‑91 (MQ=255)
    aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                       >  2:1616792/1‑91 (MQ=255)
    aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTATTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATGAatc                                                                                       <  1:3070944/91‑1 (MQ=255)
    aGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATATCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatc                                                                                       >  2:3184506/1‑91 (MQ=255)
     gTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatca                                                                                      <  2:1837036/91‑1 (MQ=255)
      tGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAatcat                                                                                     >  1:716274/1‑91 (MQ=255)
        gAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAAtcatca                                                                                   >  1:2728488/1‑91 (MQ=255)
         aaGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAg                                                                                  >  2:1906368/1‑91 (MQ=255)
           gTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGtt                                                                                <  2:1702709/91‑1 (MQ=255)
             tCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTgg                                                                              <  1:1450705/91‑1 (MQ=255)
             tCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTgg                                                                              <  1:2539286/91‑1 (MQ=255)
              cATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGt                                                                             >  2:3041361/1‑91 (MQ=255)
               aTGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTa                                                                            >  2:536903/1‑91 (MQ=255)
                 gAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGt                                                                          >  1:3469647/1‑91 (MQ=255)
                  aaGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtt                                                                         >  1:3225600/1‑91 (MQ=255)
                  aaGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtt                                                                         >  1:996159/1‑91 (MQ=255)
                   aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                        <  1:3694659/91‑1 (MQ=255)
                   aGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGttt                                                                        >  2:2019190/1‑91 (MQ=255)
                    ggTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtttt                                                                       <  2:1880289/91‑1 (MQ=255)
                    ggTGGACAACCAATTTATCAATATTTCTTTTCAGCAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtttt                                                                       >  2:3396286/1‑91 (MQ=255)
                     ttGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGttttt                                                                      <  2:3853740/90‑1 (MQ=18)
                      tGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGtttttt                                                                     >  2:2777059/1‑91 (MQ=255)
                        gACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAt                                                                   >  1:2146653/1‑91 (MQ=255)
                         aCAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAtt                                                                  <  1:1072607/91‑1 (MQ=17)
                         aCAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAACTGATGATTTATCATTGGTAGTTTTTTAtt                                                                  <  1:2681907/91‑1 (MQ=18)
                           aaCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAtttt                                                                <  1:1594614/91‑1 (MQ=255)
                           aaCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAtttt                                                                <  1:2932250/91‑1 (MQ=255)
                           aaCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTAtttt                                                                >  1:442885/1‑91 (MQ=255)
                            aCCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTg                                                               <  1:962872/91‑1 (MQ=255)
                             ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                              <  2:3342831/91‑1 (MQ=255)
                             ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                              <  2:2001241/91‑1 (MQ=255)
                             ccAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGc                                                              <  1:1959414/91‑1 (MQ=255)
                              cAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCt                                                             <  2:1791674/91‑1 (MQ=255)
                               aaTTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTa                                                            >  2:1735672/1‑91 (MQ=255)
                                 tttATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTAtt                                                          <  2:3276818/91‑1 (MQ=255)
                                      cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                     <  2:826824/91‑1 (MQ=255)
                                      cAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGt                                                     >  2:2636004/1‑91 (MQ=255)
                                       aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                                    <  1:573893/91‑1 (MQ=255)
                                       aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                                    <  2:436151/91‑1 (MQ=255)
                                       aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                                    >  1:1719349/1‑91 (MQ=255)
                                       aaTATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTg                                                    >  1:2775217/1‑91 (MQ=255)
                                        atatTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGa                                                   <  1:163005/91‑1 (MQ=255)
                                         tatTTCTTTTCAGTAGAATAAAGATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTATTTTATTTTGCTATTTTAGTGat                                                  <  2:3512055/91‑1 (MQ=255)
                                          atTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTgtta                                                 <  2:2190959/91‑4 (MQ=255)
                                           tttCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGatat                                                <  2:274377/91‑1 (MQ=255)
                                           tttCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGatat                                                <  2:273440/91‑1 (MQ=255)
                                             tCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTg                                              >  2:3492851/1‑91 (MQ=255)
                                             tCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTg                                              >  1:185153/1‑91 (MQ=255)
                                              cTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGc                                             >  1:3783622/1‑91 (MQ=255)
                                                tttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCgg                                           >  2:379735/1‑91 (MQ=255)
                                                 ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                          <  2:1392618/91‑1 (MQ=255)
                                                 ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                          <  1:2777059/91‑1 (MQ=255)
                                                 ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                          >  1:3696674/1‑91 (MQ=255)
                                                 ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                          <  2:2582502/91‑1 (MQ=255)
                                                 ttCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCggg                                          <  2:2585536/91‑1 (MQ=255)
                                                  tCAGTAGAATAGAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGt                                         <  1:3513735/91‑1 (MQ=255)
                                                    aGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTa                                       >  2:2488359/1‑91 (MQ=255)
                                                    aGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTa                                       <  1:821821/91‑1 (MQ=255)
                                                      tAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaaa                                     <  1:3855038/91‑1 (MQ=255)
                                                       aGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaaaa                                    <  2:2314266/91‑1 (MQ=255)
                                                       aGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTaaaa                                    <  2:2047336/91‑1 (MQ=255)
                                                        gAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAg                                   >  1:2265623/1‑91 (MQ=255)
                                                        gAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAg                                   <  1:3376188/91‑1 (MQ=255)
                                                             aaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCg                              >  2:364143/1‑91 (MQ=255)
                                                             aaaaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCg                              >  1:341898/1‑91 (MQ=255)
                                                               aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                            >  2:926744/1‑91 (MQ=255)
                                                               aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                            <  1:2398602/91‑1 (MQ=255)
                                                               aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                            <  1:3456434/91‑1 (MQ=255)
                                                               aaTTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGtt                            <  1:980065/91‑1 (MQ=255)
                                                                 tttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGctct                          >  1:3123436/1‑91 (MQ=255)
                                                                 tttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTtct                          >  1:3188914/1‑91 (MQ=255)
                                                                 tttAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTtct                          >  2:3107932/1‑91 (MQ=255)
                                                                   taaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCg                        >  1:1238317/1‑91 (MQ=255)
                                                                   taaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCg                        >  1:554669/1‑91 (MQ=255)
                                                                    aaaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTGTTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGa                       <  1:3271713/91‑1 (MQ=255)
                                                                     aaataaaaAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAg                      >  2:2364634/1‑91 (MQ=255)
                                                                          aaaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCt                 >  1:457647/1‑91 (MQ=255)
                                                                           aaaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTa                >  1:2034142/1‑91 (MQ=255)
                                                                            aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa               >  2:2303288/1‑91 (MQ=18)
                                                                            aaaCTACCAACTGATGATTTATCATTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa               >  1:3143115/1‑91 (MQ=18)
                                                                            aaaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTaa               <  2:86586/91‑1 (MQ=255)
                                                                             aaCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAAc              >  1:1713417/1‑91 (MQ=255)
                                                                              aCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACa             >  2:115424/1‑91 (MQ=255)
                                                                              aCTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACa             >  2:507824/1‑91 (MQ=255)
                                                                               cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa            >  1:710420/1‑91 (MQ=255)
                                                                               cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa            >  2:2847506/1‑91 (MQ=255)
                                                                               cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa            >  1:3253683/1‑91 (MQ=255)
                                                                               cTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACaa            >  2:164844/1‑91 (MQ=255)
                                                                                tgccAATGATAAATCATCAGTTGGGAGTTTTTTATTTTGCTATTTTAGTGTTATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAAt           <  2:3599512/89‑1 (MQ=255)
                                                                                 aCCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATg          >  1:64980/1‑91 (MQ=255)
                                                                                  ccAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGt         <  1:3151271/91‑1 (MQ=255)
                                                                                  ccAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGt         >  2:1201519/1‑91 (MQ=255)
                                                                                     aTGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCAt      >  1:2500144/1‑91 (MQ=255)
                                                                                       tGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATg     <  1:1880934/91‑1 (MQ=255)
                                                                                       tGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATg     <  1:423071/91‑1 (MQ=255)
                                                                                       tGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATg     <  2:2614359/91‑1 (MQ=255)
                                                                                         aTAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGtt   <  2:3011240/91‑1 (MQ=255)
                                                                                          tAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTc  <  2:130295/91‑1 (MQ=255)
                                                                                           |                                                                                         
TAGAAGTGGAAGTTCATGAAGGTGGACAACCAATTTATCAATATTTCTTTTCAGTAGAATAAAAATTTAAAATAAAAAACTACCAATGATAAATCATCAGTTGGTAGTTTTTTATTTTGCTATTTTAGTGATATTGCGGGTTAAAAGTATCGTTCTCGAGTTGCTAACAATGTCATGTTC  >  USA300TCH1516_ALE/1239472‑1239651

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: