Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,502,542 T→G Q648H (CAA→CAC fryA ← putative PTS enzyme: Hpr, enzyme I and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,502,5420TG93.1% 387.6 / 16.0 130Q648H (CAA→CACfryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base T (4/5);  new base G (59/62);  total (63/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 = 6.62e-01

Only 100 of 130 total aligned reads displayed.

CTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCA  >  NC_000913/2502445‑2502638
                                                                                                 |                                                                                                
cttcCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCAGACACGCCTCGCTATCCAGGTg                                                                                                >  1:12009/1‑100 (MQ=255)
 ttcCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGa                                                                                               >  2:3202903/1‑100 (MQ=255)
 ttcCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGa                                                                                               >  1:17627/1‑100 (MQ=255)
  tcCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGAc                                                                                              <  1:1505118/100‑1 (MQ=255)
  tcCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGAc                                                                                              <  2:1882993/100‑1 (MQ=255)
     ggCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGaa                                                                                           >  2:2699656/1‑100 (MQ=255)
       cGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGc                                                                                         >  2:2490872/1‑100 (MQ=255)
        gTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCt                                                                                        <  2:3324408/100‑1 (MQ=255)
        gTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCt                                                                                        >  2:1132714/1‑100 (MQ=255)
         tAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTg                                                                                       <  2:2080937/100‑1 (MQ=255)
          aaaGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTgg                                                                                      >  2:3261849/1‑100 (MQ=255)
           aaGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGc                                                                                     >  2:3872694/1‑100 (MQ=255)
           aaGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGc                                                                                     <  1:2805022/100‑1 (MQ=255)
           aaGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGc                                                                                     <  1:2803671/100‑1 (MQ=255)
           aaGGCGGTGAGTAACGCTTCAATTGCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGc                                                                                     <  2:891560/100‑1 (MQ=255)
              gCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCttt                                                                                  <  1:2640502/100‑1 (MQ=255)
              gCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCttt                                                                                  <  2:3547091/100‑1 (MQ=255)
                ggTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCTTTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTc                                                                                <  2:748063/100‑1 (MQ=255)
                 gTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCa                                                                               <  1:989044/100‑1 (MQ=255)
                 gTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCa                                                                               >  2:1518564/1‑100 (MQ=255)
                   gAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCAcc                                                                             <  2:2368610/100‑1 (MQ=255)
                   gAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCAcc                                                                             <  2:2369302/100‑1 (MQ=255)
                    aGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCAccg                                                                            >  2:2078397/1‑100 (MQ=255)
                      tAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCAccgcc                                                                          <  2:3278871/100‑1 (MQ=255)
                      tAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCAccgcc                                                                          <  1:3883642/100‑1 (MQ=255)
                        aCGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCgg                                                                        <  1:703898/100‑1 (MQ=255)
                        aCGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCgg                                                                        <  2:3357241/100‑1 (MQ=255)
                        aCGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCgg                                                                        <  1:216792/100‑1 (MQ=255)
                         cGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGa                                                                       >  1:3882977/1‑100 (MQ=255)
                           cTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAAt                                                                     <  1:2827047/100‑1 (MQ=255)
                           cTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAAt                                                                     <  2:746906/100‑1 (MQ=255)
                             tCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATAc                                                                   >  2:2724131/1‑100 (MQ=255)
                             tCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATAc                                                                   >  2:1431465/1‑100 (MQ=255)
                               aaTTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATAcgc                                                                 <  1:4045251/100‑1 (MQ=255)
                               aaTTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATAcgc                                                                 >  1:3698180/1‑100 (MQ=255)
                                  ttCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCggg                                                              <  1:1690345/100‑1 (MQ=255)
                                  ttCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCggg                                                              <  1:547554/100‑1 (MQ=255)
                                  ttCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCggg                                                              <  1:986619/100‑1 (MQ=255)
                                  ttCCTGCGCACCGCGGCCCCCACATGCCTGACGTAACGGTTCCCGACACGACGCCCCGTCCAGGGGGCACAGCTGGCTTTTCACCGCCGGTATACGCggg                                                              <  2:2903106/100‑1 (MQ=255)
                                    ccTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCt                                                            <  1:83579/100‑1 (MQ=255)
                                    ccTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCt                                                            <  2:338571/100‑1 (MQ=255)
                                     cTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCta                                                           >  1:364669/1‑100 (MQ=255)
                                     cTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCta                                                           >  1:638802/1‑100 (MQ=255)
                                     cTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCta                                                           >  1:640404/1‑100 (MQ=255)
                                      tGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtac                                                          >  1:3337886/1‑100 (MQ=255)
                                       gcgcACTGCGGCATTCACATGCCGCGCGCGCTAGTTCACGACACGCCTCGCTATCCCGGTGACGAAGTTGGCTTTTCACCGCCGGAATTCGCGGGCtcct                                                         <  2:452947/100‑1 (MQ=255)
                                         gcACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactca                                                       <  2:1660132/100‑1 (MQ=255)
                                            cTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactcatac                                                    <  1:2616665/100‑1 (MQ=255)
                                             tGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactcatact                                                   <  2:305522/100‑1 (MQ=255)
                                             tGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactcatact                                                   <  2:3966424/100‑1 (MQ=255)
                                              gCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactcatactc                                                  <  1:3796448/100‑1 (MQ=255)
                                               cGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactcatactca                                                 <  1:3427676/100‑1 (MQ=255)
                                                ggCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAg                                                >  1:454093/1‑100 (MQ=255)
                                                 gCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGc                                               >  2:1357375/1‑100 (MQ=255)
                                                 gCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGc                                               <  1:889340/100‑1 (MQ=255)
                                                 gCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCGCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCCGc                                               <  1:1335109/100‑1 (MQ=255)
                                                  cATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCt                                              >  1:2606679/1‑100 (MQ=255)
                                                   aTTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTc                                             >  1:1471630/1‑100 (MQ=255)
                                                   aTTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTc                                             >  2:2938960/1‑100 (MQ=255)
                                                   aTTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTc                                             >  2:1784403/1‑100 (MQ=255)
                                                     tCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCAt                                           >  1:3125346/1‑100 (MQ=255)
                                                      cacaTGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATc                                          >  1:1911025/1‑100 (MQ=255)
                                                       acaTGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATcc                                         <  2:3867396/100‑1 (MQ=255)
                                                       acaTGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATcc                                         <  1:2248721/100‑1 (MQ=255)
                                                        caagCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTAAAGAAGCTGACTCTCCAACGCCGGAATAAGCGGGCTAATCATACGCAGCACATCCa                                        >  2:2136326/4‑100 (MQ=255)
                                                        caTGCCTGACTGGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATCCTCAGCTCATCCa                                        >  1:612553/1‑100 (MQ=255)
                                                         aTGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAg                                       <  1:2766065/100‑1 (MQ=255)
                                                          tGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAgg                                      >  2:1923324/1‑100 (MQ=255)
                                                           gCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGgc                                     >  1:1575219/1‑100 (MQ=255)
                                                           gCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGgc                                     <  2:3034322/100‑1 (MQ=255)
                                                            ccTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGgcc                                    <  2:82167/100‑1 (MQ=255)
                                                             cTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGgccc                                   >  2:1908768/1‑100 (MQ=255)
                                                                aCGCGCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGgcccagc                                <  2:470726/100‑1 (MQ=255)
                                                                 cgcgCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGgcccagcc                               >  2:3242654/1‑100 (MQ=255)
                                                                 cgcgCCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGgcccagcc                               <  1:1061304/100‑1 (MQ=255)
                                                                    gcCAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCaa                            >  1:2287309/1‑100 (MQ=255)
                                                                     ccAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAg                           <  1:863573/100‑1 (MQ=255)
                                                                      cAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGc                          >  1:2710832/1‑100 (MQ=255)
                                                                      cAGTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGc                          <  1:495162/100‑1 (MQ=255)
                                                                        gTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAg                        >  1:869291/1‑100 (MQ=255)
                                                                        gTTCCCGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAg                        <  1:2129142/100‑1 (MQ=255)
                                                                            ccGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGc                    >  2:2851733/1‑100 (MQ=255)
                                                                            ccGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGc                    >  1:3225309/1‑100 (MQ=255)
                                                                            ccGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTAAGCTCATCCAGGCCCAGCCCAAGCAGTAGc                    >  1:3887988/1‑100 (MQ=255)
                                                                             cGACACGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCg                   >  1:3721528/1‑100 (MQ=255)
                                                                              gACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAGGCGGGCTACTCATCCTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCgg                  >  2:3938219/1‑100 (MQ=255)
                                                                               acacGCCGCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGc                 <  1:888235/100‑1 (MQ=255)
                                                                                cacGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCa                <  1:3283809/100‑1 (MQ=255)
                                                                                cacGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCa                <  1:4040803/100‑1 (MQ=255)
                                                                                  cGCCTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGa              >  1:3206203/1‑100 (MQ=255)
                                                                                    ccTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATa            <  1:2697198/100‑1 (MQ=255)
                                                                                    ccTCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATa            >  2:2385089/1‑100 (MQ=255)
                                                                                      tCGCTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAAc          <  2:633397/100‑1 (MQ=255)
                                                                                         cTATCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGc       >  1:1350123/1‑100 (MQ=255)
                                                                                           aTCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCtt     <  2:258121/100‑1 (MQ=255)
                                                                                            tCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCttt    <  2:1679304/100‑1 (MQ=255)
                                                                                            tCCAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCttt    >  2:673076/1‑100 (MQ=255)
                                                                                             ccAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCGCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTc   <  1:2571687/100‑1 (MQ=255)
                                                                                              cAGGTGACGAAGCTGGCTTTTCATCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCa  >  1:4031133/1‑100 (MQ=255)
                                                                                              cAGGTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCa  >  2:77785/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCATACTCAGCTCATCCAGGCCCAGCCCAAGCAGTAGCGGCAGATAACGGCTTTCA  >  NC_000913/2502445‑2502638

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