Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,700,782 G→C 5.5% P625R (CCC→CGC) ‡ kdpD ← sensor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,700,7820GC 5.5% 947.8 / 32.8 274P625R (CCC→CGC) ‡kdpDsensor protein
Reads supporting (aligned to +/- strand):  ref base G (143/116);  new base C (3/12);  total (146/128)
Fisher's exact test for biased strand distribution p-value = 1.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 279 total aligned reads displayed.

GATGCTGCCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGCCGT  >  NC_002947/4700634‑4700929
                                                                                                                                                    |                                                                                                                                                   
 aTGCTGCCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGc                                                                                                                                                   >  1:2536374/1‑150 (MQ=255)
     tgcCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGc                                                                                                                                               <  2:1432743/150‑1 (MQ=255)
     tgcCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGc                                                                                                                                               <  2:2288194/150‑1 (MQ=255)
     tgcCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGc                                                                                                                                               >  1:3211459/1‑150 (MQ=255)
        cACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTa                                                                                                                                            >  1:1434543/1‑150 (MQ=255)
         aCGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCAGGCCAGCCCCTGGCCCAGGCCCTGGcccgc                                                                                                                                           <  1:3649596/150‑4 (MQ=255)
           gCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAggg                                                                                                                                         >  2:1621938/1‑150 (MQ=255)
            cATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGCCCTAGGCCAGCCCCTGGCCCAggcc                                                                                                                                        <  1:1382992/150‑3 (MQ=255)
             aTGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCg                                                                                                                                       >  1:1063594/1‑150 (MQ=255)
             aTGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCg                                                                                                                                       >  1:1538368/1‑150 (MQ=255)
              tGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGCCCTAGGCCAGCCCCTGGCCCAggccct                                                                                                                                      <  1:925476/150‑5 (MQ=255)
               ggCTGTCAACGGAGTGCGCAAGTCATGTGATACGGAAGACAGCAGGGCGCTTCGAAGTTGCTCGGTTTCGCCGTTCAGGCGTGCGGCTTCCATCTGTTCAGACAGGCTTGCGCGGGCCAGTGCCTGGGCCAGGGGATGGCCTAGGGCGGt                                                                                                                                     <  1:489505/150‑1 (MQ=255)
               ggCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGt                                                                                                                                     <  1:3694625/150‑1 (MQ=255)
               ggCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTAGGCCAGCCCCTGGCCCAggccctg                                                                                                                                     <  2:306861/150‑6 (MQ=255)
                gCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTAGGCCAGCCCCTGGCCCAggccctgg                                                                                                                                    <  2:1335431/150‑7 (MQ=255)
                   gTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGc                                                                                                                                 <  2:3987251/150‑1 (MQ=255)
                    tCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCa                                                                                                                                <  1:3474156/150‑1 (MQ=255)
                    tCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCa                                                                                                                                >  2:3270849/1‑150 (MQ=255)
                    tCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCa                                                                                                                                <  1:2752605/150‑1 (MQ=255)
                      aaCGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGGTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAgg                                                                                                                              >  2:2627211/1‑150 (MQ=255)
                       aCGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGgc                                                                                                                             >  1:4310068/1‑150 (MQ=255)
                       aCGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGgc                                                                                                                             <  2:3599721/150‑1 (MQ=255)
                           cGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCg                                                                                                                         <  1:2070254/150‑1 (MQ=255)
                              gcgcAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCg                                                                                                                      >  2:450268/1‑150 (MQ=255)
                               cgcAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGt                                                                                                                     >  1:2896888/1‑150 (MQ=255)
                               cgcAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGt                                                                                                                     <  1:3936158/150‑1 (MQ=255)
                                gcAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTg                                                                                                                    <  2:4034502/150‑1 (MQ=255)
                                gcAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTg                                                                                                                    <  1:780020/150‑1 (MQ=255)
                                 cAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTgg                                                                                                                   <  1:2169319/150‑1 (MQ=255)
                                     tCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCg                                                                                                               >  2:2467881/1‑150 (MQ=255)
                                        tGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTg                                                                                                            >  1:2338229/1‑150 (MQ=255)
                                        tGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTg                                                                                                            >  1:1864096/1‑150 (MQ=255)
                                          ggATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAg                                                                                                          >  1:96953/1‑150 (MQ=255)
                                          ggATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAg                                                                                                          >  1:833413/1‑150 (MQ=255)
                                          ggATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAg                                                                                                          >  1:3107342/1‑150 (MQ=255)
                                             tACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCgg                                                                                                       >  1:1211463/1‑150 (MQ=255)
                                             tACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCgg                                                                                                       <  1:2927635/150‑1 (MQ=255)
                                               cGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCt                                                                                                     >  1:282100/1‑150 (MQ=255)
                                               cGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCt                                                                                                     <  2:623004/150‑1 (MQ=255)
                                                ggAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTg                                                                                                    >  1:4496718/1‑150 (MQ=255)
                                                ggAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTg                                                                                                    >  2:221417/1‑150 (MQ=255)
                                                    gCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCg                                                                                                <  1:2862609/150‑1 (MQ=255)
                                                      cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                              >  1:1799930/1‑150 (MQ=255)
                                                      cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                              <  1:74020/150‑1 (MQ=255)
                                                      cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                              <  1:3183774/150‑1 (MQ=255)
                                                      cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                              <  2:3786239/150‑1 (MQ=255)
                                                      cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                              <  1:1646180/150‑1 (MQ=255)
                                                       agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTg                                                                                             >  2:1251488/1‑150 (MQ=255)
                                                       agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTg                                                                                             >  1:1521949/1‑150 (MQ=255)
                                                         tagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTgcg                                                                                           >  2:2548591/2‑150 (MQ=255)
                                                            ggCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCgg                                                                                        >  1:621348/1‑150 (MQ=255)
                                                            ggCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCgg                                                                                        >  2:726027/1‑150 (MQ=255)
                                                            ggCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCgg                                                                                        >  1:3012962/1‑150 (MQ=255)
                                                            ggCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCgg                                                                                        <  2:3513468/150‑1 (MQ=255)
                                                                   cgcAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGcaca                                                                                 <  1:728049/150‑1 (MQ=255)
                                                                   cgcAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGcaca                                                                                 <  1:2694125/150‑1 (MQ=255)
                                                                    gcAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCacac                                                                                <  2:2921732/150‑1 (MQ=255)
                                                                      aGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACAccc                                                                              >  2:1847308/1‑150 (MQ=255)
                                                                          gCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCcagca                                                                          <  1:660627/150‑1 (MQ=255)
                                                                            tCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCagcagc                                                                        >  2:3174702/1‑150 (MQ=255)
                                                                            tCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCagcagc                                                                        >  2:4276286/1‑150 (MQ=255)
                                                                                tttCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCa                                                                    >  1:4381626/1‑150 (MQ=255)
                                                                                  tCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGc                                                                  >  1:1124661/1‑150 (MQ=255)
                                                                                    gCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCgg                                                                <  2:1044236/150‑1 (MQ=255)
                                                                                                  cGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCg                                                  <  1:2025028/150‑1 (MQ=255)
                                                                                                    gCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGcc                                                >  2:142011/1‑150 (MQ=255)
                                                                                                       tCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAAt                                             <  2:1137096/150‑1 (MQ=255)
                                                                                                        ccAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATg                                            >  1:108353/1‑150 (MQ=255)
                                                                                                        ccAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATg                                            >  1:592862/1‑150 (MQ=255)
                                                                                                        ccAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATg                                            <  1:4271997/150‑1 (MQ=255)
                                                                                                        ccAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATg                                            >  1:4198575/1‑150 (MQ=255)
                                                                                                        ccAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATg                                            <  1:2927351/150‑1 (MQ=255)
                                                                                                          aGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGc                                          >  2:2375426/1‑150 (MQ=255)
                                                                                                          aGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGc                                          <  2:1078631/150‑1 (MQ=255)
                                                                                                           gCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGcc                                         >  1:1371068/1‑150 (MQ=255)
                                                                                                           gCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGcc                                         >  1:1563398/1‑150 (MQ=255)
                                                                                                            cTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGcca                                        <  1:4147376/150‑1 (MQ=255)
                                                                                                                 cAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccacc                                   >  1:236765/1‑150 (MQ=255)
                                                                                                                  aGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccacca                                  <  1:757594/150‑1 (MQ=255)
                                                                                                                  aGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccacca                                  >  1:4457542/1‑150 (MQ=255)
                                                                                                                  aGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccacca                                  >  1:4095448/1‑150 (MQ=255)
                                                                                                                    ccAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccaccacc                                >  1:3270276/1‑150 (MQ=255)
                                                                                                                    ccAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccaccacc                                >  2:773540/1‑150 (MQ=255)
                                                                                                                     cAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccaccacca                               >  2:1878326/1‑150 (MQ=255)
                                                                                                                       ggCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGc                             <  2:3848489/150‑1 (MQ=255)
                                                                                                                       ggCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGc                             <  1:544337/150‑1 (MQ=255)
                                                                                                                           tGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGcc                         >  1:2984203/1‑150 (MQ=255)
                                                                                                                                gggCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAg                    <  2:2965234/150‑1 (MQ=255)
                                                                                                                                  gCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGc                  <  1:3529077/150‑1 (MQ=255)
                                                                                                                                  gCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATGAGGc                  >  2:4397355/1‑150 (MQ=255)
                                                                                                                                    cAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAg                <  2:1563398/150‑1 (MQ=255)
                                                                                                                                    accgcccTAGGCCAGCCCCTGGCCCAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCCCCACCACCAGCGGCCATTAGGccg                >  2:3387235/6‑148 (MQ=255)
                                                                                                                                     aGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCCTTAGGCAgg               >  1:193887/1‑150 (MQ=255)
                                                                                                                                      gggCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAggg              <  1:1793996/150‑1 (MQ=255)
                                                                                                                                         ccTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTAt           >  2:488188/1‑150 (MQ=255)
                                                                                                                                         ccTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTAt           <  1:1621938/150‑1 (MQ=255)
                                                                                                                                          cTGGGCCAGGGGCGGGCCTAGGGCGGTCAGCAGCCGCCGGCGGGGGACGCTGAGCGGCTGCCCGTCGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGACAACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTaac          <  2:566369/150‑3 (MQ=255)
                                                                                                                                            gggCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGc        >  1:1567612/1‑150 (MQ=255)
                                                                                                                                                cAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGcc    <  2:278124/150‑1 (MQ=255)
                                                                                                                                                  ggggCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAACGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGCCGt  <  1:382903/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATGCTGCCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGCCGT  >  NC_002947/4700634‑4700929

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