Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 583,939 G→A 6.0% intergenic (+23/‑249) rplA → / → rplJ 50S ribosomal protein L1/50S ribosomal protein L10

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE583,9390GA 6.0% 663.1 / 15.6 204intergenic (+23/‑249)rplA/rplJ50S ribosomal protein L1/50S ribosomal protein L10
Reads supporting (aligned to +/- strand):  ref base G (99/91);  new base A (5/8);  total (104/100)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

TCAAATCTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGA  >  USA300TCH1516_ALE/583849‑584028
                                                                                          |                                                                                         
tCAAATCTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAgg                                                                                           <  1:3682795/91‑1 (MQ=255)
tCAAATCTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAgg                                                                                           >  1:1465767/1‑91 (MQ=255)
tCAAATCTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAgg                                                                                           <  1:2596351/91‑1 (MQ=255)
tCAAATCTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAgg                                                                                           <  2:2092904/91‑1 (MQ=255)
   aaTCTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTg                                                                                        >  1:2207213/1‑91 (MQ=255)
   aaTCTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTg                                                                                        >  1:2720981/1‑91 (MQ=255)
    aTCTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTga                                                                                       >  1:606045/1‑91 (MQ=255)
      cTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTgaaa                                                                                     >  1:305472/1‑91 (MQ=255)
      cTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTgaaa                                                                                     >  1:1853869/1‑91 (MQ=255)
       tGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGCTACTGCAAGTTTCAAATAATACATGATATAAACAATTACAGGCTgaaag                                                                                    >  2:158982/1‑91 (MQ=255)
       tGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTgaaag                                                                                    >  2:2243899/1‑91 (MQ=255)
       tGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTgaaag                                                                                    <  2:1675212/91‑1 (MQ=255)
       tGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTgaaag                                                                                    <  1:2050049/91‑1 (MQ=255)
        gTTGCTGTAACTGCAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTgaaaga                                                                                   <  1:1703626/91‑1 (MQ=255)
         ttGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTgaaagaa                                                                                  >  1:555942/1‑91 (MQ=255)
            cTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTgaaagatatt                                                                               <  2:2185906/91‑5 (MQ=255)
              gTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATc                                                                             <  2:644243/91‑1 (MQ=255)
              gTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATc                                                                             <  2:643779/91‑1 (MQ=255)
              gTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATc                                                                             <  2:1544536/91‑1 (MQ=255)
                aaCTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCtt                                                                           >  2:1573446/1‑91 (MQ=255)
                  cTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTc                                                                         >  2:2398412/1‑91 (MQ=255)
                    acaacaATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGACTGAAAGATATTTCTTTCAg                                                                       <  1:922909/91‑1 (MQ=38)
                     caacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGt                                                                      <  1:2300264/91‑1 (MQ=255)
                      aacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTc                                                                     >  1:1422309/1‑91 (MQ=255)
                       acaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCt                                                                    >  1:2901796/1‑91 (MQ=255)
                        caaTGGGTCCTGGTGTTAAAATTGATACTGCAAGTTTCAAATAGTAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTg                                                                   <  1:3881661/91‑1 (MQ=255)
                        caaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTg                                                                   >  2:2690637/1‑91 (MQ=255)
                           tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTgaa                                                                                      >  2:158449/1‑69 (MQ=255)
                           tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTgaa                                                                                      <  1:158449/69‑1 (MQ=255)
                           tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaa                                                                >  1:1240271/1‑91 (MQ=255)
                           tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaa                                                                <  1:2951138/91‑1 (MQ=255)
                           tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGACTGAAAGATATTTCTTTCAGCCTGTaa                                                                <  2:1997376/91‑1 (MQ=37)
                            gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                               <  2:3792785/91‑1 (MQ=255)
                            gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                               <  1:1639966/91‑1 (MQ=255)
                              gTCCTGGAGTTAAAATTGATACTGCACGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                             >  2:2220992/1‑91 (MQ=255)
                              gTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                             >  2:6904/1‑91 (MQ=255)
                              gTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                             <  1:2636296/91‑1 (MQ=255)
                              gTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                             >  2:2563597/1‑91 (MQ=255)
                               tCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAat                                                            >  2:1057695/1‑91 (MQ=255)
                               tCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAat                                                            >  2:2813448/1‑91 (MQ=255)
                               tCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAat                                                            >  1:823838/1‑91 (MQ=255)
                                    gAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATAtt                                                       <  2:423103/91‑1 (MQ=255)
                                     aGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTg                                                      >  2:193646/1‑91 (MQ=255)
                                       ttAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTACAAATATATTGAc                                                    >  2:3226025/1‑91 (MQ=255)
                                        tAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACa                                                   >  2:3541472/1‑91 (MQ=255)
                                          aaaTTGATACTGAAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAAt                                                 >  1:2996190/1‑91 (MQ=255)
                                            aTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATaa                                               >  1:2203557/1‑91 (MQ=255)
                                            aTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATaa                                               >  2:2431095/1‑91 (MQ=255)
                                            aTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATaa                                               >  2:2329036/1‑91 (MQ=255)
                                             ttGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAg                                              >  1:1492464/1‑91 (MQ=255)
                                              tGATACTGCAAGTTTCGAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGt                                             <  1:2132228/91‑1 (MQ=255)
                                              tGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGt                                             <  1:1126386/91‑1 (MQ=255)
                                                aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                           >  1:2223847/1‑91 (MQ=255)
                                                 tACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAAt                                          <  2:3790784/91‑1 (MQ=255)
                                                  aCTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATCGACAATAAGTAAtt                                         <  2:3294404/91‑1 (MQ=255)
                                                   cTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAAttt                                        >  1:3540810/1‑91 (MQ=255)
                                                   cTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAAttt                                        <  1:2453203/91‑1 (MQ=255)
                                                         gTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGt                                  >  1:3326901/1‑91 (MQ=255)
                                                          tttCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGtt                                 <  2:2305077/91‑1 (MQ=255)
                                                            tCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTtat                               <  2:1217195/91‑1 (MQ=255)
                                                            tCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTtat                               <  1:2385657/91‑1 (MQ=255)
                                                            tCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTtat                               <  2:359754/91‑1 (MQ=255)
                                                              aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                             >  2:2627453/1‑91 (MQ=255)
                                                               aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATAtt                            <  2:3038147/91‑1 (MQ=255)
                                                               aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATAtt                            >  1:3596073/1‑91 (MQ=255)
                                                               aataatAAATGATATAAACAATTACACGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATAtt                            >  2:2209850/1‑91 (MQ=255)
                                                                 taataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTAc                          <  2:2045959/91‑1 (MQ=255)
                                                                  aataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTTACAATAAGTAATTTCCAAGTTATATTACt                         >  2:2290768/1‑91 (MQ=255)
                                                                  aataaATGATACAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATAGTGACAATAAGTAATTTCCAAGTTATATTACt                         <  2:2969284/91‑1 (MQ=255)
                                                                   ataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACtt                        <  1:1293905/91‑1 (MQ=255)
                                                                    taaATGATATAAACAATTACAGGCTAAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTa                       <  2:1910597/91‑1 (MQ=255)
                                                                     aaaTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTAt                      >  1:1389277/1‑91 (MQ=255)
                                                                      aaTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTAtt                     >  2:350129/1‑91 (MQ=255)
                                                                       aTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtg                    >  1:3148080/1‑91 (MQ=255)
                                                                        tGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgt                   >  2:3683287/1‑91 (MQ=255)
                                                                        tGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgt                   >  2:3683342/1‑91 (MQ=255)
                                                                        tGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgt                   <  2:769397/91‑1 (MQ=255)
                                                                         gATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgtt                  >  1:3408998/1‑90 (MQ=255)
                                                                          atatAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGa                 >  1:886571/1‑91 (MQ=255)
                                                                          atatAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGa                 >  1:886338/1‑91 (MQ=255)
                                                                          atatAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGa                 >  1:414022/1‑91 (MQ=255)
                                                                             taAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGatta              >  1:607645/1‑91 (MQ=255)
                                                                              aaaCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGCCAATAAGTAATTTCCAAGTTATATTACTTATTGTGattat             >  1:2491922/1‑91 (MQ=255)
                                                                               aaCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGattatt            >  1:3766356/1‑91 (MQ=255)
                                                                                aCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGAGAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTAttt           <  1:3355214/91‑1 (MQ=255)
                                                                                 cAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTAtttt          <  1:3506817/91‑1 (MQ=255)
                                                                                   tttACAGACTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTAc        >  2:899914/2‑91 (MQ=25)
                                                                                   tttACAGACTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTAc        >  1:968793/2‑91 (MQ=25)
                                                                                   aTTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTAc        >  1:1086021/1‑91 (MQ=255)
                                                                                    ttACAGACTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTAcc       >  2:3393524/1‑91 (MQ=37)
                                                                                     tACAGACTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCt      >  1:3383954/1‑91 (MQ=37)
                                                                                      aCAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTa     >  2:302605/1‑91 (MQ=255)
                                                                                      aCAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTa     >  1:59563/1‑91 (MQ=255)
                                                                                       cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa    <  2:1853869/91‑1 (MQ=255)
                                                                                       cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa    >  2:1084695/1‑91 (MQ=255)
                                                                                       cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa    >  1:648487/1‑91 (MQ=255)
                                                                                       cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa    <  2:785071/91‑1 (MQ=255)
                                                                                        aGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAg   <  1:1665206/91‑1 (MQ=255)
                                                                                        aGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAg   <  1:1059842/91‑1 (MQ=255)
                                                                                         ggCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa  <  2:2278633/91‑1 (MQ=255)
                                                                                          |                                                                                         
TCAAATCTGTTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGA  >  USA300TCH1516_ALE/583849‑584028

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: