Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 631,632:1 +T 5.0% intergenic (+258/‑244) gph_1 → / → proP Phosphoglycolate phosphatase/Proline/betaine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE631,6321.T 5.0% 953.8 / 38.5 218intergenic (+258/‑244)gph_1/proPPhosphoglycolate phosphatase/Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base . (120/87);  new base T (4/7);  total (124/94)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 224 total aligned reads displayed.

TAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTAT  >  USA300TCH1516_ALE/631542‑631722
                                                                                           |                                                                                          
 aaatTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTCAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGa                                                                                           <  1:3359748/91‑1 (MQ=255)
  aatTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGaa                                                                                          >  1:1591519/1‑91 (MQ=255)
   atTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAg                                                                                         >  1:940654/1‑91 (MQ=255)
     tAAAGTAATGACGAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAACTGGCGATAGCATATCACTTACATGAAGtt                                                                                       <  2:3384381/91‑1 (MQ=255)
     tAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGtt                                                                                       >  2:2985805/1‑91 (MQ=255)
     tAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGtt                                                                                       >  2:3448687/1‑91 (MQ=255)
     tAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGtt                                                                                       <  1:1917237/91‑1 (MQ=255)
             tGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTa                                                                               <  2:2683832/91‑1 (MQ=255)
              gACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATAGCGATAGCACATAACTTCATGTAAGTGATATGCTAt                                                                              <  1:1254360/91‑1 (MQ=17)
                 aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGc                                                                           <  1:2419875/91‑1 (MQ=255)
                 aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCACATAACTTCATGTAAGTGATATGCTATCGc                                                                           <  1:68279/91‑1 (MQ=18)
                 aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATAGCGATAGCACATAACTTCATGTAAGTGATATGCTATCGc                                                                           <  1:2877287/91‑1 (MQ=17)
                 aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATAGCGATAGCACATAACTTCATGTAAGTGATATGCTATCGc                                                                           <  2:2750245/91‑1 (MQ=17)
                 aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATAGCGATAGCACATAACTTCATGTAAGTGATATGCTATCGc                                                                           <  1:2179083/91‑1 (MQ=17)
                  aaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCt                                                                          <  2:852176/91‑1 (MQ=255)
                  aaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCt                                                                          >  2:3526667/1‑91 (MQ=255)
                  aaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCt                                                                          >  2:2366172/1‑91 (MQ=255)
                   aGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTa                                                                         >  1:800495/1‑91 (MQ=255)
                   aGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTa                                                                         >  1:3373012/1‑91 (MQ=255)
                    gTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAt                                                                        <  1:1177163/91‑1 (MQ=255)
                    gTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAt                                                                        >  1:1886403/1‑91 (MQ=255)
                     ttATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAtt                                                                       >  2:1944332/1‑91 (MQ=255)
                       aTTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAtttt                                                                     <  1:519016/91‑1 (MQ=255)
                        ttGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAttttt                                                                    <  1:2199720/91‑1 (MQ=255)
                        ttGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAttttt                                                                    >  2:2405134/1‑91 (MQ=255)
                        ttGAAGTAAATTGAGTATAAACATTTAAATACCATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAttttt                                                                    >  2:3808757/1‑91 (MQ=255)
                          gAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAg                                                                  <  1:3702464/91‑1 (MQ=255)
                          gAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAg                                                                  <  1:2727695/91‑1 (MQ=255)
                          gAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAg                                                                  >  1:1478363/1‑91 (MQ=255)
                            aGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGtt                                                                >  1:473958/1‑91 (MQ=255)
                              tAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACAAGAAGTTATGTGCTATCGCTATTTTTAGTtat                                                              <  1:3804/91‑1 (MQ=255)
                               aaaTTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTtata                                                             <  1:2908465/91‑1 (MQ=255)
                                 aTTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAAt                                                           >  1:1369632/1‑91 (MQ=255)
                                  ttGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAAtt                                                          >  2:1357459/1‑91 (MQ=255)
                                    gAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTcc                                                        >  1:2693416/1‑91 (MQ=255)
                                    gAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTcc                                                        >  1:2072892/1‑91 (MQ=255)
                                      gTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaa                                                      >  2:1121519/1‑91 (MQ=255)
                                      gTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaa                                                      >  2:3165107/1‑91 (MQ=255)
                                        ataAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaaaa                                                    >  2:2117731/1‑91 (MQ=255)
                                        ataAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaaaa                                                    >  1:719857/1‑91 (MQ=255)
                                        ataAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATCTTTAGTTATAATTCCaaaa                                                    >  2:1444766/1‑91 (MQ=255)
                                          aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                                  <  1:2985805/91‑1 (MQ=255)
                                          aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                                  >  2:194837/1‑91 (MQ=255)
                                          aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                                  <  2:1711941/91‑1 (MQ=255)
                                          aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                                  <  1:757511/91‑1 (MQ=255)
                                           aaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGt                                                 >  1:1520389/1‑91 (MQ=255)
                                             cATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTa                                               >  2:3240267/1‑91 (MQ=255)
                                                 tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                                           <  2:2512522/91‑1 (MQ=255)
                                                 tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                                           <  1:971938/91‑1 (MQ=255)
                                                 tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                                           <  1:973537/91‑1 (MQ=255)
                                                 tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                                           <  1:3676169/91‑1 (MQ=255)
                                                 tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                                           <  2:607266/91‑1 (MQ=255)
                                                  aaaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGt                                          >  1:3593548/1‑91 (MQ=255)
                                                   aaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGtt                                         >  2:1825566/1‑91 (MQ=255)
                                                   aaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGtt                                         >  1:2389354/1‑91 (MQ=255)
                                                   aaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGtt                                         >  1:2339315/1‑91 (MQ=255)
                                                    aTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTc                                        >  1:1648801/1‑91 (MQ=255)
                                                    aTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTc                                        >  1:2570616/1‑91 (MQ=255)
                                                     tACGATGTCGAAAATGGCGATAGCATATCACTTTCATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCg                                       >  2:2831215/1‑91 (MQ=255)
                                                     tACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCg                                       <  1:2092592/91‑1 (MQ=255)
                                                     tACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCg                                       >  1:3354962/1‑91 (MQ=255)
                                                     tACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCg                                       <  1:608165/91‑1 (MQ=255)
                                                       cGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCgat                                     >  1:472513/1‑91 (MQ=255)
                                                         aTGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCgatga                                   <  1:1693055/91‑1 (MQ=255)
                                                           gTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGAtt                                 <  2:1121260/91‑1 (MQ=255)
                                                           gTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGAtt                                 >  2:2061046/1‑91 (MQ=255)
                                                            tCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGAttt                                >  1:1317241/1‑91 (MQ=255)
                                                             cGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTa                               >  1:3599115/1‑91 (MQ=255)
                                                               aaaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAg                             >  1:1407821/1‑91 (MQ=255)
                                                               aaaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAg                             <  2:3404459/91‑1 (MQ=255)
                                                                aaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                            >  2:2290256/1‑91 (MQ=255)
                                                                aaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                            <  2:162735/91‑1 (MQ=255)
                                                                aaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                            <  2:2122767/91‑1 (MQ=255)
                                                                aaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                            <  2:705628/91‑1 (MQ=255)
                                                                aaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                            <  2:160288/91‑1 (MQ=255)
                                                                 aaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGaa                           >  2:2137830/1‑91 (MQ=255)
                                                                 aaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGaa                           >  2:1217257/1‑91 (MQ=255)
                                                                  aTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAat                          <  1:3592692/91‑1 (MQ=255)
                                                                   tagCGATAGCACATAACTTCATGTAAGTGATATGCTATCGCCATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAatt                         >  1:115836/3‑91 (MQ=14)
                                                                    agCGATAGCACATAACTTCATGTAAGTGATATGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAatta                        >  2:955962/2‑91 (MQ=18)
                                                                     gCGATAGCATATCACTTACATGAAGATATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATGGTACGATGATTTCAGAattat                       >  2:2597272/1‑91 (MQ=255)
                                                                      cGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattatt                      >  1:877093/1‑91 (MQ=255)
                                                                       gATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattattt                     >  2:3858911/1‑90 (MQ=255)
                                                                       gATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattatta                     >  1:1079929/1‑91 (MQ=255)
                                                                        aTAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattattat                    <  2:2198235/91‑1 (MQ=255)
                                                                        aTAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattattat                    >  1:426796/1‑91 (MQ=255)
                                                                           gCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTACGATGATTTAAGAATTATTATTGt                 >  2:3609117/1‑91 (MQ=255)
                                                                            cATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGtt                <  2:3111943/91‑1 (MQ=255)
                                                                             atatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGttt               <  1:1063875/91‑1 (MQ=255)
                                                                             atatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGttt               >  2:1003329/1‑91 (MQ=255)
                                                                              tatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTa              >  1:2156531/1‑91 (MQ=255)
                                                                               atCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTaa             <  1:3155409/91‑1 (MQ=255)
                                                                                 cACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAAtt           >  1:2926543/1‑91 (MQ=255)
                                                                                 cACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAAtt           <  2:2333030/91‑1 (MQ=255)
                                                                                   cTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCa         <  1:3739879/91‑1 (MQ=255)
                                                                                     tACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCaaa       >  1:2471077/1‑91 (MQ=255)
                                                                                     tACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCaaa       <  1:2214299/91‑1 (MQ=255)
                                                                                     tACATGAAGTGATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCaaa       <  2:3550078/91‑1 (MQ=255)
                                                                                        aTGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGt    >  1:2675914/1‑91 (MQ=255)
                                                                                          gAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTAt  >  1:347021/1‑91 (MQ=255)
                                                                                           |                                                                                          
TAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATTCAAATGTAT  >  USA300TCH1516_ALE/631542‑631722

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