Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 631,620 T→C 5.8% intergenic (+246/‑256) 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,6200TC 5.8% 675.9 / 20.7 207intergenic (+246/‑256)gph_1/proPPhosphoglycolate phosphatase/Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base T (111/84);  new base C (5/7);  total (116/91)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 220 total aligned reads displayed.

TTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTA  >  USA300TCH1516_ALE/631530‑631710
                                                                                          |                                                                                           
tttttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCat                                                                                             >  2:3409686/1‑91 (MQ=255)
 ttttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCata                                                                                            >  1:652689/1‑91 (MQ=255)
 ttttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCata                                                                                            <  1:1486049/91‑1 (MQ=255)
 ttttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCata                                                                                            <  2:3859463/91‑1 (MQ=255)
  tttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCatat                                                                                           >  2:599307/1‑91 (MQ=255)
    tGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATATGATGTCGAAAATGGCGATAGCATATCa                                                                                         >  1:389632/1‑91 (MQ=255)
     gTAAAAATAAATTAAAGTAATGACAAAGTTATTGGAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAc                                                                                        <  2:1598765/91‑1 (MQ=255)
     gTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAc                                                                                        <  1:3197915/91‑1 (MQ=255)
     gTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAc                                                                                        <  1:1267721/91‑1 (MQ=255)
      tAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACt                                                                                       >  2:2021188/1‑91 (MQ=255)
       aaaaaTAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACtt                                                                                      <  2:3663114/91‑1 (MQ=255)
        aaaaTAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTa                                                                                     >  1:926938/1‑91 (MQ=255)
          aataaatTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACa                                                                                   <  2:424316/91‑1 (MQ=255)
          aataaatTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACa                                                                                   >  1:2333030/1‑91 (MQ=255)
             aaatTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTCAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGa                                                                               <  1:3359652/91‑1 (MQ=255)
              aatTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGaa                                                                              >  2:2680621/1‑91 (MQ=255)
               atTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAg                                                                             >  1:940654/1‑91 (MQ=255)
                ttAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGt                                                                            <  1:331601/91‑1 (MQ=255)
                 tAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGtt                                                                           >  2:3448687/1‑91 (MQ=255)
                  aaaGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTa                                                                          >  2:2134657/1‑91 (MQ=255)
                   aaGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTAt                                                                         >  1:1516615/1‑91 (MQ=255)
                     gTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTAtgt                                                                       <  2:25746/91‑1 (MQ=255)
                         tGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTa                                                                   <  2:2683913/91‑1 (MQ=255)
                          gACAAAGTTATTGAAGTAACTTGAGTATATACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTAt                                                                  >  2:3386998/1‑91 (MQ=255)
                             aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGc                                                               <  1:2420474/91‑1 (MQ=255)
                             aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCACATAACTTCATGTAAGTGATATGCTATCGc                                                               <  1:68279/91‑1 (MQ=18)
                             aaaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATAGCGATAGCACATAACTTCATGTAAGTGATATGCTATCGc                                                               <  1:2179083/91‑1 (MQ=17)
                              aaGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCt                                                              >  2:3526667/1‑91 (MQ=255)
                               aGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTa                                                             >  1:3373011/1‑91 (MQ=255)
                               aGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTa                                                             >  1:800495/1‑91 (MQ=255)
                                gTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAt                                                            <  1:1183277/91‑1 (MQ=255)
                                 ttATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAtt                                                           >  2:1947147/1‑91 (MQ=255)
                                 ttATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAtt                                                           <  1:1590582/91‑1 (MQ=255)
                                   aTTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAtttt                                                         <  1:519016/91‑1 (MQ=255)
                                    ttGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTAttttt                                                        >  2:2405134/1‑91 (MQ=255)
                                     tGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTa                                                       >  1:18911/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)
                                             aTTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAAt                                               >  2:757511/1‑91 (MQ=255)
                                              ttGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAAtt                                              <  1:464009/91‑1 (MQ=255)
                                              ttGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAAtt                                              >  2:1357459/1‑91 (MQ=255)
                                                gAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTcc                                            >  1:2560376/1‑91 (MQ=255)
                                                gAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTcc                                            >  1:2072892/1‑91 (MQ=255)
                                                  gTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaa                                          >  2:3165107/1‑91 (MQ=255)
                                                  gTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaa                                          >  2:3322733/1‑91 (MQ=255)
                                                   tataAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaaa                                         >  2:3061012/1‑91 (MQ=255)
                                                   tataAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaaa                                         >  1:1788528/1‑91 (MQ=255)
                                                    ataAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaaaa                                        >  2:2118528/1‑91 (MQ=255)
                                                    ataAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATCTTTAGTTATAATTCCaaaa                                        >  2:1444766/1‑91 (MQ=255)
                                                     taAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCaaaaa                                       >  1:2705325/1‑91 (MQ=255)
                                                      aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                      <  2:63036/91‑1 (MQ=255)
                                                      aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                      >  1:285471/1‑91 (MQ=255)
                                                      aaaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAg                                      >  2:2177586/1‑91 (MQ=255)
                                                       aaCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGt                                     <  2:3546635/91‑1 (MQ=255)
                                                        aCATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGtt                                    <  2:2073265/91‑1 (MQ=255)
                                                         cATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTa                                   >  2:3155409/1‑91 (MQ=255)
                                                             tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                               <  2:2512522/91‑1 (MQ=255)
                                                             tAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCg                               <  1:3676169/91‑1 (MQ=255)
                                                              aaaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGt                              >  1:3593548/1‑91 (MQ=255)
                                                               aaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGtt                             >  2:1120234/1‑91 (MQ=255)
                                                               aaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGtt                             <  1:220144/91‑1 (MQ=255)
                                                               aaTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGtt                             <  1:2374831/91‑1 (MQ=255)
                                                                aTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTc                            <  2:960433/91‑1 (MQ=255)
                                                                aTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTc                            >  1:2688532/1‑91 (MQ=255)
                                                                aTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTc                            >  2:1617397/1‑91 (MQ=255)
                                                                aTACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTc                            >  1:1648801/1‑91 (MQ=255)
                                                                 tACGATGTCGAAAATGGCGATAGCATATCACTTTCATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCg                           >  2:2831215/1‑91 (MQ=255)
                                                                 tACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCg                           <  1:2093569/91‑1 (MQ=255)
                                                                 tACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCg                           <  1:608165/91‑1 (MQ=255)
                                                                  aCGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCga                          >  2:943158/1‑91 (MQ=255)
                                                                     aTGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCgatga                       <  1:1693055/91‑1 (MQ=255)
                                                                     aTGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCgatga                       >  1:3111536/1‑91 (MQ=255)
                                                                       gTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGAtt                     <  2:1121479/91‑1 (MQ=255)
                                                                       gTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGAtt                     >  2:2061046/1‑91 (MQ=255)
                                                                        tCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGAttt                    <  1:722456/91‑1 (MQ=255)
                                                                         cGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTa                   >  1:2389846/1‑91 (MQ=255)
                                                                         cGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTa                   >  1:3599115/1‑91 (MQ=255)
                                                                          gAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTaa                  >  2:2206512/1‑91 (MQ=255)
                                                                          gAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTaa                  <  1:3526667/91‑1 (MQ=255)
                                                                            aaaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGa                >  2:2290256/1‑91 (MQ=255)
                                                                             aaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGaa               >  1:3788377/1‑91 (MQ=255)
                                                                             aaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGaa               >  2:1217257/1‑91 (MQ=255)
                                                                             aaTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGaa               <  1:972272/91‑1 (MQ=255)
                                                                              aTGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAat              <  1:3592692/91‑1 (MQ=255)
                                                                               tagCGATAGCACATAACTTCATGTAAGTGATATGCTATCGCCATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAatt             >  1:115836/3‑91 (MQ=14)
                                                                               tGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAatt             <  1:1335375/91‑1 (MQ=255)
                                                                                  cGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattatt          >  1:877093/1‑91 (MQ=255)
                                                                                  cGATAGCACATAACTTCATGTAAGTGATATGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattatt          >  2:1714087/1‑91 (MQ=18)
                                                                                   gATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattatta         >  2:3858855/1‑91 (MQ=255)
                                                                                   gATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattatta         >  1:1079929/1‑91 (MQ=255)
                                                                                    aTAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattattat        >  1:426796/1‑91 (MQ=255)
                                                                                    aTAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAattattat        >  2:1644816/1‑91 (MQ=255)
                                                                                       gCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGt     >  2:3742037/1‑91 (MQ=255)
                                                                                         atatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGttt   <  1:1472522/91‑1 (MQ=255)
                                                                                         atatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGttt   >  2:1003329/1‑91 (MQ=255)
                                                                                         atatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGttt   <  2:1909106/91‑1 (MQ=255)
                                                                                         atatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGttt   <  1:1063875/91‑1 (MQ=255)
                                                                                          tatCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTa  >  1:2156531/1‑91 (MQ=255)
                                                                                          |                                                                                           
TTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACTTACATGAAGTTATGTGCTATCGCTATTTTTAGTTATAATTCCAAAAAGTTAATCGTTCGATGATTTAAGAATTATTATTGTTTA  >  USA300TCH1516_ALE/631530‑631710

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