Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 631,623 C→A 5.2% intergenic (+249/‑253) 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,6230CA 5.2% 747.3 / 23.4 214intergenic (+249/‑253)gph_1/proPPhosphoglycolate phosphatase/Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base C (114/88);  new base A (4/7);  total (119/95)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 217 total aligned reads displayed.

TTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATT  >  USA300TCH1516_ALE/631533‑631713
                                                                                          |                                                                                           
ttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATc                                                                                             <  2:3157474/91‑1 (MQ=255)
 tGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATATGATGTCGAAAATGGCGATAGCATATCa                                                                                            >  1:389632/1‑91 (MQ=255)
  gTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAc                                                                                           <  1:3197915/91‑1 (MQ=255)
  gTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAc                                                                                           <  1:1267721/91‑1 (MQ=255)
     aaaaTAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTa                                                                                        <  2:3681514/91‑1 (MQ=255)
       aataaatTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACa                                                                                      <  2:424669/91‑1 (MQ=255)
       aataaatTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACa                                                                                      >  1:2333030/1‑91 (MQ=255)
          aaatTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTCAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGa                                                                                  <  1:3359652/91‑1 (MQ=255)
            atTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAg                                                                                >  1:940591/1‑91 (MQ=255)
             ttAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGt                                                                               <  1:331601/91‑1 (MQ=255)
              tAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGtt                                                                              <  1:1917237/91‑1 (MQ=255)
               aaaGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTa                                                                             >  2:2134657/1‑91 (MQ=255)
                      tGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGACAATGGCGATAGCATATGAATTACATGGAGTTAAGTGCTa                                                                      >  2:2568334/1‑91 (MQ=255)
                      tGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTa                                                                      <  2:2683832/91‑1 (MQ=255)
                       gACAAAGTTATTGAAGTAACTTGAGTATATACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTAt                                                                     >  2:3386998/1‑91 (MQ=255)
                       gACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCACATAACTTCATGTAAGTGATATGCTAt                                                                     <  1:539799/91‑1 (MQ=18)
                         cAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCg                                                                   <  2:1558905/91‑1 (MQ=255)
                          aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGc                                                                  <  1:2420474/91‑1 (MQ=255)
                          aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATAGCGATAGCACATAACTTCATGTAAGTGATATGCTATCGc                                                                  <  1:492416/91‑1 (MQ=17)
                          aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATAGCGATAGCACATAACTTCATGTAAGTGATATGCTATCGc                                                                  <  2:2750245/91‑1 (MQ=17)
                           aaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCt                                                                 >  2:3526667/1‑91 (MQ=255)
                           aaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCt                                                                 >  2:2366172/1‑91 (MQ=255)
                           aaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCt                                                                 <  2:852176/91‑1 (MQ=255)
                            aGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTa                                                                >  1:800495/1‑91 (MQ=255)
                            aGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTa                                                                >  1:3373011/1‑91 (MQ=255)
                             gTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAt                                                               >  1:1886403/1‑91 (MQ=255)
                             gTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAt                                                               <  1:1183277/91‑1 (MQ=255)
                              ttATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAtt                                                              >  2:1944332/1‑91 (MQ=255)
                              ttATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAtt                                                              <  1:1590582/91‑1 (MQ=255)
                                 ttGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAttttt                                                           <  1:2199720/91‑1 (MQ=255)
                                 ttGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAttttt                                                           >  1:1720403/1‑91 (MQ=255)
                                 ttGAAGTAAATTGAGTATAAACATTTAAATACCATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAttttt                                                           >  2:3808757/1‑91 (MQ=255)
                                  tGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTa                                                          >  1:22952/1‑91 (MQ=255)
                                    aaGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGt                                                        >  1:1404749/1‑91 (MQ=255)
                                     aGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGtt                                                       >  1:473353/1‑91 (MQ=255)
                                        aaaTTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTtata                                                    <  1:2908465/91‑1 (MQ=255)
                                          aTTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAAt                                                  <  1:2885870/91‑1 (MQ=255)
                                           ttGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAAtt                                                 >  2:1357459/1‑91 (MQ=255)
                                           ttGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAAtt                                                 <  1:464009/91‑1 (MQ=255)
                                             gAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTcc                                               >  2:3410314/1‑91 (MQ=255)
                                             gAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTcc                                               >  1:2072892/1‑91 (MQ=255)
                                             gAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTCCTATCGCTATCTTTAGTTATAATTcc                                               >  2:3360193/1‑91 (MQ=255)
                                                 ataAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaaaa                                           >  1:719857/1‑91 (MQ=255)
                                                 ataAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaaaa                                           >  2:2118528/1‑91 (MQ=255)
                                                  taAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaaaaa                                          >  1:2705325/1‑91 (MQ=255)
                                                   aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                         >  2:194837/1‑91 (MQ=255)
                                                   aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                         >  2:2177586/1‑91 (MQ=255)
                                                   aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                         <  1:2985805/91‑1 (MQ=255)
                                                   aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                         <  1:757511/91‑1 (MQ=255)
                                                   aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                         >  1:285471/1‑91 (MQ=255)
                                                    aaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGt                                        >  1:1520389/1‑91 (MQ=255)
                                                    aaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGt                                        <  2:3546635/91‑1 (MQ=255)
                                                      cATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTa                                      >  2:3240267/1‑91 (MQ=255)
                                                       aTTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTaa                                     <  1:1486124/91‑1 (MQ=255)
                                                          tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                                  <  1:971938/91‑1 (MQ=255)
                                                          tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                                  <  1:973537/91‑1 (MQ=255)
                                                          tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                                  <  2:607266/91‑1 (MQ=255)
                                                          tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                                  <  1:3676169/91‑1 (MQ=255)
                                                           aaaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGt                                 <  1:3728579/91‑1 (MQ=255)
                                                           aaaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGt                                 >  1:3593548/1‑91 (MQ=255)
                                                            aaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGtt                                <  1:2374831/91‑1 (MQ=255)
                                                            aaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGtt                                >  2:1825566/1‑91 (MQ=255)
                                                            aaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGtt                                >  2:1120234/1‑91 (MQ=255)
                                                             aTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTc                               >  1:1648801/1‑91 (MQ=255)
                                                             aTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTc                               >  1:2570616/1‑91 (MQ=255)
                                                             aTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTc                               >  1:2688532/1‑91 (MQ=255)
                                                              tACGATGTCGAAAATGGCGATAGCATATCACTTTCATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCg                              >  2:2831215/1‑91 (MQ=255)
                                                              tACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCg                              <  1:2093569/91‑1 (MQ=255)
                                                                cGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCgat                            >  2:2452380/1‑91 (MQ=255)
                                                                    gTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGAtt                        >  2:2060999/1‑91 (MQ=255)
                                                                    gTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGAtt                        >  2:2062465/1‑91 (MQ=255)
                                                                    gTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGAtt                        <  2:1121479/91‑1 (MQ=255)
                                                                     tCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGAttt                       >  1:1317241/1‑91 (MQ=255)
                                                                      cGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTa                      >  1:3599115/1‑91 (MQ=255)
                                                                      cGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTa                      >  1:2389846/1‑91 (MQ=255)
                                                                       gAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTaa                     >  2:2206512/1‑91 (MQ=255)
                                                                       gAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTaa                     <  1:3526667/91‑1 (MQ=255)
                                                                         aaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                   <  2:162735/91‑1 (MQ=255)
                                                                         aaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                   <  2:160288/91‑1 (MQ=255)
                                                                         aaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                   >  2:2290256/1‑91 (MQ=255)
                                                                         aaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                   <  2:705628/91‑1 (MQ=255)
                                                                          aaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGaa                  >  2:1217257/1‑91 (MQ=255)
                                                                           aTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAat                 >  2:1435823/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)
                                                                               cGATAGCACATAACTTCATGTAAGTGATATGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattatt             >  2:1714087/1‑91 (MQ=18)
                                                                                gATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattattt            >  2:3858911/1‑90 (MQ=255)
                                                                                gATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattatta            >  1:655138/1‑91 (MQ=255)
                                                                                gATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattatta            >  1:1079929/1‑91 (MQ=255)
                                                                                 aTAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattattat           >  1:426796/1‑91 (MQ=255)
                                                                                    gCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTACGATGATTTAAGAATTATTATTGt        >  2:3609117/1‑91 (MQ=255)
                                                                                     cATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGtt       <  2:3111943/91‑1 (MQ=255)
                                                                                      atatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGttt      >  2:1003329/1‑91 (MQ=255)
                                                                                      atatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGttt      <  1:1472522/91‑1 (MQ=255)
                                                                                      atatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGttt      <  1:1063875/91‑1 (MQ=255)
                                                                                       tatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTa     >  1:2156531/1‑91 (MQ=255)
                                                                                        atCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTaa    <  1:3155409/91‑1 (MQ=255)
                                                                                          |                                                                                           
TTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTAATT  >  USA300TCH1516_ALE/631533‑631713

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: