Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 631,536 T→C 12.4% intergenic (+162/‑340) 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,5360TC12.4% 360.7 / 36.4 129intergenic (+162/‑340)gph_1/proPPhosphoglycolate phosphatase/Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base T (54/59);  new base C (8/8);  total (62/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 132 total aligned reads displayed.

TCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACT  >  USA300TCH1516_ALE/631446‑631626
                                                                                          |                                                                                          
tCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGt                                                                                            >  1:4013779/1‑91 (MQ=255)
tCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGt                                                                                            >  1:586754/1‑91 (MQ=255)
 cGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTa                                                                                           >  1:1905659/1‑91 (MQ=255)
  gCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTaa                                                                                          >  1:2263400/1‑91 (MQ=255)
      tAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaat                                                                                      >  1:906583/1‑91 (MQ=255)
      tAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaat                                                                                      <  1:3411810/91‑1 (MQ=255)
      tAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaat                                                                                      <  1:1737361/91‑1 (MQ=255)
       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaata                                                                                     <  1:3415072/91‑1 (MQ=255)
       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                     <  1:92409/91‑1 (MQ=255)
       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                     <  1:2209248/91‑1 (MQ=255)
       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                     <  1:2209085/91‑1 (MQ=255)
       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                     <  1:2207449/91‑1 (MQ=255)
       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                     <  1:1788762/91‑1 (MQ=255)
        gATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaataa                                                                                    >  1:3472141/1‑91 (MQ=255)
          ttAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaataaat                                                                                  <  1:3738551/91‑1 (MQ=255)
             ggAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTaa                                                                               >  1:306412/1‑91 (MQ=255)
              gAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTaaa                                                                              <  1:184676/91‑1 (MQ=255)
              gAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTaaa                                                                              >  1:66689/1‑91 (MQ=255)
              gAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTaaa                                                                              <  1:2916317/91‑1 (MQ=255)
               aaTTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAg                                                                             <  1:159222/91‑1 (MQ=255)
               aaTTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAg                                                                             <  1:157842/91‑1 (MQ=255)
                 tttAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTa                                                                           <  1:1608232/91‑1 (MQ=255)
                   tAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAAt                                                                         >  1:1730088/1‑91 (MQ=255)
                    aaaaTAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATg                                                                        >  1:2330357/1‑91 (MQ=255)
                       ataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACa                                                                     <  1:415622/91‑1 (MQ=255)
                        taaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACaa                                                                    <  1:3424864/91‑1 (MQ=255)
                         aaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACaaa                                                                   <  1:801952/91‑1 (MQ=255)
                         aaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACaaa                                                                   >  1:3889147/1‑91 (MQ=255)
                         aaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACaaa                                                                   >  1:3889229/1‑91 (MQ=255)
                           atatTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGt                                                                 <  1:1981789/91‑1 (MQ=255)
                             atttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTa                                                               >  1:3402030/1‑91 (MQ=255)
                              tttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTAt                                                              >  1:2331215/1‑91 (MQ=255)
                               ttattaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTAtt                                                             >  1:3152959/1‑91 (MQ=255)
                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                          <  1:2786331/91‑1 (MQ=255)
                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                          >  1:261941/1‑91 (MQ=255)
                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                          <  1:2784944/91‑1 (MQ=255)
                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                          <  1:2472297/91‑1 (MQ=255)
                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                          <  1:2472405/91‑1 (MQ=255)
                                   taAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAg                                                         <  1:39002/91‑1 (MQ=255)
                                   taAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAg                                                         <  1:660157/91‑1 (MQ=255)
                                     aaCATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTa                                                       >  1:746549/1‑91 (MQ=255)
                                      aCATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTaa                                                      >  1:1060081/1‑91 (MQ=255)
                                       cATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTaaa                                                     >  1:318561/1‑91 (MQ=255)
                                        aTTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAAt                                                    <  1:259688/91‑1 (MQ=255)
                                         ttATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAAtt                                                   >  1:2496218/1‑91 (MQ=255)
                                          tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTg                                                  >  1:1259092/1‑91 (MQ=255)
                                          tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTg                                                  >  1:741162/1‑91 (MQ=255)
                                          tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTg                                                  >  1:741000/1‑91 (MQ=255)
                                          tATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTg                                                  <  1:777472/91‑1 (MQ=255)
                                           aTGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGa                                                 >  1:2370914/1‑91 (MQ=255)
                                             gAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGt                                               <  1:1186433/91‑1 (MQ=255)
                                              aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                              <  1:1793941/91‑1 (MQ=255)
                                                 ttttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATaa                                           <  1:2684820/91‑1 (MQ=255)
                                                 ttttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATaa                                           <  1:540708/91‑1 (MQ=255)
                                                 ttttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATaa                                           <  1:597333/91‑1 (MQ=255)
                                                   ttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAAc                                         >  1:1480258/1‑91 (MQ=255)
                                                   ttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAAc                                         >  1:1480765/1‑91 (MQ=255)
                                                     aaaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAt                                       >  1:3406382/1‑91 (MQ=255)
                                                     aaaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAt                                       >  1:108904/1‑91 (MQ=255)
                                                     aaaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAt                                       >  1:109046/1‑91 (MQ=255)
                                                      aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAtt                                      >  1:2102003/1‑91 (MQ=255)
                                                      aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAtt                                      >  1:2035442/1‑91 (MQ=255)
                                                      aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAtt                                      >  1:767878/1‑91 (MQ=255)
                                                      aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAtt                                      >  1:1925409/1‑91 (MQ=255)
                                                      aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAtt                                      >  1:3559666/1‑91 (MQ=255)
                                                       agagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAttt                                     <  1:4077953/91‑1 (MQ=255)
                                                       agagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAttt                                     <  1:2384454/91‑1 (MQ=255)
                                                       agagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAttt                                     >  1:1722043/1‑91 (MQ=255)
                                                       agagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAttt                                     >  1:2543151/1‑91 (MQ=255)
                                                        gagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTa                                    <  1:1817571/91‑1 (MQ=255)
                                                        gagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTa                                    <  1:2089829/91‑1 (MQ=255)
                                                         agTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTaa                                   >  1:2675780/1‑91 (MQ=255)
                                                         agTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTaa                                   <  1:474323/91‑1 (MQ=255)
                                                          gTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTaaa                                  <  1:2067123/91‑1 (MQ=255)
                                                          gTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTaaa                                  <  1:2769930/91‑1 (MQ=255)
                                                            aaTGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATa                                <  1:1146015/91‑1 (MQ=255)
                                                             aTGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATAc                               <  1:3623591/91‑1 (MQ=255)
                                                               gTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGa                             <  1:555056/91‑1 (MQ=255)
                                                                tCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGAt                            >  1:2325575/1‑91 (MQ=255)
                                                                 cTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATg                           <  1:2484010/91‑1 (MQ=255)
                                                                        gTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGaaaa                    <  1:3079059/91‑1 (MQ=255)
                                                                        gTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGaaaa                    >  1:2422992/1‑91 (MQ=255)
                                                                           tttAATTTATTTTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATgg                 >  1:2655737/3‑91 (MQ=38)
                                                                           gATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATgg                 >  1:3728475/1‑91 (MQ=255)
                                                                           gATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATgg                 >  1:3727758/1‑91 (MQ=255)
                                                                             tAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCg               <  1:3681964/91‑1 (MQ=255)
                                                                             tAATTTATTTTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCg               >  1:277367/1‑91 (MQ=255)
                                                                              aaTTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGa              <  1:1269521/91‑1 (MQ=255)
                                                                               atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt             <  1:2495756/91‑1 (MQ=255)
                                                                               atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt             <  1:1893041/91‑1 (MQ=255)
                                                                                tttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATa            >  1:4134213/1‑91 (MQ=255)
                                                                                tttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATa            >  1:1945249/1‑91 (MQ=255)
                                                                                  tatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGc          >  1:3873696/1‑91 (MQ=255)
                                                                                  tatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGc          <  1:4000302/91‑1 (MQ=255)
                                                                                     ttttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCata       <  1:2451312/91‑1 (MQ=255)
                                                                                      tttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCatat      <  1:1439805/91‑1 (MQ=255)
                                                                                       ttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATc     >  1:883247/1‑91 (MQ=255)
                                                                                        tGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCa    >  1:2244604/1‑91 (MQ=255)
                                                                                        tGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCa    >  1:2244485/1‑91 (MQ=255)
                                                                                          tAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACt  >  1:3020567/1‑91 (MQ=255)
                                                                                          |                                                                                          
TCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACT  >  USA300TCH1516_ALE/631446‑631626

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