Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 583,948 A→T 5.7% intergenic (+32/‑240) 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,9480AT 5.7% 651.0 / 19.3 194intergenic (+32/‑240)rplA/rplJ50S ribosomal protein L1/50S ribosomal protein L10
Reads supporting (aligned to +/- strand):  ref base A (90/93);  new base T (8/3);  total (98/96)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 223 total aligned reads displayed.

TTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGT  >  USA300TCH1516_ALE/583858‑584038
                                                                                          |                                                                                          
 cggtGTAACTACAATAATGGGTCCTGGAGTTCACATTGATTGCTCGAGTTTCAAGTAATAACTGATATATACCGTCACCGGCTgcgagaac                                                                                           <  2:2440146/88‑2 (MQ=255)
     gTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATc                                                                                       <  1:708829/91‑1 (MQ=255)
     gTAACTACAACAATGCGGCCTGGAGTTAAAATTGATACTTCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATc                                                                                       <  2:1645957/91‑1 (MQ=255)
       aaCTACAACAATGGGTCCTGGAGTTAAGATTGATACTGCGAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCtt                                                                                     >  1:264223/1‑91 (MQ=255)
         cTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTc                                                                                   <  1:2164589/91‑1 (MQ=255)
         cTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGACTGAAAGATAtttctttc                                                                                   <  1:2713084/91‑1 (MQ=39)
           acaacaATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATACTAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAg                                                                                 >  1:3952430/1‑91 (MQ=255)
           acaacaATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAg                                                                                 >  2:1495934/1‑91 (MQ=255)
            caacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGt                                                                                <  1:2982842/91‑1 (MQ=255)
            caacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGt                                                                                >  2:2770020/1‑91 (MQ=255)
             aacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTc                                                                               <  2:738573/91‑1 (MQ=255)
             aacaaTGGGTCCTGGAGTTAAAATTGATACTACAAGTTTCAAATTATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTc                                                                               <  2:3546826/91‑1 (MQ=255)
                aaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGt                                                                            <  1:2016564/91‑1 (MQ=255)
                  tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaa                                                                          >  2:1296529/1‑91 (MQ=255)
                   gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                         >  1:1704131/1‑91 (MQ=255)
                   gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                         >  1:3021174/1‑91 (MQ=255)
                   gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                         >  2:1557348/1‑91 (MQ=255)
                   gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                         >  2:886863/1‑91 (MQ=255)
                    ggTCCTGGAGTTAAAATTAATACTGAAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaa                                                                        <  1:1821451/91‑1 (MQ=255)
                    ggTCCTGGAGTTAAAAGTGATACTGCAAGTTTCAAATAATAAATGAGATAAACAATTACAGGCTGAAAGAAATATCTTACAGTCTGTaaaa                                                                        <  2:2650127/91‑1 (MQ=255)
                     gTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                                       <  1:1835397/91‑1 (MQ=255)
                     gTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaaa                                                                       <  2:2973085/91‑1 (MQ=255)
                        cTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatat                                                                    <  1:732481/91‑1 (MQ=255)
                        cTGGAGTTAAAATAGATACTGCAAGTTTCAAATAATAAAAGATATAAACAATTACAGGCTGAAAGAAATATCTTAAAGTCTATAAAAatat                                                                    <  2:3711398/91‑1 (MQ=255)
                         tGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatata                                                                   >  2:926412/1‑91 (MQ=255)
                          ggAGTTATAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatatat                                                                  <  2:1999033/91‑1 (MQ=255)
                          ggAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatatat                                                                  >  1:3825647/1‑91 (MQ=255)
                          ggAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatatat                                                                  <  2:2838341/91‑1 (MQ=255)
                           gAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATAtt                                                                 >  1:35902/1‑91 (MQ=255)
                             gTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGa                                                               >  2:3859120/1‑91 (MQ=255)
                                   aTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATaa                                                         <  2:841490/91‑1 (MQ=255)
                                     tGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGt                                                       <  1:2471329/91‑1 (MQ=255)
                                     tGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGt                                                       <  1:3708493/91‑1 (MQ=255)
                                      gATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTa                                                      >  2:709218/1‑91 (MQ=255)
                                      gATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTa                                                      >  1:2730652/1‑91 (MQ=255)
                                       aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                     >  1:2205695/1‑91 (MQ=255)
                                       aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                     >  2:2154375/1‑91 (MQ=255)
                                       aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                     >  2:564397/1‑91 (MQ=255)
                                       aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                     >  1:1345463/1‑91 (MQ=255)
                                       aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                     >  1:1624563/1‑91 (MQ=255)
                                        tACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAAt                                                    <  2:3327090/91‑1 (MQ=255)
                                         aCTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAAtt                                                   >  1:3806098/1‑91 (MQ=255)
                                         aCTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAAtt                                                   >  2:1245235/1‑91 (MQ=255)
                                          cTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAAttt                                                  >  2:485855/1‑91 (MQ=255)
                                            gCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTcc                                                >  2:3365256/1‑91 (MQ=255)
                                              aaGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCaa                                              <  1:2077491/91‑1 (MQ=255)
                                               aGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAg                                             <  1:232180/91‑1 (MQ=255)
                                               aGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAg                                             <  2:1221400/91‑1 (MQ=255)
                                                gTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGt                                            <  2:2053063/91‑1 (MQ=255)
                                                  ttCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATAGCTTTCAGTCTGTAAAAAGATATTGACAATAAGTAATTTCCAAGTta                                          <  2:2294369/91‑1 (MQ=255)
                                                    cAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTtata                                        <  1:147006/91‑1 (MQ=255)
                                                     aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                       <  1:998032/91‑1 (MQ=255)
                                                     aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                       <  1:1088531/91‑1 (MQ=255)
                                                     aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                       <  1:907972/91‑1 (MQ=255)
                                                     aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                       <  1:1922583/91‑1 (MQ=255)
                                                     aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                       >  2:2222934/1‑91 (MQ=255)
                                                      aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATAtt                                      >  1:878961/1‑91 (MQ=255)
                                                            aaaTGATATAAACAATTACAGGCTGAAAGAAATATCATTCAGTCTGTGAAAATATATTGACAATAATCAATTTCCACGTTATATTACTTAt                                >  2:3673738/1‑91 (MQ=255)
                                                             aaTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTAtt                               >  1:128644/1‑91 (MQ=255)
                                                              aTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtg                              <  2:1139872/91‑1 (MQ=255)
                                                              aTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtg                              >  1:536753/1‑91 (MQ=255)
                                                              aTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtg                              <  2:1827330/91‑1 (MQ=255)
                                                                  tataAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGat                          >  1:929757/1‑91 (MQ=255)
                                                                  tataAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGat                          >  2:1581310/1‑91 (MQ=255)
                                                                   ataAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGatt                         >  1:2327646/1‑91 (MQ=255)
                                                                   ataAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGatt                         <  2:199718/91‑1 (MQ=255)
                                                                       aCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTAttt                     <  1:2406866/91‑1 (MQ=255)
                                                                       aCAATTACAGGCGGAAAGAAATATCTTTCAGGCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTAttt                     <  2:711140/91‑1 (MQ=255)
                                                                         aaTTACAGGCTGACAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCACGTTATATTACTTATTGTGATTCTTTTa                   >  2:2900372/1‑91 (MQ=255)
                                                                         aaTTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTa                   <  1:3490408/91‑1 (MQ=255)
                                                                          aTTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTAc                  <  1:2562247/91‑1 (MQ=255)
                                                                          aTTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTAc                  >  2:2589932/1‑91 (MQ=255)
                                                                            tACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCt                <  1:1417873/91‑1 (MQ=255)
                                                                            tACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCt                >  1:2522618/1‑91 (MQ=255)
                                                                            tACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCt                <  1:530248/91‑1 (MQ=255)
                                                                             aCAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTa               >  1:2000610/1‑91 (MQ=255)
                                                                             aCAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTa               <  2:2266468/91‑1 (MQ=255)
                                                                              cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa              >  2:2811419/1‑91 (MQ=255)
                                                                              cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa              <  2:140318/91‑1 (MQ=255)
                                                                                ggCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa            >  1:2444600/1‑91 (MQ=255)
                                                                                gacTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa            >  2:648656/3‑91 (MQ=38)
                                                                                  cTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACa          <  1:1178631/91‑1 (MQ=255)
                                                                                  cTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACa          <  2:3178985/91‑1 (MQ=255)
                                                                                  cTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACa          >  1:1372571/1‑91 (MQ=255)
                                                                                   tGAGAGAAATATCTTTTAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAg         <  1:3182818/91‑1 (MQ=255)
                                                                                   tGAAAGATATTTCTTTCAGCCTGTAACAATATATTGACAATAAGTCATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAg         >  2:2853953/1‑91 (MQ=38)
                                                                                   tGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAg         >  1:1192963/1‑91 (MQ=38)
                                                                                   tGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAg         <  1:489718/91‑1 (MQ=255)
                                                                                    gaaagataTTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGt        >  1:2778896/1‑91 (MQ=39)
                                                                                    gaaagaaaTATCTCTCAGTCTGTAAAAATATATTGACAATAAGGAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGt        <  1:3955634/91‑1 (MQ=255)
                                                                                     aaagaaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGGGATTGCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTa       <  2:3810919/91‑1 (MQ=255)
                                                                                      aagataTTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAg      >  2:307708/1‑91 (MQ=255)
                                                                                      aagaaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAg      <  2:2430917/91‑1 (MQ=255)
                                                                                        gaaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGa    <  1:3646452/91‑1 (MQ=255)
                                                                                          aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGCTTATTTTACCTAAGAACGTAGGAGt  >  2:2736657/1‑91 (MQ=255)
                                                                                          aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt  >  1:2282273/1‑91 (MQ=255)
                                                                                          aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt  <  1:3380290/91‑1 (MQ=255)
                                                                                          aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt  >  2:234722/1‑91 (MQ=255)
                                                                                          aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt  <  1:1143643/91‑1 (MQ=255)
                                                                                          aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTACTTATTTTACCTAAGACAGTCGGAGt  >  2:3422778/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTGCTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGT  >  USA300TCH1516_ALE/583858‑584038

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