Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 253,983 A→G 9.5% intergenic (+196/+166) asbF → / ← USA300HOU_RS01125 3‑dehydroshikimate dehydratase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE253,9830AG 9.5% 717.6 / 45.5 243intergenic (+196/+166)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (122/96);  new base G (8/15);  total (132/111)
Fisher's exact test for biased strand distribution p-value = 7.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 253 total aligned reads displayed.

CAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTC  >  USA300TCH1516_ALE/253893‑254073
                                                                                          |                                                                                          
 aaGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAagag                                                                                           <  2:486464/91‑1 (MQ=255)
 aaGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAagag                                                                                           >  1:1293872/1‑91 (MQ=255)
  tgTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATTTAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                          >  1:3580910/2‑91 (MQ=255)
  aGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                          >  2:903740/1‑91 (MQ=255)
  aGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                          >  1:1242111/1‑91 (MQ=255)
   gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGta                                                                                         >  1:2220240/1‑91 (MQ=255)
   gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGta                                                                                         <  1:2793269/91‑1 (MQ=255)
     cATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtata                                                                                       <  1:1901499/91‑1 (MQ=255)
      aTTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatat                                                                                      >  1:968139/1‑91 (MQ=255)
      aTTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatat                                                                                      >  1:802258/1‑91 (MQ=255)
      aTTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatat                                                                                      <  1:2736964/91‑1 (MQ=255)
       ttgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatata                                                                                     >  1:1312388/1‑91 (MQ=255)
       ttgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatata                                                                                     <  2:308952/91‑1 (MQ=255)
         gtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAc                                                                                   >  1:560900/1‑91 (MQ=255)
          tttgtAAAGAAGGTGTACTTTATGTAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                  <  2:2839822/91‑1 (MQ=255)
          tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGGTtatatactc                                                                                  <  1:1208220/91‑4 (MQ=255)
          tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGGTtatatactc                                                                                  <  1:1854553/91‑4 (MQ=255)
          tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                  <  2:313742/91‑1 (MQ=255)
          tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                  <  2:2197397/91‑1 (MQ=255)
          tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                  <  2:3836560/91‑1 (MQ=255)
          tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                  <  1:1129276/91‑1 (MQ=255)
            tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                <  2:1091178/91‑1 (MQ=255)
            tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                <  2:3665497/91‑1 (MQ=255)
             gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCttt                                                                               <  1:3835372/91‑1 (MQ=255)
              tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtttt                                                                              <  1:3106763/91‑1 (MQ=255)
               aaaGGAGGTGTACTTTATATAACTATATAGCGATGGTCATACCCATCCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                             >  2:2048330/1‑91 (MQ=255)
               aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                             <  1:2293851/91‑1 (MQ=255)
               aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                             >  2:1667829/1‑91 (MQ=255)
               aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                             >  2:1731353/1‑91 (MQ=255)
               aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                             >  1:3154564/1‑91 (MQ=255)
               aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                             <  2:3205917/91‑1 (MQ=255)
               aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                             <  1:3493129/91‑1 (MQ=255)
                aagaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCa                                                                            <  2:551134/91‑1 (MQ=255)
                 agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCaa                                                                           <  1:3499799/91‑1 (MQ=255)
                   aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAtg                                                                         <  1:1464236/91‑2 (MQ=25)
                    agGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAtgg                                                                        <  2:3920007/91‑3 (MQ=21)
                    agGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                                        <  1:3815048/91‑1 (MQ=255)
                     ggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGt                                                                       >  1:692346/1‑91 (MQ=255)
                        gtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAg                                                                    <  2:3860964/91‑1 (MQ=255)
                           cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                 <  2:2958649/91‑1 (MQ=255)
                              tatataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTatat                                                              <  2:3309184/91‑1 (MQ=255)
                               atctaAGTATATAGCGATGGTCATACCCATCCACAGTACCAATCCTCACCATTGCAAACCGCCTATAACCTTTTCAACCGTGCGGTATATg                                                             >  2:442158/4‑91 (MQ=255)
                               atataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATg                                                             >  2:1258507/1‑91 (MQ=255)
                               atataAGTATATAGCGATAGTCATACCCATTCACAGTAACAATACTCACCATTGAAAAGAGGATATAACCTTTGCAATAGTGAGCTATATg                                                             >  2:2648114/1‑91 (MQ=255)
                                   aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                         <  2:2769104/91‑1 (MQ=255)
                                    aGTATATAGCGATGGTCATACCCATTCACAGTACCAATCCTCACCATTGCAAAGAGTACATAACCTTTTCAATAGTGAGGTATATGataat                                                        >  2:144389/1‑91 (MQ=255)
                                    aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                        >  2:3067340/1‑91 (MQ=255)
                                    aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                        >  1:3845178/1‑91 (MQ=255)
                                    aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTAAATGataat                                                        >  2:3696684/1‑91 (MQ=255)
                                      tatataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAtaataa                                                      <  1:1521936/91‑1 (MQ=255)
                                       atataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa                                                     <  1:93274/91‑1 (MQ=255)
                                          taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATTGTGAGGTATATGATAATaaaaaa                                                  <  1:400023/91‑1 (MQ=255)
                                          taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaa                                                  >  1:1272329/1‑91 (MQ=255)
                                            gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                >  1:1342297/1‑91 (MQ=255)
                                            gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                >  2:1839859/1‑91 (MQ=255)
                                            gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                <  1:745282/91‑1 (MQ=255)
                                             cGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAg                                               >  1:3268451/1‑91 (MQ=255)
                                              gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATAATAATAAAAAAAAGc                                              <  2:767540/91‑1 (MQ=255)
                                                tGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCt                                            >  1:1001048/1‑91 (MQ=255)
                                                tGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCt                                            >  1:3687262/1‑91 (MQ=255)
                                                  gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                          <  1:1481924/91‑1 (MQ=255)
                                                  gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                          >  1:85980/1‑91 (MQ=255)
                                                   tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgtt                                         <  1:1982019/91‑1 (MQ=255)
                                                    cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                        >  1:2538084/1‑91 (MQ=255)
                                                     aTACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttgt                                       <  1:1398631/91‑1 (MQ=255)
                                                      tACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTc                                      <  2:1170613/91‑1 (MQ=255)
                                                        cccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcac                                    <  1:2862222/91‑1 (MQ=255)
                                                           aTTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAAAAAAAAAAAAAGCCTGTTGTCACAAt                                 <  2:2909136/91‑1 (MQ=255)
                                                            ttCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                <  1:2455469/91‑1 (MQ=255)
                                                              cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                              <  2:1574105/91‑1 (MQ=255)
                                                               acaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTc                             <  1:456553/91‑1 (MQ=255)
                                                                caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                            <  1:1556345/91‑1 (MQ=255)
                                                                caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                            >  2:1060431/1‑91 (MQ=255)
                                                                caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTCCAATAGCGCGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGCCa                            >  2:622840/1‑91 (MQ=255)
                                                                caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTGTTCAATAGTGCGGTACATGATAATAAAAAAAAGCCTGTTGTCACAATGGCCa                            >  2:3819451/1‑91 (MQ=255)
                                                                 aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                           >  1:1188798/1‑91 (MQ=255)
                                                                 aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                           >  1:3162711/1‑91 (MQ=255)
                                                                  gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                          >  2:3542832/1‑91 (MQ=255)
                                                                  gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                          >  1:3264027/1‑91 (MQ=255)
                                                                    aaCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa                        >  2:1627004/1‑91 (MQ=255)
                                                                    aaCAATCCACACCATAGAAAAGAGTATATAACCTATTAAATAGTGAGGAAAAAGAAAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGa                        <  2:2219415/91‑1 (MQ=255)
                                                                     aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                       <  1:3931816/91‑1 (MQ=255)
                                                                     aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                       <  1:1626344/91‑1 (MQ=255)
                                                                     aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                       >  1:767386/1‑91 (MQ=255)
                                                                      tataccTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca                      >  2:1039041/5‑91 (MQ=17)
                                                                      tataccTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca                      >  2:1930911/5‑91 (MQ=17)
                                                                       aaTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTAAATGAAAAAAAAAAAAAGCCTGTTGTCACAATGGTCATAGacac                     <  2:715797/91‑1 (MQ=255)
                                                                        aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                    <  2:395789/91‑1 (MQ=255)
                                                                          ccTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGAc                  >  1:581222/1‑91 (MQ=255)
                                                                          ccTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGAc                  <  1:203218/91‑1 (MQ=255)
                                                                           cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                 >  1:2099861/1‑91 (MQ=255)
                                                                              aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc              >  1:2575786/1‑91 (MQ=255)
                                                                              aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc              >  2:2312320/1‑91 (MQ=255)
                                                                                    gAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAACAGCCTGTTGTCACAATGGTCATCGACACGACATAATTTcaa        >  2:1015031/1‑91 (MQ=21)
                                                                                    gAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa        >  1:1637194/1‑91 (MQ=37)
                                                                                     aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg       <  1:2938895/91‑1 (MQ=255)
                                                                                      aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAgg      <  2:2426737/91‑1 (MQ=255)
                                                                                      aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAgg      >  2:2602955/1‑91 (MQ=255)
                                                                                       aaGAGTATATAACCTTTTCAATAGTGAGGTATAAGAAAAAAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGt     <  2:1525293/91‑1 (MQ=255)
                                                                                          agTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTc  >  1:1802513/1‑91 (MQ=255)
                                                                                          |                                                                                          
CAAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTC  >  USA300TCH1516_ALE/253893‑254073

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: