Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,722,313 G→A 9.2% 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 9.2% 333.6 / 18.5 109Q3* (CAA→TAA) yrrKPutative pre‑16S rRNA nuclease
Reads supporting (aligned to +/- strand):  ref base G (52/47);  new base A (2/8);  total (54/55)
Fisher's exact test for biased strand distribution p-value = 9.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 123 total aligned reads displayed.

TGTCTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGATTTCT  >  USA300TCH1516_ALE/1722223‑1722403
                                                                                          |                                                                                          
tgtCTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTAtgttg                                                                                            >  1:388391/1‑91 (MQ=255)
tgtCTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTAtgttg                                                                                            >  1:483864/1‑91 (MQ=255)
   cTAATCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTaa                                                                                         <  1:1058113/91‑1 (MQ=255)
      aTCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACAt                                                                                      <  1:3036609/91‑1 (MQ=255)
      aTCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACAt                                                                                      <  1:3227908/91‑1 (MQ=255)
       tCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATa                                                                                     <  1:84693/91‑1 (MQ=255)
       tCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATa                                                                                     <  1:84470/91‑1 (MQ=255)
       tCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATa                                                                                     <  1:3008453/91‑1 (MQ=255)
       tCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATa                                                                                     <  1:373539/91‑1 (MQ=255)
       tCCTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATa                                                                                     <  1:3441549/91‑1 (MQ=255)
        ccTTGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTACAACATaa                                                                                    <  1:2622365/91‑1 (MQ=255)
          ttGTGCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAA‑TTTATGTTACAACATaaaat                                                                                 <  1:1788416/91‑3 (MQ=255)
           tgtgCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAatt                                                                                 >  1:1261700/1‑91 (MQ=255)
            gtgCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAatta                                                                                <  1:692131/91‑1 (MQ=255)
             tgCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTAGTAACATAAagtat                                                                               >  1:3740530/1‑91 (MQ=255)
              gCCGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAattatt                                                                              >  1:2199089/1‑91 (MQ=255)
               ccGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAAttattt                                                                             <  1:1130031/91‑1 (MQ=255)
               ccGTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAAttattt                                                                             <  1:251592/91‑1 (MQ=255)
                 gTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAAttatttat                                                                           >  1:3726679/1‑91 (MQ=255)
                 gTCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAAttatttat                                                                           <  1:2129454/91‑1 (MQ=255)
                  tCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATtatttatt                                                                          <  1:644701/91‑1 (MQ=255)
                  tCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCATAAAATTTTATGTTGTAACATAAATtatttatt                                                                          >  1:3677974/1‑91 (MQ=255)
                  tCCAACCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCATAAAATTTTATGTTGTAACATAAATtatttatt                                                                          >  1:3678511/1‑91 (MQ=255)
                   ccAACCCATTATATCCCTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTatttattt                                                                         >  1:3522385/1‑91 (MQ=255)
                     aaCCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTg                                                                       >  1:3380735/1‑91 (MQ=255)
                      aCCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGc                                                                      <  1:3840869/91‑1 (MQ=255)
                      aCCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGc                                                                      <  1:3840147/91‑1 (MQ=255)
                      aCCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGc                                                                      >  1:1890125/1‑91 (MQ=255)
                      aCCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGc                                                                      <  1:3840254/91‑1 (MQ=255)
                      aCCCATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGc                                                                      <  1:3840290/91‑1 (MQ=255)
                       cccATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGct                                                                     <  1:2582962/91‑1 (MQ=255)
                        ccATTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGctc                                                                    >  1:4117461/1‑91 (MQ=255)
                          aTTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCtt                                                                  >  1:24770/1‑91 (MQ=255)
                          aTTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCtt                                                                  <  1:849651/91‑1 (MQ=255)
                          aTTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCtt                                                                  <  1:3659386/91‑1 (MQ=255)
                          aTTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCtt                                                                  >  1:679949/1‑91 (MQ=255)
                          aTTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCtt                                                                  <  1:1096577/91‑1 (MQ=255)
                          aTTATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATCTATTTTGCTCtt                                                                  <  1:3188521/91‑1 (MQ=255)
                           ttATATCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCttt                                                                 >  1:3390435/1‑91 (MQ=255)
                            tataTCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCtttt                                                                <  1:2214562/91‑1 (MQ=255)
                              tatCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTaa                                                              >  1:336639/1‑91 (MQ=255)
                              tatCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTACTTATTTTGCTCTTTTaa                                                              >  1:3557363/1‑91 (MQ=255)
                                tCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAACATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAAt                                                            >  1:3769468/1‑91 (MQ=255)
                                tCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTCTGTTGTAACATAAATTATTTATTTTGCTCTTTTAAAt                                                            <  1:3506248/91‑1 (MQ=255)
                                tCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAAt                                                            >  1:1109192/1‑91 (MQ=255)
                                tCACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAAt                                                            >  1:3770605/1‑91 (MQ=255)
                                 cACTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATa                                                           <  1:3252984/91‑1 (MQ=255)
                                   cTAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGt                                                         >  1:518730/1‑91 (MQ=255)
                                    tAATTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTa                                                        <  1:890269/91‑1 (MQ=255)
                                      aTTGCAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAg                                                      >  1:476823/1‑91 (MQ=255)
                                          cAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAAc                                                  <  1:3372665/91‑1 (MQ=255)
                                          cAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAAc                                                  >  1:1059128/1‑91 (MQ=255)
                                          cAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAAc                                                  <  1:2343601/91‑1 (MQ=255)
                                          cAATTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAAc                                                  <  1:2342692/91‑1 (MQ=255)
                                            aTTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACaa                                                >  1:1826281/1‑91 (MQ=255)
                                            aTTCCTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACaa                                                >  1:1826430/1‑91 (MQ=255)
                                               ccTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGct                                             >  1:1271914/1‑91 (MQ=255)
                                               ccTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGct                                             <  1:1623417/91‑1 (MQ=255)
                                               ccTACCGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGct                                             >  1:2297812/1‑91 (MQ=255)
                                                  aCCGTTCTACTACCGACGTCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCtt                                          <  1:3559459/91‑1 (MQ=255)
                                                   ccGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTc                                         <  1:1995917/91‑1 (MQ=255)
                                                    cGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTcc                                        >  1:1223351/1‑91 (MQ=255)
                                                    cGTTCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTcc                                        >  1:3991366/1‑91 (MQ=255)
                                                       tCTACTACCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCata                                     <  1:2417356/91‑1 (MQ=255)
                                                         tactacCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCataat                                   <  1:1084514/91‑1 (MQ=255)
                                                         tactacCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCataat                                   <  1:61909/91‑1 (MQ=255)
                                                         tactacCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCataat                                   >  1:1912239/1‑91 (MQ=255)
                                                          actacCGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCataata                                  >  1:1605041/1‑91 (MQ=255)
                                                           ctacCGACATCGAGTCCTAAAATCTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCAtaataa                                 >  1:587417/1‑91 (MQ=255)
                                                               cGACATCGAGTCNTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACAtc                             <  1:48469/91‑1 (MQ=255)
                                                               cGACATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACAtc                             >  1:1625660/1‑91 (MQ=255)
                                                                 aCATCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACAtctc                           <  1:3029218/91‑1 (MQ=255)
                                                                   aTCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTa                         >  1:3261194/1‑91 (MQ=255)
                                                                   aTCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTa                         >  1:4049671/1‑91 (MQ=255)
                                                                    tCGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTAt                        >  1:735204/1‑91 (MQ=255)
                                                                     cGAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATc                       >  1:692972/1‑91 (MQ=255)
                                                                      gAGTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCa                      <  1:1940801/91‑1 (MQ=255)
                                                                        gTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAat                    >  1:2762351/1‑91 (MQ=255)
                                                                        gTCCTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAat                    >  1:1452857/1‑91 (MQ=255)
                                                                           cTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATg                 <  1:1767350/91‑1 (MQ=255)
                                                                           cTAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATg                 <  1:2159921/91‑1 (MQ=255)
                                                                            tAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGa                <  1:1733902/91‑1 (MQ=255)
                                                                            tAAAATTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGa                >  1:401406/1‑91 (MQ=255)
                                                                              aaaTTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACg              >  1:730512/1‑91 (MQ=255)
                                                                               taatttATGTTACAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGa             >  1:2999979/4‑91 (MQ=255)
                                                                               aaTTTTATGTTGCAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGa             >  1:3211553/1‑91 (MQ=255)
                                                                                aatttATGTTACAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGaa            >  1:236081/3‑91 (MQ=255)
                                                                                aTTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGaa            >  1:1439557/1‑91 (MQ=255)
                                                                                aTTTTATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGaa            >  1:1439556/1‑91 (MQ=255)
                                                                                  tttATGTTGTAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAAtt          <  1:3791722/91‑1 (MQ=255)
                                                                                     aTGTTGTAACATCAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGTTCTTCCATAATAACATCTCTATCAATATGACGAATGTGa       >  1:3207258/1‑91 (MQ=255)
                                                                                      tgttgtAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGAt      <  1:365605/91‑1 (MQ=255)
                                                                                      tgttgtAACATAAATTATTTATTTTACTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGAt      <  1:3873332/91‑1 (MQ=255)
                                                                                       gttgtAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGAtt     >  1:3100464/1‑91 (MQ=255)
                                                                                       gttgtAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGAtt     >  1:1415354/1‑91 (MQ=255)
                                                                                        ttgtAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGAttt    >  1:1388158/1‑91 (MQ=255)
                                                                                          gtAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGATTTCt  >  1:2216835/1‑91 (MQ=255)
                                                                                          gtAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGATTTCt  <  1:20498/91‑1 (MQ=255)
                                                                                          gtAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGATTTCt  >  1:2145923/1‑91 (MQ=255)
                                                                                          gtAACATAAATTATTTATTTTGCTCTTTTAAATAGTAAGAAACAAGCTCTTCCATAATAACATCTCTATCAATATGACGAATTTGATTTCt  >  1:2146207/1‑91 (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: ‑