Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,722,313 G→A 7.8% Q3* (CAA→TAA)  yrrK ← Putative pre‑16S rRNA nuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,722,3130GA 7.8% 644.1 / 30.7 207Q3* (CAA→TAA) yrrKPutative pre‑16S rRNA nuclease
Reads supporting (aligned to +/- strand):  ref base G (97/93);  new base A (4/12);  total (102/105)
Fisher's exact test for biased strand distribution p-value = 6.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

Only 100 of 213 total aligned reads displayed.

TGTCTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGATTTCT  >  USA300TCH1516_ALE/1722223‑1722403
                                                                                          |                                                                                          
  tCTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTa                                                                                          <  2:3725571/91‑1 (MQ=255)
   cTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTaa                                                                                         >  1:2516109/1‑91 (MQ=255)
   cTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTaa                                                                                         >  1:2771865/1‑91 (MQ=255)
   cTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTaa                                                                                         >  2:2674775/1‑91 (MQ=255)
   cTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTaa                                                                                         >  1:1766334/1‑91 (MQ=255)
   cTAATCCATGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTaa                                                                                         <  2:3454968/91‑1 (MQ=255)
    tAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAAc                                                                                        <  2:1717837/91‑1 (MQ=255)
       tCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATa                                                                                     <  1:2487109/91‑1 (MQ=255)
       tCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATa                                                                                     <  1:3391880/91‑1 (MQ=255)
       tCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATa                                                                                     <  2:1877058/91‑1 (MQ=255)
       tCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATa                                                                                     <  2:3339245/91‑1 (MQ=255)
       tCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATa                                                                                     <  1:1050828/91‑1 (MQ=255)
         cTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATaaa                                                                                   >  1:714599/1‑91 (MQ=255)
           tgtgCCGTCCAGCCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATaaaat                                                                                 <  1:3518774/91‑3 (MQ=255)
           tgtgCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATaaa                                                                                   >  2:127975/1‑89 (MQ=255)
           tgtgCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAatt                                                                                 >  1:919663/1‑91 (MQ=255)
           tgtgCCGTCCAACCCATATATCCCCAATTGCCATTCCCACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTTACATAAattt                                                                                >  1:1001620/1‑90 (MQ=255)
            gtgCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAatta                                                                                >  2:2249548/1‑91 (MQ=255)
            gtgCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAatta                                                                                <  2:2879012/91‑1 (MQ=255)
            gtgCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAatta                                                                                >  1:744518/1‑91 (MQ=255)
               ccGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAAttattt                                                                             >  2:803484/1‑91 (MQ=255)
               ccGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAAttattt                                                                             >  1:2115879/1‑91 (MQ=255)
                cGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAAttattta                                                                            <  1:2139141/91‑1 (MQ=255)
                  tCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATtatttatt                                                                          <  2:1011386/91‑1 (MQ=255)
                  tCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATtatttatt                                                                          >  2:1215894/1‑91 (MQ=255)
                   ccAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTatttattt                                                                         <  1:2323709/91‑1 (MQ=255)
                   ccAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTatttattt                                                                         <  1:3682491/91‑1 (MQ=255)
                       cccATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGct                                                                     >  1:2839204/1‑91 (MQ=255)
                        ccATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGctc                                                                    <  1:2596637/91‑1 (MQ=255)
                         cATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGctct                                                                   <  2:1438465/91‑1 (MQ=255)
                          aTTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCtt                                                                  >  1:2864680/1‑91 (MQ=255)
                           ttATATCACTAATTGTAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTCACATAAATTATTTATTTTGCTCttt                                                                 >  1:3627950/1‑91 (MQ=255)
                            tataTCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCtttt                                                                <  2:3583919/91‑1 (MQ=255)
                            tataTCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCtttt                                                                <  2:1662716/91‑1 (MQ=255)
                             atatCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTa                                                               <  1:3142816/91‑1 (MQ=255)
                             atatCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTa                                                               <  1:1216939/91‑1 (MQ=255)
                               atCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTaaa                                                             >  2:2771338/1‑91 (MQ=255)
                                tCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAAt                                                            <  2:3871072/91‑1 (MQ=255)
                                  aCTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAg                                                          >  1:3752672/1‑91 (MQ=255)
                                   cTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTCTTTTGCTCTTTTAAATAGt                                                         >  2:2204990/1‑91 (MQ=255)
                                   cTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGt                                                         <  1:1864251/91‑1 (MQ=255)
                                   cTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGt                                                         >  2:3879756/1‑91 (MQ=255)
                                     aaTTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTaa                                                       <  1:1707964/91‑1 (MQ=255)
                                     aaTTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTaa                                                       <  2:1088207/91‑1 (MQ=255)
                                     aaTTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTaa                                                       <  2:3675504/91‑1 (MQ=255)
                                     aaTTGCAAGTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTaa                                                       <  2:2937456/91‑1 (MQ=255)
                                       ttGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGa                                                     >  1:166081/1‑91 (MQ=255)
                                       ttGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGa                                                     >  1:2462978/1‑91 (MQ=255)
                                       ttGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGa                                                     >  1:2464091/1‑91 (MQ=255)
                                         gCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGaaa                                                   <  1:2229399/91‑1 (MQ=255)
                                          cAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAAc                                                  >  2:1658146/1‑91 (MQ=255)
                                           atttCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACa                                                 >  1:2082251/3‑91 (MQ=255)
                                           aaTTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACa                                                 <  1:58864/91‑1 (MQ=255)
                                           aaTTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACa                                                 <  1:1809107/91‑1 (MQ=255)
                                           aaTTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATGTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACa                                                 <  1:3500148/91‑1 (MQ=255)
                                            aTTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACaa                                                >  1:258010/1‑91 (MQ=255)
                                            aTTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACaa                                                <  1:1330692/91‑1 (MQ=255)
                                             ttCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAg                                               >  2:2702912/1‑91 (MQ=255)
                                             ttCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAg                                               <  1:321093/91‑1 (MQ=255)
                                               ccTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTCATTTTGCTCTTTTAAATAGTAAGAAACAAGct                                             <  2:3804259/91‑1 (MQ=255)
                                                  aCCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCtt                                          >  1:1474744/1‑91 (MQ=255)
                                                   ccGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTc                                         >  1:1164308/1‑91 (MQ=255)
                                                     gTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCa                                       >  2:243362/1‑91 (MQ=255)
                                                      ttCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCat                                      <  1:1580811/91‑1 (MQ=255)
                                                      ttCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCat                                      <  2:2043596/91‑1 (MQ=255)
                                                       tCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCata                                     <  2:512087/91‑1 (MQ=255)
                                                       tCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCata                                     <  1:634343/91‑1 (MQ=255)
                                                        ctactaCCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCataa                                    >  1:1038960/1‑91 (MQ=255)
                                                        ctactaCCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCataa                                    <  2:1292906/91‑1 (MQ=255)
                                                          actacCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCataata                                  >  2:698873/1‑91 (MQ=255)
                                                          actacCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCataata                                  >  2:698712/1‑91 (MQ=255)
                                                           ctacCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCAtaataa                                 <  1:553651/91‑1 (MQ=255)
                                                           ctacCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCAtaataa                                 <  2:3223075/91‑1 (MQ=255)
                                                            tacCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAAc                                <  1:2928741/91‑1 (MQ=255)
                                                            tacCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAAc                                >  2:2373481/1‑91 (MQ=255)
                                                             acCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACa                               >  2:1982188/1‑91 (MQ=255)
                                                              ccGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACAt                              >  2:370679/1‑91 (MQ=255)
                                                               cGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACAtc                             <  1:991443/91‑1 (MQ=255)
                                                                 aCATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACAtctc                           >  1:378405/1‑91 (MQ=255)
                                                                 aCATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACAtctc                           <  2:2398410/91‑1 (MQ=255)
                                                                   aTCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTa                         >  2:1752654/1‑91 (MQ=255)
                                                                    tCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTAt                        <  1:252438/91‑1 (MQ=255)
                                                                     cGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTAtt                       >  2:1048912/1‑90 (MQ=255)
                                                                      gAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCa                      <  1:502997/91‑1 (MQ=255)
                                                                       aGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCaa                     >  1:3774936/1‑91 (MQ=255)
                                                                       aGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCaa                     >  2:2314466/1‑91 (MQ=255)
                                                                        gTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGc                                              <  1:3867150/65‑1 (MQ=255)
                                                                        gTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAat                    >  1:3006196/1‑91 (MQ=255)
                                                                        gTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAat                    >  1:44818/1‑91 (MQ=255)
                                                                          ccTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAGTAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAatat                  <  2:2995826/91‑1 (MQ=255)
                                                                               taatttATGTTACAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGa             >  1:495170/4‑91 (MQ=255)
                                                                                aatttATGTTACAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGaa            >  2:2569068/3‑91 (MQ=255)
                                                                                aTTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGaa            <  2:1137485/91‑1 (MQ=255)
                                                                                  tttATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAAtt          <  1:2429405/91‑1 (MQ=255)
                                                                                  tttATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAAtt          >  1:719422/1‑91 (MQ=255)
                                                                                  tttATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAAtt          <  1:429654/91‑1 (MQ=255)
                                                                                   ttATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAAttt         <  2:2839204/91‑1 (MQ=255)
                                                                                    tATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTg        <  2:575518/91‑1 (MQ=255)
                                                                                      tgttgtAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGAt      >  1:2666271/1‑91 (MQ=255)
                                                                                       gttgtAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGAtt     <  2:1782220/91‑1 (MQ=255)
                                                                                          |                                                                                          
TGTCTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGATTTCT  >  USA300TCH1516_ALE/1722223‑1722403

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: