Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 253,984 G→T 9.5% intergenic (+197/+165) 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,9840GT 9.5% 819.1 / 50.4 243intergenic (+197/+165)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (124/95);  new base T (8/15);  total (132/111)
Fisher's exact test for biased strand distribution p-value = 5.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 257 total aligned reads displayed.

AAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCT  >  USA300TCH1516_ALE/253894‑254074
                                                                                          |                                                                                          
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)
 aGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGt                                                                                           >  1:469445/1‑91 (MQ=255)
  gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGta                                                                                          <  1:3689657/91‑1 (MQ=255)
  gTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGta                                                                                          <  1:2795328/91‑1 (MQ=255)
   tCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtat                                                                                         <  2:3241592/91‑1 (MQ=255)
     aTTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatat                                                                                       >  1:778566/1‑91 (MQ=255)
     aTTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatat                                                                                       >  1:968139/1‑91 (MQ=255)
      ttgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGGTtatata                                                                                      <  1:236489/91‑1 (MQ=255)
      ttgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatata                                                                                      >  1:1312388/1‑91 (MQ=255)
      ttgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatata                                                                                      <  2:308952/91‑1 (MQ=255)
      ttgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGtatata                                                                                      >  1:2561127/1‑91 (MQ=255)
       tgtttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATaa                                                                                     >  2:1678741/1‑91 (MQ=255)
         tttgtAAAGAAGGTGTACTTTATGTAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                   <  2:2839822/91‑1 (MQ=255)
         tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGGTtatatactc                                                                                   <  1:1854553/91‑4 (MQ=255)
         tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGGTtatatactc                                                                                   <  1:1208220/91‑4 (MQ=255)
         tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                   <  2:2197397/91‑1 (MQ=255)
         tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                   <  1:1129276/91‑1 (MQ=255)
           tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                 <  2:3665497/91‑1 (MQ=255)
            gtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCCTACCCATCCACACTAACAATCCTCACCATTGAAAAGCGTATATAACCttc                                                                                >  2:3834521/1‑90 (MQ=255)
             tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCtttt                                                                               <  2:1547865/91‑1 (MQ=37)
             tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtttt                                                                               <  1:3106763/91‑1 (MQ=255)
              aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                              >  2:1783173/1‑91 (MQ=255)
              aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                              <  1:3493129/91‑1 (MQ=255)
              aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                              >  1:407/1‑91 (MQ=255)
              aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                              <  2:3205917/91‑1 (MQ=255)
               aagaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCa                                                                             <  2:551134/91‑1 (MQ=255)
                agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCaa                                                                            <  1:3534217/91‑1 (MQ=25)
                agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCaa                                                                            <  1:222089/91‑1 (MQ=25)
                agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCaa                                                                            <  2:3903558/91‑1 (MQ=25)
                  aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAtg                                                                          <  1:1464236/91‑2 (MQ=25)
                  aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                                          >  2:3132607/1‑91 (MQ=255)
                   agGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                                         <  1:3815048/91‑1 (MQ=255)
                   agGTGTACTTTAGATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAtgg                                                                         <  2:1635880/91‑3 (MQ=21)
                    ggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTAACAGTAACAATCCTCACCATTGAAACGAGTATATAACCTTTTCAATAGt                                                                        >  2:3006756/1‑91 (MQ=255)
                     gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTg                                                                       >  1:3938669/1‑91 (MQ=255)
                      tgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGa                                                                      <  1:565478/91‑1 (MQ=255)
                       gtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAg                                                                     <  2:3860964/91‑1 (MQ=255)
                        tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAgg                                                                    >  1:753780/1‑91 (MQ=255)
                          cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                  >  1:1912619/1‑91 (MQ=255)
                          cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                  <  2:2958649/91‑1 (MQ=255)
                             tatataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTatat                                                               <  2:1328875/91‑1 (MQ=255)
                              atctaAGTATATAGCGATGGTCATACCCATCCACAGTACCAATCCTCACCATTGCAAACCGCCTATAACCTTTTCAACCGTGCGGTATATg                                                              >  2:442158/4‑91 (MQ=255)
                                  aaGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataa                                                          <  2:2769104/91‑1 (MQ=255)
                                   aGTATATAGCGATGGTCATACCCATTCACAGTACCAATCCTCACCATTGCAAAGAGTACATAACCTTTTCAATAGTGAGGTATATGataat                                                         >  2:144389/1‑91 (MQ=255)
                                   aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                         >  2:1135667/1‑91 (MQ=255)
                                   aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                         >  1:3845178/1‑91 (MQ=255)
                                    gTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataata                                                        <  1:1502655/91‑1 (MQ=255)
                                     tatataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAtaataa                                                       <  1:1521936/91‑1 (MQ=255)
                                     tatataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAtaataa                                                       <  1:2770952/91‑1 (MQ=255)
                                      atataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaa                                                      <  1:93274/91‑1 (MQ=255)
                                         taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaa                                                   >  1:1272329/1‑91 (MQ=255)
                                         taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaa                                                   <  1:815851/91‑1 (MQ=255)
                                         taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaa                                                   <  1:1449564/91‑1 (MQ=255)
                                           gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                 >  1:1342297/1‑91 (MQ=255)
                                            cGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAg                                                >  1:3268451/1‑91 (MQ=255)
                                             gATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATAATAATAAAAAAAAGc                                               <  2:767540/91‑1 (MQ=255)
                                               tGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCt                                             >  1:1001048/1‑91 (MQ=255)
                                                 gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                           >  1:1556739/1‑91 (MQ=255)
                                                 gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                           <  1:1803378/91‑1 (MQ=255)
                                                  tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgtt                                          <  1:1982019/91‑1 (MQ=255)
                                                     tACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTc                                       >  1:2974777/1‑91 (MQ=255)
                                                     tACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTc                                       <  2:1170613/91‑1 (MQ=255)
                                                       cccATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTcac                                     <  1:2862612/91‑1 (MQ=255)
                                                          aTTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAAAAAAAAAAAAAGCCTGTTGTCACAAt                                  <  2:2909136/91‑1 (MQ=255)
                                                           ttCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATg                                 >  1:714582/1‑91 (MQ=255)
                                                             cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                               <  1:3842050/91‑1 (MQ=255)
                                                             cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                               <  2:1574105/91‑1 (MQ=255)
                                                               caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                             >  2:225243/1‑91 (MQ=255)
                                                               caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                             >  2:1060431/1‑91 (MQ=255)
                                                               caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                             <  1:675515/91‑1 (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:2479701/1‑91 (MQ=255)
                                                                aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                            >  1:3162711/1‑91 (MQ=255)
                                                                 gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                           >  1:884439/1‑91 (MQ=255)
                                                                 gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                           >  1:3264027/1‑91 (MQ=255)
                                                                 gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTAAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                           >  2:3542539/1‑91 (MQ=255)
                                                                    aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                        <  1:1659342/91‑1 (MQ=255)
                                                                    aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                        <  1:1626344/91‑1 (MQ=255)
                                                                     tataccTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca                       >  2:1930911/5‑91 (MQ=17)
                                                                     tataccTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca                       >  2:1039041/5‑91 (MQ=17)
                                                                      aaTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGacac                      <  2:2077378/91‑1 (MQ=255)
                                                                      aaTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTAAATGAAAAAAAAAAAAAGCCTGTTGTCACAATGGTCATAGacac                      <  2:715797/91‑1 (MQ=255)
                                                                       aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                     >  2:3695769/1‑91 (MQ=255)
                                                                       aTCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACg                     <  2:395789/91‑1 (MQ=255)
                                                                         ccTCACCATTGAAGAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAGAAAAAGCCTGTTGTCACGATGGTCATAGACACGAc                   <  1:203268/91‑1 (MQ=255)
                                                                          cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                  >  1:2099861/1‑91 (MQ=255)
                                                                             aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc               >  1:2575861/1‑91 (MQ=255)
                                                                             aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc               >  2:2312320/1‑91 (MQ=255)
                                                                                  tGAAAAGAGTATATAACCTGTGCAATATGGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaa          <  2:1844845/91‑1 (MQ=255)
                                                                                   gAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAACAGCCTGTTGTCACAATGGTCATCGACACGACATAATTTcaa         >  2:1015031/1‑91 (MQ=21)
                                                                                     aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAgg       <  2:2426737/91‑1 (MQ=255)
                                                                                     aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAgg       >  2:2604188/1‑91 (MQ=255)
                                                                                      aaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGt      >  1:315808/1‑91 (MQ=255)
                                                                                          gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt  >  2:2598440/1‑91 (MQ=255)
                                                                                          |                                                                                          
AAGTCATTGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCT  >  USA300TCH1516_ALE/253894‑254074

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: