Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 583,960 T→C 6.7% intergenic (+44/‑228) 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,9600TC 6.7% 654.9 / 22.1 210intergenic (+44/‑228)rplA/rplJ50S ribosomal protein L1/50S ribosomal protein L10
Reads supporting (aligned to +/- strand):  ref base T (107/88);  new base C (11/3);  total (119/91)
Fisher's exact test for biased strand distribution p-value = 9.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 212 total aligned reads displayed.

CAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACT  >  USA300TCH1516_ALE/583870‑584049
                                                                                          |                                                                                         
caacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGt                                                                                           <  1:2300264/91‑1 (MQ=255)
 aacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTc                                                                                          <  1:1123070/91‑1 (MQ=255)
 aacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTc                                                                                          >  1:1422309/1‑91 (MQ=255)
  acaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCt                                                                                         >  1:2917182/1‑91 (MQ=255)
   caaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTg                                                                                        >  2:2690637/1‑91 (MQ=255)
      tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaa                                                                                     <  1:2951138/91‑1 (MQ=255)
       gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                                    <  1:3401819/91‑1 (MQ=255)
       gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                                    <  2:3792785/91‑1 (MQ=255)
       gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                                    <  1:1639966/91‑1 (MQ=255)
       gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                                    >  1:3072026/1‑91 (MQ=255)
       gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGACTGAAAGATATTTCTTTCAGCCTGTaat                                                                                    <  1:606382/91‑2 (MQ=25)
         gTCCTGGAGTTAAAATTGATACTGCAGGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                                                  >  2:2220671/1‑91 (MQ=255)
         gTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGACTGAAAGATATTTCTTTCAGCCTGTaattg                                                                                  <  2:1077700/91‑4 (MQ=25)
          tCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAat                                                                                 >  2:2813448/1‑91 (MQ=255)
            cTGGAGTTAAAATTGATACTGCAAGTTTCAAATGATAAATGGTATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatat                                                                               <  2:3881583/91‑1 (MQ=255)
              ggAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatatat                                                                             >  1:3250693/1‑91 (MQ=255)
              ggAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatatat                                                                             <  1:2309861/91‑1 (MQ=255)
                aGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTg                                                                           >  2:193137/1‑91 (MQ=255)
                  ttAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTACAAATATATTGAc                                                                         >  2:3226025/1‑91 (MQ=255)
                   tAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACa                                                                        >  2:3540504/1‑91 (MQ=255)
                   tAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACa                                                                        >  2:3541472/1‑91 (MQ=255)
                    aaaaTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACaa                                                                       >  1:989022/1‑91 (MQ=255)
                     aaaTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAAt                                                                      >  1:2424609/1‑91 (MQ=255)
                     aaaTTGATACTGAAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAAt                                                                      >  1:2996190/1‑91 (MQ=255)
                       aTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATaa                                                                    >  2:2329036/1‑91 (MQ=255)
                       aTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATaa                                                                    >  1:2203557/1‑91 (MQ=255)
                        ttGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAg                                                                   >  1:1492464/1‑91 (MQ=255)
                        ttGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAg                                                                   <  2:2589567/91‑1 (MQ=255)
                         tGATACTGCAAGTTTCGAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGt                                                                  <  1:2132228/91‑1 (MQ=255)
                           aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                                >  2:1611111/1‑91 (MQ=255)
                           aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                                >  1:395704/1‑91 (MQ=255)
                           aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                                >  1:2223847/1‑91 (MQ=255)
                                     tttCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGtt                                                      <  2:2305077/91‑1 (MQ=255)
                                     tttCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGtt                                                      <  2:2234309/91‑1 (MQ=255)
                                     tttCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGtt                                                      <  2:2108016/91‑1 (MQ=255)
                                       tCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTtat                                                    <  2:359754/91‑1 (MQ=255)
                                        cAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTtata                                                   <  2:1874055/91‑1 (MQ=255)
                                          aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTg                                                                           <  2:2541171/65‑1 (MQ=255)
                                          aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTg                                                                           >  1:2541171/1‑65 (MQ=255)
                                          aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATAtt                                                 >  1:1731519/1‑91 (MQ=255)
                                          aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATAtt                                                 <  2:3038147/91‑1 (MQ=255)
                                          aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATAtt                                                 >  1:1235419/1‑91 (MQ=255)
                                          aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATAtt                                                 <  1:2467004/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)
                                             aataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACt                                              <  2:1939445/91‑1 (MQ=255)
                                             aataaATGATACAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATAGTGACAATAAGTAATTTCCAAGTTATATTACt                                              <  2:2969284/91‑1 (MQ=255)
                                                aaaTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTAt                                           >  1:1389277/1‑91 (MQ=255)
                                                 aaTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTAtt                                          >  2:350129/1‑91 (MQ=255)
                                                  aTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtg                                         >  1:2353360/1‑91 (MQ=255)
                                                   tGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgt                                        >  2:3683287/1‑91 (MQ=255)
                                                   tGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgt                                        >  2:3683342/1‑91 (MQ=255)
                                                    gATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgtg                                       >  2:3333946/1‑91 (MQ=255)
                                                    gATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgtg                                       >  1:3542663/1‑91 (MQ=255)
                                                     atatAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGa                                      >  1:886338/1‑91 (MQ=255)
                                                     atatAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGa                                      >  2:1222401/1‑91 (MQ=255)
                                                     atatAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGa                                      >  1:886571/1‑91 (MQ=255)
                                                        taAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGatta                                   >  1:608015/1‑91 (MQ=255)
                                                         aaaCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGattat                                  <  1:3868796/91‑1 (MQ=255)
                                                          aaCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGattatt                                 >  1:3766356/1‑91 (MQ=255)
                                                           aCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTAttt                                <  1:946127/91‑1 (MQ=255)
                                                           aCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTAttt                                >  2:548570/1‑91 (MQ=255)
                                                              tttACAGACTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTAc                             >  2:899914/2‑91 (MQ=25)
                                                                tACAGACTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCt                           >  1:3383954/1‑91 (MQ=37)
                                                                 aCAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTa                          >  2:302605/1‑91 (MQ=255)
                                                                 aCAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTa                          >  1:59563/1‑91 (MQ=255)
                                                                  cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa                         <  1:416530/91‑1 (MQ=255)
                                                                  cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa                         <  2:785071/91‑1 (MQ=255)
                                                                   aGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAg                        <  1:1070614/91‑1 (MQ=255)
                                                                   aGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAg                        <  1:1665206/91‑1 (MQ=255)
                                                                   aGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAg                        >  1:1284270/1‑91 (MQ=255)
                                                                    ggCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa                       <  2:2278633/91‑1 (MQ=255)
                                                                    ggCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa                       >  2:1552034/1‑91 (MQ=255)
                                                                    ggCTGAAAGAAATATCTATCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa                       <  1:3330108/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)
                                                                         aaagaaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTa                  <  1:264533/91‑1 (MQ=255)
                                                                            gaaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGa               >  1:899067/1‑91 (MQ=255)
                                                                            gaaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGa               <  1:3677736/91‑1 (MQ=255)
                                                                             aaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAg              <  2:2429577/91‑1 (MQ=255)
                                                                              tatttCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt             >  2:619000/5‑91 (MQ=255)
                                                                              aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACGTAAGACAGTAGGAGt             >  1:1446137/1‑91 (MQ=255)
                                                                              aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt             >  1:2714126/1‑91 (MQ=255)
                                                                              aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt             >  1:348070/1‑91 (MQ=255)
                                                                              aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt             >  1:73315/1‑91 (MQ=255)
                                                                              aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt             >  1:1442752/1‑91 (MQ=255)
                                                                              aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt             >  2:177588/1‑91 (MQ=255)
                                                                                tatCTTTCAGTCTGTTAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTACTGTGATTATTTTACCTAAGACAGTAGGAGtta           >  1:3308247/1‑91 (MQ=255)
                                                                                tatCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGtta           >  1:2270237/1‑91 (MQ=255)
                                                                                tatCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGtta           >  2:1306137/1‑91 (MQ=255)
                                                                                 atCTTCCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGttat          <  2:3675059/91‑1 (MQ=255)
                                                                                  tCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGttatt         >  2:1463410/1‑91 (MQ=255)
                                                                                   cTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGttattt        >  2:1593709/1‑91 (MQ=255)
                                                                                     ttCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGttataaat      <  1:1552034/91‑5 (MQ=255)
                                                                                        aGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAAc   <  2:417305/91‑1 (MQ=255)
                                                                                        aGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAAc   <  1:1057695/91‑1 (MQ=255)
                                                                                          |                                                                                         
CAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTATTTATAACT  >  USA300TCH1516_ALE/583870‑584049

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: