Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 253,991 A→C 6.9% intergenic (+204/+158) 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,9910AC 6.9% 803.7 / 36.3 241intergenic (+204/+158)asbF/USA300HOU_RS011253‑dehydroshikimate dehydratase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (124/99);  new base C (8/10);  total (132/109)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

TGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATA  >  USA300TCH1516_ALE/253901‑254081
                                                                                           |                                                                                          
  tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                          <  1:1129276/91‑1 (MQ=255)
  tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                          <  2:2197397/91‑1 (MQ=255)
  tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                          <  2:313742/91‑1 (MQ=255)
  tttgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAAcc                                                                                          <  2:3836560/91‑1 (MQ=255)
    tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                        <  2:3665497/91‑1 (MQ=255)
    tgtAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtt                                                                                        <  2:536961/91‑1 (MQ=255)
      tAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCtttt                                                                                      <  1:3106763/91‑1 (MQ=255)
      tAAAGAAGGTGTACTTTATAGAAGTATATAGCGATGGTCATACCCATTCACAGTAGCAATCCTCACTATTGAAAAGGTTATATACTCtttt                                                                                      <  2:1547909/91‑1 (MQ=25)
       aaaGGAGGTGTACTTTATATAACTATATAGCGATGGTCATACCCATCCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  2:2048330/1‑91 (MQ=255)
       aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  2:1667829/1‑91 (MQ=255)
       aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  1:3154564/1‑91 (MQ=255)
       aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  2:1731353/1‑91 (MQ=255)
       aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     <  1:2293851/91‑1 (MQ=255)
       aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     <  2:3205917/91‑1 (MQ=255)
       aaaGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTc                                                                                     >  1:534285/1‑91 (MQ=255)
         agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCaa                                                                                   <  1:222089/91‑1 (MQ=25)
         agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCaa                                                                                   <  2:3903558/91‑1 (MQ=25)
         agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCaa                                                                                   <  1:3499799/91‑1 (MQ=255)
         agaagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCaa                                                                                   <  1:1508950/91‑1 (MQ=255)
           aagGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATa                                                                                 >  2:724955/1‑91 (MQ=255)
            agGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                                                <  1:3815048/91‑1 (MQ=255)
            agGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAg                                                                                <  1:1453870/91‑1 (MQ=255)
            agGTGTACTTTAGATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACTATTGAAAAGGTTATATACTCTTTTCAAtgg                                                                                <  2:1635880/91‑3 (MQ=21)
             ggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGt                                                                               >  1:692346/1‑91 (MQ=255)
             ggTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTAACAGTAACAATCCTCACCATTGAAACGAGTATATAACCTTTTCAATAGt                                                                               >  2:3006756/1‑91 (MQ=255)
              gtgtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTg                                                                              >  1:3938669/1‑91 (MQ=255)
                gtACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAg                                                                            <  2:3860964/91‑1 (MQ=255)
                 tACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAgg                                                                           >  1:753780/1‑91 (MQ=255)
                   cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                         <  2:2958649/91‑1 (MQ=255)
                   cTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGta                                                                         >  1:1912619/1‑91 (MQ=255)
                   acttATATAAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGt                                                                          >  1:1247476/3‑91 (MQ=255)
                      tatataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTatat                                                                      >  2:271707/1‑91 (MQ=255)
                      tatataAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTatat                                                                      <  2:1328875/91‑1 (MQ=255)
                          taAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGata                                                                  <  1:1340952/91‑1 (MQ=255)
                             aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                >  2:1135667/1‑91 (MQ=255)
                             aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                >  2:3067573/1‑91 (MQ=255)
                             aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                >  1:3845178/1‑91 (MQ=255)
                             aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGataat                                                                >  1:1307334/1‑91 (MQ=255)
                             aGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTAAATGataat                                                                >  2:3696684/1‑91 (MQ=255)
                               tatataGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGAtaataa                                                              <  1:2770952/91‑1 (MQ=255)
                                   taGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATTGTGAGGTATATGATAATaaaaaa                                                          <  1:400023/91‑1 (MQ=255)
                                     gCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATaaaaaaaa                                                        >  2:1839859/1‑91 (MQ=255)
                                      cGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAg                                                       >  1:3268451/1‑91 (MQ=255)
                                         tGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCt                                                    >  1:1001048/1‑91 (MQ=255)
                                           gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                                  <  1:1481924/91‑1 (MQ=255)
                                           gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                                  <  1:1803378/91‑1 (MQ=255)
                                           gTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgt                                                  >  1:3212239/1‑91 (MQ=255)
                                            tCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgtt                                                 <  1:1982019/91‑1 (MQ=255)
                                             cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                                >  2:1591554/1‑91 (MQ=255)
                                             cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                                >  1:996918/1‑91 (MQ=255)
                                             cATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttg                                                >  1:1571862/1‑91 (MQ=255)
                                              aTACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCtgttgt                                               <  1:1398631/91‑1 (MQ=255)
                                               tACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGCGGTATATGATAATAAAAAAAGCCTGTTGTca                                             >  2:474127/1‑91 (MQ=255)
                                               tACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTc                                              <  2:1170613/91‑1 (MQ=255)
                                                       cacaGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGt                                      <  1:3842050/91‑1 (MQ=255)
                                                         caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                    >  2:1059595/1‑91 (MQ=255)
                                                         caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                    <  1:1556345/91‑1 (MQ=255)
                                                         caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                    <  1:675515/91‑1 (MQ=255)
                                                         caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCa                                    >  2:225243/1‑91 (MQ=255)
                                                         caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTCCAATAGCGCGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGCCa                                    >  2:622840/1‑91 (MQ=255)
                                                         caGTAACAATCCTCACCATTGAAAAGAGTATATAACCTGTTCAATAGTGCGGTACATGATAATAAAAAAAAGCCTGTTGTCACAATGGCCa                                    >  2:3819451/1‑91 (MQ=255)
                                                          aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                   >  1:2479701/1‑91 (MQ=255)
                                                          aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                   >  1:1188798/1‑91 (MQ=255)
                                                          aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                   >  2:395607/1‑91 (MQ=255)
                                                          aGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCAt                                   >  1:3162711/1‑91 (MQ=255)
                                                           gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                                  >  2:3542832/1‑91 (MQ=255)
                                                           gTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATa                                  >  1:884439/1‑91 (MQ=255)
                                                            tAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAg                                 >  1:1630185/1‑91 (MQ=255)
                                                              aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                               <  1:1659342/91‑1 (MQ=255)
                                                              aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                               >  2:430866/1‑91 (MQ=255)
                                                              aCAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGac                               >  1:767386/1‑91 (MQ=255)
                                                               tataccTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca                              >  2:1930911/5‑91 (MQ=17)
                                                               tataccTCACTATTGAAAAGGTTATATACTCTTTTCAATGGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGaca                              >  2:1039041/5‑91 (MQ=17)
                                                                   ccTCACCATTGAAGAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAGAAAAAGCCTGTTGTCACGATGGTCATAGACACGAc                          <  1:203268/91‑1 (MQ=255)
                                                                   ccTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGAc                          >  1:2831715/1‑91 (MQ=255)
                                                                    cTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACa                         >  1:2099861/1‑91 (MQ=255)
                                                                     tCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCACAGACACGACAt                        >  1:3575164/1‑91 (MQ=255)
                                                                       aCCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATAc                      >  1:2575786/1‑91 (MQ=255)
                                                                        ccATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACt                     >  1:3464019/1‑91 (MQ=255)
                                                                         cATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACtt                    <  2:2287138/91‑1 (MQ=255)
                                                                             gAAAAGAGTGTATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaa                <  1:3790852/91‑1 (MQ=255)
                                                                              aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg               <  1:2938895/91‑1 (MQ=255)
                                                                              aaaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAg               <  2:812646/91‑1 (MQ=255)
                                                                               aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTaaacg              >  2:2878882/1‑89 (MQ=255)
                                                                               aaaGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAgg              <  2:2426737/91‑1 (MQ=255)
                                                                                aaGAGTATATAACCTTTTCAATAGTGAGGTATAAGAAAAAAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGt             <  2:1525293/91‑1 (MQ=255)
                                                                                    gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt         >  2:2913722/1‑91 (MQ=255)
                                                                                    gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt         >  1:928638/1‑91 (MQ=255)
                                                                                    gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt         >  1:2895063/1‑91 (MQ=255)
                                                                                    gTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCt         >  2:2598440/1‑91 (MQ=255)
                                                                                      atataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGa       >  2:2822941/1‑91 (MQ=255)
                                                                                       tataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaa      >  2:660049/1‑91 (MQ=255)
                                                                                        gtaACCTTTTCAATAGTGAGGTATATGATAATAAAATAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaat     <  2:3910272/90‑1 (MQ=255)
                                                                                        ataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaat     >  2:3466272/1‑91 (MQ=255)
                                                                                        ataACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaat     >  2:264506/1‑91 (MQ=255)
                                                                                         taACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata    <  1:1720042/91‑1 (MQ=255)
                                                                                         taACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata    <  2:3274291/91‑1 (MQ=255)
                                                                                         taACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaata    <  1:561870/91‑1 (MQ=255)
                                                                                           aCCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatata  >  2:1419835/1‑91 (MQ=255)
                                                                                           aCCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGaatata  <  2:1003716/91‑1 (MQ=255)
                                                                                           |                                                                                          
TGTTTGTAAAGAAGGTGTACTTTATATAAGTATATAGCGATGGTCATACCCATTCACAGTAACAATCCTCACCATTGAAAAGAGTATATAACCTTTTCAATAGTGAGGTATATGATAATAAAAAAAAGCCTGTTGTCACAATGGTCATAGACACGACATACTTTAAAGGTTTCTGAATATA  >  USA300TCH1516_ALE/253901‑254081

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: