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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,700,7830GC 5.5% 949.4 / 32.8 274P625A (CCC→GCC) ‡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.79e-01

Only 100 of 281 total aligned reads displayed.

GATGCTGCCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGCCGTGAC  >  NC_002947/4700634‑4700932
                                                                                                                                                     |                                                                                                                                                     
 aTGCTGCCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGc                                                                                                                                                      >  1:927873/1‑150 (MQ=255)
 aTGCTGCCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGc                                                                                                                                                      >  1:2536374/1‑150 (MQ=255)
     tgcCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGc                                                                                                                                                  <  1:80365/150‑1 (MQ=255)
     tgcCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGc                                                                                                                                                  <  1:1405257/150‑1 (MQ=255)
     tgcCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGc                                                                                                                                                  >  2:4380721/1‑150 (MQ=255)
     tgcCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGc                                                                                                                                                  >  1:4322935/1‑150 (MQ=255)
       ccACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGAGGGCCAGGGCCTGGGCCAGGGGCTGGCCt                                                                                                                                                >  1:4228224/1‑150 (MQ=255)
        cACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTa                                                                                                                                               >  1:1434543/1‑150 (MQ=255)
         aCGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCTGTCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCAGGCCAGCCCCTGGCCCAGGCCCTGGcccgc                                                                                                                                              <  1:515423/150‑4 (MQ=255)
         aCGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCAGGCCAGCCCCTGGCCCAGGCCCTGGcccgc                                                                                                                                              <  1:3649596/150‑4 (MQ=255)
             aTGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCg                                                                                                                                          >  1:1063594/1‑150 (MQ=255)
             aTGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCg                                                                                                                                          >  1:2872089/1‑150 (MQ=255)
             aTGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCg                                                                                                                                          >  1:1538368/1‑150 (MQ=255)
                 cGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGa                                                                                                                                      <  2:1849931/150‑1 (MQ=255)
                   gTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGcctggccc                                                                                                                                    <  1:2512696/150‑8 (MQ=255)
                   gTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGc                                                                                                                                    <  2:3987251/150‑1 (MQ=255)
                    tCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCa                                                                                                                                   <  2:1789591/150‑1 (MQ=255)
                      aaCGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCTGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAgg                                                                                                                                 >  1:566369/1‑150 (MQ=255)
                      aaCGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAgg                                                                                                                                 >  2:3204663/1‑150 (MQ=255)
                           cGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCg                                                                                                                            <  1:2070254/150‑1 (MQ=255)
                             tGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGGTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGc                                                                                                                          <  1:2756042/150‑1 (MQ=255)
                              gcgcAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCg                                                                                                                         >  2:450268/1‑150 (MQ=255)
                               cgcAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGt                                                                                                                        >  1:2896888/1‑150 (MQ=255)
                                 cAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTgg                                                                                                                      <  1:2169319/150‑1 (MQ=255)
                                     tCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCg                                                                                                                  >  2:2467881/1‑150 (MQ=255)
                                     tCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCg                                                                                                                  >  1:2992367/1‑150 (MQ=255)
                                      cATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGc                                                                                                                 <  1:2352657/150‑1 (MQ=255)
                                       aTGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCt                                                                                                                >  2:1076756/1‑150 (MQ=255)
                                       aTGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCt                                                                                                                >  1:2959299/1‑150 (MQ=255)
                                          ggATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAg                                                                                                             >  1:3107342/1‑150 (MQ=255)
                                             tACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCgg                                                                                                          >  1:1211463/1‑150 (MQ=255)
                                              aCGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGc                                                                                                         >  1:917867/1‑150 (MQ=255)
                                               cGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCt                                                                                                        <  2:623004/150‑1 (MQ=255)
                                                ggAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTg                                                                                                       >  1:4496757/1‑150 (MQ=255)
                                                 gAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGc                                                                                                      >  1:3042776/1‑150 (MQ=255)
                                                     ccAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGt                                                                                                  >  1:3276078/1‑150 (MQ=255)
                                                      cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                                 <  1:74150/150‑1 (MQ=255)
                                                      cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                                 <  1:1646180/150‑1 (MQ=255)
                                                      cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                                 >  1:1799930/1‑150 (MQ=255)
                                                      cagcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGtt                                                                                                 <  1:3183774/150‑1 (MQ=255)
                                                       agcagGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTg                                                                                                >  1:1521949/1‑150 (MQ=255)
                                                            ggCGCTGCGCAGTTGCTCGGTTTCTCAGTGCACGCGTGCGGCTTCCAGCTGTTCAGCAAGGCGTGCTCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGTGCGGTGAGCAGGCGCCGGAGTGGGTCGCTGAGCTGCTGCCCTTTGCGCgg                                                                                           >  1:3884201/1‑150 (MQ=255)
                                                            ggCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCgg                                                                                           >  2:726027/1‑150 (MQ=255)
                                                            ggCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCgg                                                                                           <  1:597208/150‑1 (MQ=255)
                                                              cgcTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGgc                                                                                         <  2:1126189/150‑1 (MQ=255)
                                                                 tGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCCGGGCCTGGGCCAGGGGCTGGCCTCGGGCGGTGAGCGGGCGCCGGCGTGGGTCGCTGGGCGGCTGCCCGTTGCGCGGgcgcc                                                                                      >  2:3273949/1‑149 (MQ=255)
                                                                   cgcAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGcaca                                                                                    <  1:2694125/150‑1 (MQ=255)
                                                                    gcAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCacac                                                                                   <  2:2921732/150‑1 (MQ=255)
                                                                     cAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACAcc                                                                                  >  2:3913608/1‑150 (MQ=255)
                                                                            tCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCagcagc                                                                           >  1:4036411/1‑150 (MQ=255)
                                                                             cGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAgcg                                                                          >  1:3517881/1‑150 (MQ=255)
                                                                                tttCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCa                                                                       >  1:1578985/1‑150 (MQ=255)
                                                                                tttCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCa                                                                       >  1:3221657/1‑150 (MQ=255)
                                                                                  tCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGc                                                                     <  2:443061/150‑1 (MQ=255)
                                                                                  tCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGc                                                                     >  1:1124661/1‑150 (MQ=255)
                                                                                           aGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCg                                                            >  2:2576645/1‑150 (MQ=255)
                                                                                            ggCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGt                                                           <  2:463775/150‑1 (MQ=255)
                                                                                                  cGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCg                                                     <  1:2025028/150‑1 (MQ=255)
                                                                                                      ttCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCaa                                                 >  1:711928/1‑150 (MQ=255)
                                                                                                      ttCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCaa                                                 >  1:75837/1‑150 (MQ=255)
                                                                                                        ccAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATg                                               <  2:2637129/150‑1 (MQ=255)
                                                                                                        ccAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATg                                               <  2:1408575/150‑1 (MQ=255)
                                                                                                        ccAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATg                                               >  1:108353/1‑150 (MQ=255)
                                                                                                        ccAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATg                                               <  1:2927351/150‑1 (MQ=255)
                                                                                                          aGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGc                                             >  2:3630697/1‑150 (MQ=255)
                                                                                                          aGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGc                                             >  1:4150666/1‑150 (MQ=255)
                                                                                                           gCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGcc                                            >  1:1371068/1‑150 (MQ=255)
                                                                                                           gCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGcc                                            >  1:1563398/1‑150 (MQ=255)
                                                                                                                    ccAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccaccacc                                   >  2:773540/1‑150 (MQ=255)
                                                                                                                    ccAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGccaccaccacc                                   >  2:4083120/1‑150 (MQ=255)
                                                                                                                    ccAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAACGGccaccaccacc                                   >  2:382903/1‑150 (MQ=255)
                                                                                                                       ggCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGc                                <  1:3708953/150‑1 (MQ=255)
                                                                                                                       ggCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGc                                <  2:1946365/150‑1 (MQ=255)
                                                                                                                          gTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGc                             <  2:2991218/150‑1 (MQ=255)
                                                                                                                           tGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGcc                            >  2:1222109/1‑150 (MQ=255)
                                                                                                                               cgGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTa                        <  1:2841056/150‑1 (MQ=255)
                                                                                                                               cgGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTa                        <  2:3526579/150‑1 (MQ=255)
                                                                                                                                gggCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAg                       <  2:1431159/150‑1 (MQ=255)
                                                                                                                                gggCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAg                       <  2:157319/150‑1 (MQ=255)
                                                                                                                                  gCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGc                     >  1:4099357/1‑150 (MQ=255)
                                                                                                                                  gCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGc                     <  1:3529077/150‑1 (MQ=255)
                                                                                                                                    cAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAg                   >  1:4405457/1‑150 (MQ=255)
                                                                                                                                    cAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAg                   <  1:4478256/150‑1 (MQ=255)
                                                                                                                                    accgcccTAGGCCAGCCCCTGGCCCAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCCCCACCACCAGCGGCCATTAGGccg                   >  2:3387235/6‑148 (MQ=255)
                                                                                                                                     aGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCCTTAGGCAgg                  >  1:193887/1‑150 (MQ=255)
                                                                                                                                     aGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAgg                  <  1:622379/150‑1 (MQ=255)
                                                                                                                                      gggCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAggg                 <  1:1793996/150‑1 (MQ=255)
                                                                                                                                        cccTAGGCCAGCCCCTGGCCTAGGGCGGTGAACAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTa               >  2:1488509/2‑150 (MQ=255)
                                                                                                                                         ccTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTAt              >  2:3077961/1‑150 (MQ=255)
                                                                                                                                         ccTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTAt              <  1:1621938/150‑1 (MQ=255)
                                                                                                                                          cTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATc             >  2:2322958/1‑150 (MQ=255)
                                                                                                                                          cTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGCGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATc             <  1:1877292/150‑1 (MQ=255)
                                                                                                                                          cTGGGCCAGGGGCGGGCCTAGGGCGGTCAGCAGCCGCCGGCGGGGGACGCTGAGCGGCTGCCCGTCGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGACAACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTaac             <  2:566369/150‑3 (MQ=255)
                                                                                                                                               ccAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGc        <  2:2283429/150‑1 (MQ=255)
                                                                                                                                               ccAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGc        >  1:4277044/1‑150 (MQ=255)
                                                                                                                                               ccAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGc        >  1:728888/1‑150 (MQ=255)
                                                                                                                                                cAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTCTCGCTGcc       >  2:1018402/1‑150 (MQ=255)
                                                                                                                                                cAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGcc       <  2:278124/150‑1 (MQ=255)
                                                                                                                                                  ggggCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGCCGt     >  2:1269851/1‑150 (MQ=255)
                                                                                                                                                     gCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGCCGTGAc  >  2:2385365/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGCTGCCACGCATGGCGGTCAACGGCGTGCGCAAGTCATGGGATACGGAAGCCAGCAGGGCGCTGCGCAGTTGCTCGGTTTCGCCGTGCAGGCGTGCGGCTTCCAGCTGTTCAGCCAGGCGTGCGCGGGCCAGGGCCTGGGCCAGGGGCTGGCCTAGGGCGGTGAGCAGGCGCCGGCGTGGGTCGCTGAGCGGCTGCCCGTTGCGCGGGCGCACACCCAGCAGCGCCAGCGGCTGATCGTCTACCGCCAATGGCCACCACCACCAGCGGCCATTAGGCAGGGTATCGCTGCCGTGAC  >  NC_002947/4700634‑4700932

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: ‑