Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 583,951 A→T 8.0% intergenic (+35/‑237) 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,9510AT 8.0% 671.0 / 35.3 201intergenic (+35/‑237)rplA/rplJ50S ribosomal protein L1/50S ribosomal protein L10
Reads supporting (aligned to +/- strand):  ref base A (99/86);  new base T (10/6);  total (109/92)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 215 total aligned reads displayed.

CTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTAT  >  USA300TCH1516_ALE/583861‑584041
                                                                                          |                                                                                          
cTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAata                                                                                            >  2:670932/1‑91 (MQ=255)
cTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAata                                                                                            <  1:3052936/91‑1 (MQ=255)
  gTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATc                                                                                          <  2:644221/91‑1 (MQ=255)
        acaacaATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGACTGAAAGATATTTCTTTCAg                                                                                    <  1:922909/91‑1 (MQ=38)
         caacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGt                                                                                   <  1:2300264/91‑1 (MQ=255)
          aacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTc                                                                                  <  1:3721789/91‑1 (MQ=255)
           acaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCt                                                                                 >  1:2901796/1‑91 (MQ=255)
           acaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCt                                                                                 <  2:3047943/91‑1 (MQ=255)
           acaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGACTGAAAGATATTTCTTTCAgcct                                                                                 <  1:1424086/91‑4 (MQ=37)
            caaTGGGTCCTGGTGTTAAAATTGATACTGCAAGTTTCAAATAGTAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTg                                                                                <  1:3881661/91‑1 (MQ=255)
            caaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTg                                                                                >  2:2690637/1‑91 (MQ=255)
               tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaa                                                                             >  1:1240271/1‑91 (MQ=255)
               tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaa                                                                             <  1:2951138/91‑1 (MQ=255)
               tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGACTGAAAGATATTTCTTTCAGCCTGTaa                                                                             <  2:1997376/91‑1 (MQ=37)
                gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                            <  1:1639840/91‑1 (MQ=255)
                gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                            <  2:3792785/91‑1 (MQ=255)
                  gTCCTGGAGTTAAAATTGATACTGCAGGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                                          >  2:2220671/1‑91 (MQ=255)
                  gTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                                          >  1:3393998/1‑91 (MQ=255)
                  gTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                                          >  2:2563597/1‑91 (MQ=255)
                  gTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                                          <  1:2636227/91‑1 (MQ=255)
                  gTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGACTGAAAGATATTTCTTTCAGCCTGTaattg                                                                          <  2:1077700/91‑4 (MQ=25)
                   tCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAat                                                                         >  2:2813448/1‑91 (MQ=255)
                       ggAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatatat                                                                     <  1:2309861/91‑1 (MQ=255)
                        gAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATAtt                                                                    <  2:423103/91‑1 (MQ=255)
                         aGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTg                                                                   >  2:193646/1‑91 (MQ=255)
                            tAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACa                                                                >  2:3540504/1‑91 (MQ=255)
                              aaaTTGATACTGAAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAAt                                                              >  1:2996190/1‑91 (MQ=255)
                                aTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATaa                                                            >  1:2203558/1‑91 (MQ=255)
                                aTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATaa                                                            >  1:2203355/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)
                                  tGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGt                                                          <  1:1126386/91‑1 (MQ=255)
                                    aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                        >  1:2223847/1‑91 (MQ=255)
                                    aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                        >  2:1611111/1‑91 (MQ=255)
                                     tACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAAt                                                       <  2:3790784/91‑1 (MQ=255)
                                       cTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAAttt                                                     <  1:2453203/91‑1 (MQ=255)
                                             gTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGt                                               >  1:3326901/1‑91 (MQ=255)
                                              tttCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGtt                                              <  2:2305077/91‑1 (MQ=255)
                                              tttCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGtt                                              <  2:2234309/91‑1 (MQ=255)
                                                tCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTtat                                            <  2:359754/91‑1 (MQ=255)
                                                tCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTtat                                            <  1:2385657/91‑1 (MQ=255)
                                                  aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                          >  2:2627453/1‑91 (MQ=255)
                                                   aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATAtt                                         >  1:3771945/1‑91 (MQ=255)
                                                     taataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTAc                                       <  2:2045959/91‑1 (MQ=255)
                                                     taataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTAc                                       <  2:2191541/91‑1 (MQ=255)
                                                      aataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTTACAATAAGTAATTTCCAAGTTATATTACt                                      >  2:2290768/1‑91 (MQ=255)
                                                      aataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACt                                      <  1:350942/91‑1 (MQ=255)
                                                       ataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACtt                                     <  1:1293905/91‑1 (MQ=255)
                                                       ataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACtt                                     >  1:1400015/1‑91 (MQ=255)
                                                        taaATGATATAAACAATTACAGGCTAAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTa                                    <  2:1910597/91‑1 (MQ=255)
                                                          aaTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTAtt                                  >  2:350129/1‑91 (MQ=255)
                                                           aTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtg                                 <  1:1693743/91‑1 (MQ=255)
                                                           aTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtg                                 >  1:3148080/1‑91 (MQ=255)
                                                            tGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgt                                >  2:3683342/1‑91 (MQ=255)
                                                            tGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgt                                >  2:3683287/1‑91 (MQ=255)
                                                             gATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgtt                               >  1:3408998/1‑90 (MQ=255)
                                                             gATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgtg                               >  1:3542663/1‑91 (MQ=255)
                                                             gATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgtg                               >  2:3333946/1‑91 (MQ=255)
                                                              atatAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGa                              >  2:1222401/1‑91 (MQ=255)
                                                              atatAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGa                              >  1:414022/1‑91 (MQ=255)
                                                              atatAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGa                              >  1:886116/1‑91 (MQ=255)
                                                              atatAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGat                             >  1:886375/1‑91 (MQ=255)
                                                                 taAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGatta                           >  1:607645/1‑91 (MQ=255)
                                                                  aaaCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGCCAATAAGTAATTTCCAAGTTATATTACTTATTGTGattat                          >  1:2491922/1‑91 (MQ=255)
                                                                    aCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGAGAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTAttt                        <  1:3355214/91‑1 (MQ=255)
                                                                      aaTTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTa                      <  2:3250693/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)
                                                                          aCAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTa                  >  1:59683/1‑91 (MQ=255)
                                                                          aCAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTa                  >  2:302026/1‑91 (MQ=255)
                                                                           cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa                 <  2:1853869/91‑1 (MQ=255)
                                                                           cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa                 <  2:785071/91‑1 (MQ=255)
                                                                            aGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAg                <  1:1059842/91‑1 (MQ=255)
                                                                            aGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAg                <  1:1665099/91‑1 (MQ=255)
                                                                             ggCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa               >  2:3030317/1‑91 (MQ=255)
                                                                             ggCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa               >  2:1552034/1‑91 (MQ=255)
                                                                             ggCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa               <  2:2278633/91‑1 (MQ=255)
                                                                             gacTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATTTTACTTATTGAGATTATTTTACCTAAGa               >  1:755556/3‑91 (MQ=37)
                                                                             gacTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa               >  2:1087962/3‑91 (MQ=38)
                                                                             gacTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa               >  2:706821/3‑91 (MQ=38)
                                                                             gacTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa               >  1:3465073/3‑91 (MQ=38)
                                                                                tGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAg            <  1:2113509/91‑1 (MQ=255)
                                                                                 gaaagaaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGt           <  2:576036/91‑1 (MQ=255)
                                                                                  aaagaaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTa          <  1:264533/91‑1 (MQ=255)
                                                                                    agaaaTATCTTTCAGACTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTTCTTATTGTGATTATTTTACCTAAGACAGTAgg        >  1:3302725/1‑91 (MQ=255)
                                                                                     gaaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGa       >  1:899067/1‑91 (MQ=255)
                                                                                      aaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAg      <  2:2429519/91‑1 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATCTATTGACAATAAGTAATTTTCAAGATATATTACATATTGAGATTATTTTACCTAAGACAGTAGGAGt     >  1:1445238/1‑91 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt     >  1:1446182/1‑91 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt     >  1:348070/1‑91 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt     >  1:2714126/1‑91 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt     >  1:73315/1‑91 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt     >  2:177588/1‑91 (MQ=255)
                                                                                        atatCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGtt    >  2:1434952/1‑91 (MQ=255)
                                                                                         tatCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGtta   >  1:2270237/1‑91 (MQ=255)
                                                                                         tatCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGtta   >  2:1306137/1‑91 (MQ=255)
                                                                                          atCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGttat  >  1:458134/1‑91 (MQ=255)
                                                                                          atCTTCCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGttat  <  2:3675059/91‑1 (MQ=255)
                                                                                          |                                                                                          
CTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTAT  >  USA300TCH1516_ALE/583861‑584041

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