Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 631,535 G→A 7.6% intergenic (+161/‑341) 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,5350GA 7.6% 823.2 / 37.4 265intergenic (+161/‑341)gph_1/proPPhosphoglycolate phosphatase/Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base G (110/133);  new base A (12/8);  total (124/141)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 271 total aligned reads displayed.

TTCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAC  >  USA300TCH1516_ALE/631445‑631625
                                                                                          |                                                                                          
ttCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTg                                                                                            <  1:1882467/91‑1 (MQ=255)
 tCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGt                                                                                           >  2:1762407/1‑91 (MQ=255)
 tCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGt                                                                                           <  2:947728/91‑1 (MQ=255)
 tCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGt                                                                                           <  2:2443194/91‑1 (MQ=255)
 tCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGt                                                                                           >  1:1727089/1‑91 (MQ=255)
  cGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTa                                                                                          <  1:2256507/91‑1 (MQ=255)
   gCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTaa                                                                                         >  2:3592692/1‑91 (MQ=255)
   gCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTaa                                                                                         >  2:331601/1‑91 (MQ=255)
     tttAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTaaaa                                                                                       >  2:3664533/1‑91 (MQ=255)
        aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaata                                                                                    >  2:2692140/1‑91 (MQ=255)
        aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaata                                                                                    <  2:2163563/91‑1 (MQ=255)
        aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                    <  1:1332501/91‑1 (MQ=255)
          aTTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaataaa                                                                                  <  1:2921573/91‑1 (MQ=255)
           ttAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaataaat                                                                                 >  1:2900586/1‑91 (MQ=255)
           ttAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGAAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaataaat                                                                                 <  2:338309/91‑1 (MQ=255)
            tAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAAtt                                                                                >  1:2122767/1‑91 (MQ=255)
            tAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAAtt                                                                                <  2:37758/91‑1 (MQ=255)
              ggAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTaa                                                                              <  2:3737042/91‑1 (MQ=255)
               gAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGTTAATTTATTTTTACAAAAATAAATTatc                                                                             <  2:853245/91‑3 (MQ=38)
               gAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTaaa                                                                             <  2:645936/91‑1 (MQ=255)
                  tttAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTa                                                                          >  1:587394/1‑91 (MQ=255)
                  tttAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTa                                                                          >  1:1025543/1‑91 (MQ=255)
                     aaaaTAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATg                                                                       <  2:2359063/91‑1 (MQ=255)
                      aaataaatATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGa                                                                      <  2:398771/91‑1 (MQ=255)
                      aaataaatATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGa                                                                      >  1:2450626/1‑91 (MQ=255)
                       aataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGAc                                                                     <  2:1444114/91‑1 (MQ=255)
                       aataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGAc                                                                     <  2:2244547/91‑1 (MQ=255)
                       aataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGAc                                                                     >  2:519016/1‑91 (MQ=255)
                       aataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGAc                                                                     >  2:3239915/1‑91 (MQ=255)
                       aataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGAc                                                                     >  1:3546635/1‑91 (MQ=255)
                        ataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACa                                                                    <  1:1444607/91‑1 (MQ=255)
                        aataataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACa                                                                    <  2:1384596/88‑1 (MQ=255)
                         taaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACaa                                                                   >  1:721184/1‑91 (MQ=255)
                              atttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTa                                                              >  1:3809016/1‑91 (MQ=255)
                               tttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaata                                                                                    <  1:3167503/68‑1 (MQ=255)
                               tttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTAt                                                             <  1:925354/91‑1 (MQ=255)
                               tttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTCTGTAAAAATAAATTAAAGTAATGACAAAGTTAt                                                             >  1:3181602/1‑91 (MQ=255)
                               tttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTATGTAAAAATAAATTAAAGTAATGACAAAGTTAt                                                             >  1:3181763/1‑91 (MQ=255)
                                ttattaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTAtt                                                            >  1:25746/1‑91 (MQ=255)
                                   ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                         <  1:2230204/91‑1 (MQ=255)
                                   ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                         >  2:3728579/1‑91 (MQ=255)
                                    taAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAg                                                        >  2:100773/1‑91 (MQ=255)
                                    taAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAg                                                        >  1:302777/1‑91 (MQ=255)
                                     aaaCATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGt                                                       >  2:2770670/1‑91 (MQ=255)
                                       aCATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTaa                                                     >  1:31353/1‑91 (MQ=255)
                                       aCATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTaa                                                     >  2:1934090/1‑91 (MQ=255)
                                        cATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTaaa                                                    >  1:3241676/1‑91 (MQ=255)
                                           tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTg                                                 <  2:1067670/91‑1 (MQ=255)
                                              gAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGt                                              >  2:2493410/1‑91 (MQ=255)
                                              gAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGt                                              <  2:455611/91‑1 (MQ=255)
                                               aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                             <  1:2680244/91‑1 (MQ=255)
                                               aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                             <  2:2048214/91‑1 (MQ=255)
                                               aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                             <  1:1180753/91‑1 (MQ=255)
                                               aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                             <  2:302777/91‑1 (MQ=255)
                                               aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                             <  2:3263508/91‑1 (MQ=255)
                                               aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                             <  2:3681445/91‑1 (MQ=255)
                                                aTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGtat                                            >  2:973537/1‑91 (MQ=255)
                                                 tttttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGtata                                           <  1:346809/91‑1 (MQ=255)
                                                   tttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATaaa                                         >  2:220144/1‑91 (MQ=255)
                                                   tttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATaaa                                         <  1:657348/91‑1 (MQ=255)
                                                    ttAAAGAGTAATGTCTGACTCGTTGATAATTTTTTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAAc                                        >  2:846755/1‑91 (MQ=255)
                                                      aaaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAt                                      <  1:3735636/91‑1 (MQ=255)
                                                        agagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAttt                                    >  2:2885870/1‑91 (MQ=255)
                                                           gTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTaaa                                 <  1:1235509/91‑1 (MQ=255)
                                                           gTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTaaa                                 <  1:2274534/91‑1 (MQ=255)
                                                           gTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTaaa                                 <  2:1051881/91‑1 (MQ=255)
                                                             aaTGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATa                               <  2:679011/91‑1 (MQ=255)
                                                             aaTGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATa                               >  1:419693/1‑91 (MQ=255)
                                                              aTGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAGGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATAc                              <  1:3873091/91‑1 (MQ=255)
                                                               tGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACg                             <  1:452512/91‑1 (MQ=255)
                                                                gTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGa                            >  1:904872/1‑91 (MQ=255)
                                                                 tCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGAt                           <  2:231226/91‑1 (MQ=255)
                                                                 tCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGAt                           >  2:1486124/1‑91 (MQ=255)
                                                                 tCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGAt                           >  1:1505733/1‑91 (MQ=255)
                                                                 tCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGATATTGAAGTAAATTGAGTATAAACATTTAAATACGAt                           <  1:3261280/91‑1 (MQ=255)
                                                                   tGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGt                         <  2:3461704/91‑1 (MQ=255)
                                                                   tGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGt                         <  2:1694617/91‑1 (MQ=255)
                                                                     aCTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCg                       <  2:1919909/91‑1 (MQ=255)
                                                                      cTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGa                      <  1:2172236/91‑1 (MQ=255)
                                                                      cTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGa                      <  2:1107630/91‑1 (MQ=255)
                                                                            gATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATgg                >  1:193433/1‑91 (MQ=255)
                                                                             ttAATTTATTTTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGc               >  1:101755/2‑91 (MQ=255)
                                                                              tAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCg              >  1:3671789/1‑91 (MQ=255)
                                                                              tAATTTATTTTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCg              >  1:2854019/1‑91 (MQ=255)
                                                                               aaTTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGa             >  2:2089389/1‑91 (MQ=255)
                                                                                atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt            <  1:3360701/91‑1 (MQ=255)
                                                                                atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt            >  1:2395624/1‑91 (MQ=255)
                                                                                atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt            >  2:1301706/1‑91 (MQ=255)
                                                                                atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt            >  1:63036/1‑91 (MQ=255)
                                                                                 tttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATa           >  2:2070652/1‑91 (MQ=255)
                                                                                  ttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAg          >  1:2191999/1‑91 (MQ=255)
                                                                                  ttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAg          >  2:2582063/1‑91 (MQ=255)
                                                                                   tatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGc         <  2:688797/91‑1 (MQ=255)
                                                                                    atttTTGTAAAAATAAATTAAAGTACTGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCa        <  1:3620942/91‑1 (MQ=255)
                                                                                       tttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCatat     >  2:599307/1‑91 (MQ=255)
                                                                                        ttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATc    <  2:3157474/91‑1 (MQ=255)
                                                                                          gTAAAAATAAATTAAAGTAATGACAAAGTTATTGGAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAc  <  2:1598765/91‑1 (MQ=255)
                                                                                          gTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAc  <  1:1267721/91‑1 (MQ=255)
                                                                                          gTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAc  <  1:1585023/91‑1 (MQ=255)
                                                                                          gTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAc  <  1:3197915/91‑1 (MQ=255)
                                                                                          |                                                                                          
TTCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAC  >  USA300TCH1516_ALE/631445‑631625

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