Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 631,536 T→C 9.2% 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,5360TC 9.2% 778.0 / 51.6 269intergenic (+162/‑340)gph_1/proPPhosphoglycolate phosphatase/Proline/betaine transporter
Reads supporting (aligned to +/- strand):  ref base T (108/135);  new base C (10/15);  total (118/151)
Fisher's exact test for biased strand distribution p-value = 8.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 276 total aligned reads displayed.

TCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACT  >  USA300TCH1516_ALE/631446‑631626
                                                                                          |                                                                                          
tCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGt                                                                                            >  2:2900954/1‑91 (MQ=255)
 cGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTa                                                                                           >  1:2228634/1‑91 (MQ=255)
 cGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTa                                                                                           <  2:1238242/91‑1 (MQ=255)
  gCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTaa                                                                                          >  1:1263093/1‑91 (MQ=255)
     ttAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTaaaaa                                                                                       <  1:1333192/91‑1 (MQ=255)
      tAGATTAGGAATTTAAAATAAATATGTATTAGACATTATGAATTTTTAAAGAGTAATGTCTTACTCGTTGATAATTTATTTTTACAAaaat                                                                                      <  2:2467641/91‑1 (MQ=255)
       aGATTCGGGATTTAAAATGAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                     <  2:3654939/91‑1 (MQ=255)
       aGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaata                                                                                     <  2:3071668/91‑1 (MQ=255)
        gtttAGGAATTTAAAATATATATTTATTAAACATTATGGATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaataa                                                                                    <  2:3878549/89‑1 (MQ=255)
        gATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaataa                                                                                    >  1:2009971/1‑91 (MQ=255)
         aTTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTACAAaaataaa                                                                                   <  2:2788333/91‑1 (MQ=255)
          ttAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTGGTTGATAATTTATTTTTGTAAaaataaat                                                                                  <  2:3438105/91‑1 (MQ=255)
          ttAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaataaat                                                                                  >  2:2667786/1‑91 (MQ=255)
          ttAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaataaat                                                                                  <  1:2218273/91‑1 (MQ=255)
          ttAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAaaataaat                                                                                  <  2:454703/91‑1 (MQ=255)
           tAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAtaaataaaat                                                                                 >  2:1636620/1‑89 (MQ=255)
           tAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAAtt                                                                                 >  1:3751257/1‑91 (MQ=255)
           tAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAAtt                                                                                 >  2:299626/1‑91 (MQ=255)
           tAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAAtt                                                                                 >  1:2942713/1‑91 (MQ=255)
             ggAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTaa                                                                               >  1:1617438/1‑91 (MQ=255)
             ggAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTaa                                                                               >  2:380831/1‑91 (MQ=255)
               aaTTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAg                                                                             <  1:1869271/91‑1 (MQ=255)
                aTTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGt                                                                            <  2:3782916/91‑1 (MQ=255)
                aTTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGt                                                                            >  1:3658327/1‑91 (MQ=255)
                 tttAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTa                                                                           <  2:2095540/91‑1 (MQ=255)
                 tttAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTa                                                                           <  2:1544956/91‑1 (MQ=255)
                  ttAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTaa                                                                          >  2:338074/1‑91 (MQ=255)
                    aaaaTAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATg                                                                        >  1:2578170/1‑91 (MQ=255)
                    aaaaTAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATg                                                                        <  1:3062451/91‑1 (MQ=255)
                      aataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGAc                                                                      <  2:2517175/91‑1 (MQ=255)
                      aataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATGTTTGTAAAAATAAATTAAAGTAATGAc                                                                      <  1:2131124/91‑1 (MQ=255)
                       ataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACa                                                                     >  2:2069304/1‑91 (MQ=255)
                       ataaataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACa                                                                     >  1:482231/1‑91 (MQ=255)
                          aataTTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAg                                                                  <  2:3018839/91‑1 (MQ=255)
                           atatTTATTTAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGt                                                                 >  1:444613/1‑91 (MQ=255)
                           atatTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGt                                                                 <  1:3505906/91‑1 (MQ=255)
                           atatTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGt                                                                 <  1:1796224/91‑1 (MQ=255)
                             atttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTa                                                               >  2:2172672/1‑91 (MQ=255)
                              tttattAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTAt                                                              >  1:1208824/1‑91 (MQ=255)
                               ttattaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTAtt                                                             <  1:824157/91‑1 (MQ=255)
                                tattaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTg                                                            >  2:1758008/1‑91 (MQ=255)
                                tattaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTg                                                            <  2:236744/91‑1 (MQ=255)
                                 aataAATATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGa                                                           >  2:2268700/3‑91 (MQ=255)
                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                          <  2:956130/91‑1 (MQ=255)
                                  ttaAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGaa                                                          <  1:2916857/91‑1 (MQ=255)
                                   taAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAg                                                         >  1:2963853/1‑91 (MQ=255)
                                    aaaCATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGt                                                        <  1:928163/91‑1 (MQ=255)
                                    aaaCATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGt                                                        <  2:2391375/91‑1 (MQ=255)
                                    aaaCATTATGAATTTTGAAAGAGTAATGTCTGACTCGTTGATAATTTATTTATGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGt                                                        <  2:3776655/91‑1 (MQ=255)
                                     aaCATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTa                                                       <  1:552633/91‑1 (MQ=255)
                                         ttATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAAtt                                                   <  2:702542/91‑1 (MQ=255)
                                              aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                              <  2:1829112/91‑1 (MQ=255)
                                              aaTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                              >  2:1263437/1‑91 (MQ=255)
                                              aaTTTTTAAAGAGGAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGta                                              <  1:1697024/91‑1 (MQ=255)
                                               aTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTTTTTTTTTTATTTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGtat                                             >  2:1931499/1‑91 (MQ=255)
                                               aTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGtat                                             >  1:219373/1‑91 (MQ=255)
                                               aTTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGtat                                             <  1:1589426/91‑1 (MQ=255)
                                                  tttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATaaa                                          <  2:2029183/91‑1 (MQ=255)
                                                  tttAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATaaa                                          >  1:3872620/1‑91 (MQ=255)
                                                    tAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACa                                        >  2:3097959/1‑91 (MQ=255)
                                                    tAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACa                                        <  1:1456988/91‑1 (MQ=255)
                                                     aaaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAt                                       >  2:671013/1‑91 (MQ=255)
                                                      aaGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAtt                                      >  1:3161368/1‑91 (MQ=255)
                                                       agagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACAttt                                     <  1:3322373/91‑1 (MQ=255)
                                                        gagTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTa                                    >  2:1319071/1‑91 (MQ=255)
                                                           tAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAAt                                 <  2:1624290/91‑1 (MQ=255)
                                                           tAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAAt                                 <  1:978820/91‑1 (MQ=255)
                                                            aaTTTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATa                                <  2:3737034/91‑1 (MQ=255)
                                                             aTGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATAc                               >  1:2620931/1‑91 (MQ=255)
                                                              tgccTGACTCGTTGATAATTTATTTTTGTGAAAATACATTAAAGTAATGACAACGTCATTGCAGTAAATTGAGTGCAGACATTTCAATCCg                              >  2:3529630/4‑91 (MQ=255)
                                                                tCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGAt                            <  1:2301767/91‑1 (MQ=255)
                                                                  tGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGt                          <  1:1649217/91‑1 (MQ=255)
                                                                   gACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTc                         <  2:822472/91‑1 (MQ=255)
                                                                    aCTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCg                        >  1:415967/1‑91 (MQ=255)
                                                                      tCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGaa                      >  1:642461/1‑91 (MQ=255)
                                                                        gTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGaaaa                    <  2:2004242/91‑1 (MQ=255)
                                                                         ttGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAAt                   >  2:2584954/1‑91 (MQ=255)
                                                                         ttGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAAt                   <  2:2973177/91‑1 (MQ=255)
                                                                         ttGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAAt                   <  2:703310/91‑1 (MQ=255)
                                                                         ttGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAAt                   >  1:2438180/1‑91 (MQ=255)
                                                                          tGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATg                  >  1:2250409/1‑91 (MQ=255)
                                                                             tAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCg               <  1:370034/91‑1 (MQ=255)
                                                                              aaTTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGa              >  1:717146/1‑91 (MQ=255)
                                                                               atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt             <  1:1110067/91‑1 (MQ=255)
                                                                               atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt             <  2:1076973/91‑1 (MQ=255)
                                                                               atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt             <  1:2832944/91‑1 (MQ=255)
                                                                               atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAAGTGGGTATAAACATTTAAATACGATGTCGAAAATGGCGAt             <  2:1418401/91‑1 (MQ=255)
                                                                               atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGTTAGTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt             <  2:3478383/91‑1 (MQ=255)
                                                                               atttatttTTGTAAAAATAAATTAAAGTAATGACAAAGCTATTGAAGTAAATTGAGAATAAACATTTAAATACGATGTCGAAAATGGCGaa             >  2:3291712/1‑90 (MQ=255)
                                                                               atttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGAt             >  1:3361368/1‑91 (MQ=255)
                                                                                 ttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAg           >  2:1299436/1‑91 (MQ=255)
                                                                                 ttatttTTACAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAg           >  1:3306775/1‑91 (MQ=255)
                                                                                   atttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCa         <  2:1155541/91‑1 (MQ=255)
                                                                                   atttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCa         <  1:1927134/91‑1 (MQ=255)
                                                                                   atttTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCa         >  1:1105265/1‑91 (MQ=255)
                                                                                    tttttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCat        <  2:586463/91‑1 (MQ=255)
                                                                                     ttttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGGAAATGGCGATAGCata       >  2:1698181/1‑91 (MQ=255)
                                                                                      tttGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCatat      >  2:3339558/1‑91 (MQ=255)
                                                                                         gTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCAc   <  2:857919/91‑1 (MQ=255)
                                                                                          tAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACt  <  2:482231/91‑1 (MQ=255)
                                                                                          |                                                                                          
TCGCTTTAGATTAGGAATTTAAAATAAATATTTATTAAACATTATGAATTTTTAAAGAGTAATGTCTGACTCGTTGATAATTTATTTTTGTAAAAATAAATTAAAGTAATGACAAAGTTATTGAAGTAAATTGAGTATAAACATTTAAATACGATGTCGAAAATGGCGATAGCATATCACT  >  USA300TCH1516_ALE/631446‑631626

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