Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 583,951 A→T 6.4% 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 6.4% 729.8 / 27.0 217intergenic (+35/‑237)rplA/rplJ50S ribosomal protein L1/50S ribosomal protein L10
Reads supporting (aligned to +/- strand):  ref base A (100/103);  new base T (12/2);  total (112/105)
Fisher's exact test for biased strand distribution p-value = 1.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 227 total aligned reads displayed.

CTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTAT  >  USA300TCH1516_ALE/583861‑584041
                                                                                          |                                                                                          
  gTAACTACAACAATGCGGCCTGGAGTTAAAATTGATACTTCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATc                                                                                          <  2:1645957/91‑1 (MQ=255)
     acggcaacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGGTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCttt                                                                                       <  2:2576813/87‑1 (MQ=255)
         caacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGt                                                                                   >  1:2467442/1‑91 (MQ=255)
         caacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGt                                                                                   >  2:2770020/1‑91 (MQ=255)
         caacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGt                                                                                   <  1:150959/91‑1 (MQ=255)
         caacaaTGGGCCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGt                                                                                   >  2:3409844/1‑91 (MQ=255)
          aacaaTGGGTCCTGGAGTTAAAATTGATACTGCAAGTTACAAATAATAAATGATATAAACAATTACAGGCTGAAACAAATATCTTTCAGTc                                                                                  <  2:737465/91‑1 (MQ=255)
          aacaaTGGGTCCTGGAGTTAAAATTGATACTACAAGTTTCAAATTATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTc                                                                                  <  2:3546826/91‑1 (MQ=255)
               tGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaa                                                                             >  2:1296529/1‑91 (MQ=255)
                gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                            >  2:886863/1‑91 (MQ=255)
                gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                            >  1:1704131/1‑91 (MQ=255)
                gggTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaa                                                                            >  1:3020957/1‑91 (MQ=255)
                 ggTCCTGGAGTTAAAATTAATACTGAAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTaaaa                                                                           <  1:1821451/91‑1 (MQ=255)
                   tCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAat                                                                         <  1:2337823/91‑1 (MQ=255)
                     cTGGAGTTAAAATTGATACTTCAAGTTTCAACTACTAAATGATATAAACAATTACAGGCTGAACGAAATATCTTTCAGTCTGTACAAatat                                                                       >  2:2944416/1‑91 (MQ=255)
                     cTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatat                                                                       <  1:3735302/91‑1 (MQ=255)
                     cTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatat                                                                       <  1:732264/91‑1 (MQ=255)
                      tGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatata                                                                      >  2:926412/1‑91 (MQ=255)
                       ggAGTTATAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatatat                                                                     <  2:1999005/91‑1 (MQ=255)
                       ggAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatatat                                                                     >  1:3825659/1‑91 (MQ=255)
                       ggAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAatatat                                                                     <  2:2838341/91‑1 (MQ=255)
                        gAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATAtt                                                                    >  1:513958/1‑91 (MQ=255)
                         aGTTAAAAGTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTg                                                                   <  1:1958125/91‑1 (MQ=255)
                          gTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGa                                                                  >  2:3700250/1‑91 (MQ=255)
                          gTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGa                                                                  >  2:3859120/1‑91 (MQ=255)
                           ttAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGAc                                                                 >  2:2041476/1‑91 (MQ=255)
                           ttAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGAc                                                                 >  2:2653846/1‑91 (MQ=255)
                               aaTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATa                                                             <  1:2724428/91‑1 (MQ=255)
                                aTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATaa                                                            >  1:363774/1‑91 (MQ=255)
                                aTTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATaa                                                            <  2:841490/91‑1 (MQ=255)
                                   gATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTa                                                         >  2:709218/1‑91 (MQ=255)
                                    aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                        >  2:2154375/1‑91 (MQ=255)
                                    aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                        >  1:2205695/1‑91 (MQ=255)
                                    aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                        >  1:3538604/1‑91 (MQ=255)
                                    aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                        >  1:1624563/1‑91 (MQ=255)
                                    aTACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTaa                                                        >  1:1345463/1‑91 (MQ=255)
                                       cTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAAttt                                                     >  2:485855/1‑91 (MQ=255)
                                         gCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTcc                                                   >  1:2699787/1‑91 (MQ=255)
                                         gCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTcc                                                   <  1:427057/91‑1 (MQ=255)
                                           aaGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCaa                                                 <  1:2077732/91‑1 (MQ=255)
                                            aGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAg                                                <  2:3283190/91‑1 (MQ=255)
                                             ggttCAAATAATAAATGATATAATAAATTACAAGCTGAAAGAAATATCTTACAAACTAACAAAATATATTGACAATAAGTAATTTCCAAGt                                               <  2:3026745/89‑1 (MQ=255)
                                               ttCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATAGCTTTCAGTCTGTAAAAAGATATTGACAATAAGTAATTTCCAAGTta                                             <  2:2294369/91‑1 (MQ=255)
                                                 cAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTtata                                           >  1:3823654/1‑91 (MQ=255)
                                                  aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                          >  2:514799/1‑91 (MQ=255)
                                                  aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                          >  2:2222021/1‑91 (MQ=255)
                                                  aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                          <  1:998032/91‑1 (MQ=255)
                                                  aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                          <  1:1088531/91‑1 (MQ=255)
                                                  aaaTAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTatat                                          <  1:852796/91‑1 (MQ=255)
                                                   aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATAtt                                         >  1:880036/1‑91 (MQ=255)
                                                   aataatAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCACTCTGTAAAAATATATTGACAATAAGTAATCTCCAAGTTATAtt                                         >  1:3954752/1‑91 (MQ=255)
                                                     taataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTAc                                       <  1:1288951/91‑1 (MQ=255)
                                                     taataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTAc                                       <  2:3592904/91‑1 (MQ=255)
                                                       ataaATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACtt                                     >  1:202955/1‑91 (MQ=255)
                                                         aaaTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTAt                                   >  2:1495663/1‑91 (MQ=255)
                                                           aTGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtg                                 <  2:1139872/91‑1 (MQ=255)
                                                             gATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgtg                               <  2:1149989/91‑1 (MQ=255)
                                                             gATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATtgtg                               <  2:1548717/91‑1 (MQ=255)
                                                               tataAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGat                             <  2:2606486/91‑1 (MQ=255)
                                                               tataAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGat                             >  2:1581310/1‑91 (MQ=255)
                                                                ataAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGatt                            <  2:1229931/91‑1 (MQ=255)
                                                                ataAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGatt                            >  1:2327596/1‑91 (MQ=255)
                                                                ataAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGatt                            >  2:1694204/1‑91 (MQ=255)
                                                                 taAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGatta                           <  1:704274/91‑1 (MQ=255)
                                                                   aaCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGattatt                         >  1:3481108/1‑91 (MQ=255)
                                                                   aaCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGattatt                         >  1:536599/1‑91 (MQ=255)
                                                                    aCAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTAttt                        <  1:2406866/91‑1 (MQ=255)
                                                                    aCAATTACAGGCGGAAAGAAATATCTTTCAGGCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTAttt                        <  2:711140/91‑1 (MQ=255)
                                                                        ttACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTAcc                    <  1:1878932/91‑1 (MQ=255)
                                                                         tACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCt                   <  1:1417873/91‑1 (MQ=255)
                                                                         tACAGACTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCt                   >  1:3654148/1‑91 (MQ=37)
                                                                          aCAGGCTGAATGAAATATCTTTCAGCCTGTAAAATCATATTCACAATAAGTAATTTCCAAGTCATATTACTTATTGTGATTATTTGACCTa                  <  2:2373439/91‑1 (MQ=255)
                                                                          aCAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTa                  <  2:2266468/91‑1 (MQ=255)
                                                                          aCAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTa                  >  1:2000610/1‑91 (MQ=255)
                                                                           cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa                 <  2:140318/91‑1 (MQ=255)
                                                                           cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa                 >  2:2811419/1‑91 (MQ=255)
                                                                           cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa                 >  1:2537347/1‑91 (MQ=255)
                                                                           cAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTaa                 <  2:2885500/91‑1 (MQ=255)
                                                                             ggCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa               <  1:2763314/91‑1 (MQ=255)
                                                                             gacTGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGa               >  2:648656/3‑91 (MQ=38)
                                                                              gCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGAc              >  1:2629811/1‑91 (MQ=255)
                                                                               cTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACa             <  1:1178631/91‑1 (MQ=255)
                                                                               cTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACa             <  2:3178985/91‑1 (MQ=255)
                                                                               cTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACa             >  1:1372571/1‑91 (MQ=255)
                                                                               cTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACa             <  1:2875971/91‑1 (MQ=255)
                                                                               cTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGCTATATTACTTATTGTGATTATTTTACCTAAGACa             >  1:2167328/1‑91 (MQ=255)
                                                                                tGAGAGAAATATCTTTTAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAg            <  1:3182818/91‑1 (MQ=255)
                                                                                tGAAAGATATTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAg            >  1:1192963/1‑91 (MQ=38)
                                                                                 gaaagataTTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGt           >  2:3347368/1‑91 (MQ=39)
                                                                                 gaaagataTTTCTTTCAGCCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGt           >  1:2483064/1‑91 (MQ=39)
                                                                                  aaagaaaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGGGATTGCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTa          <  2:3810919/91‑1 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGCTTATTTTACCTAAGAACGTAGGAGt     >  2:2736657/1‑91 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt     <  1:1143643/91‑1 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt     >  1:2298679/1‑91 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt     >  2:234722/1‑91 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt     <  1:3380290/91‑1 (MQ=255)
                                                                                       aaTATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGt     <  1:918647/91‑1 (MQ=255)
                                                                                        atatCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGtt    >  1:3068810/1‑91 (MQ=255)
                                                                                          atCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGttat  >  2:2878340/1‑91 (MQ=255)
                                                                                          atCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGttat  >  1:3072717/1‑91 (MQ=255)
                                                                                          |                                                                                          
CTGTAACTACAACAATGGGTCCTGGAGTTAAAATTGATACTGCAAGTTTCAAATAATAAATGATATAAACAATTACAGGCTGAAAGAAATATCTTTCAGTCTGTAAAAATATATTGACAATAAGTAATTTCCAAGTTATATTACTTATTGTGATTATTTTACCTAAGACAGTAGGAGTTAT  >  USA300TCH1516_ALE/583861‑584041

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