Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,236,193 T→A 5.8% G415G (GGT→GGA) ‡ PP_2834 → D‑galactarate dehydratase/altronate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,236,1930TA 5.8% 251.8 / 4.1 86G415G (GGT→GGA) ‡PP_2834D‑galactarate dehydratase/altronate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (39/42);  new base A (3/2);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

GACCACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGCGGCGGCAT  >  NC_002947/3236046‑3236342
                                                                                                                                                   |                                                                                                                                                     
gACCACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCa                                                                                                                                                     <  2:2037905/150‑1 (MQ=255)
  ccACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCCCTCAgg                                                                                                                                                   <  2:2331176/150‑1 (MQ=255)
   cACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGt                                                                                                                                                  >  1:1825291/1‑150 (MQ=255)
    aCTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTg                                                                                                                                                 <  1:1347176/150‑1 (MQ=255)
    aCTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTg                                                                                                                                                 <  1:168824/150‑1 (MQ=255)
      tATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGc                                                                                                                                               >  1:2030239/1‑150 (MQ=255)
        tCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCcg                                                                                                                                             >  1:182808/1‑150 (MQ=255)
         ccTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCTGTCTCGGCCACCGGTCAGGTGGCcgg                                                                                                                                            >  2:718621/1‑150 (MQ=255)
           tCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCCCGGCCACCTGACCGGTGGCcgaga                                                                                                                                          <  1:1957648/150‑4 (MQ=255)
           tCGAGAAAGCCCTCGGAGCAGTGCCCAAGGCAGGCTCTAGCAGCCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCTGGGCATCTTCTGCTTTGACACCCCCGGCTACGACCCGGTCTCGGCCACGTGACCGGTGGCcgaga                                                                                                                                          <  2:1736921/150‑4 (MQ=255)
               gaAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAg                                                                                                                                                    <  2:27028/136‑1 (MQ=255)
               gaAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAg                                                                                                                                                    >  1:27028/1‑136 (MQ=255)
               gaAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGCTGCCCGGCGCGGc                                                                                                                                      >  2:12312/1‑150 (MQ=255)
                  aTCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCaa                                                                                                                                   >  2:2273997/1‑150 (MQ=255)
                   tCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAAc                                                                                                                                  >  1:2275411/1‑150 (MQ=255)
                       tGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGa                                                                                                                              >  2:1337871/1‑150 (MQ=255)
                         ggCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATc                                                                                                                            <  1:1912813/150‑1 (MQ=255)
                          gcgcAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCg                                                                                                                           >  2:2278565/1‑150 (MQ=255)
                          gcgcAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCg                                                                                                                           <  1:2164065/150‑1 (MQ=255)
                           cggaGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACTACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGc                                                                                                                          <  2:719418/147‑1 (MQ=255)
                              agaggCCAATGCAGACTCCAACAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGAGGCTGGTGTTCATGGACACCCCCCGCTACTACCCGCTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATTGCCtt                                                                                                                       >  2:1148763/4‑150 (MQ=255)
                                   ccAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGc                                                                                                                                      >  2:204683/1‑130 (MQ=255)
                                   ccAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGc                                                                                                                                      <  1:204683/130‑1 (MQ=255)
                                          aggcTCCAGCAACCTGGTCGATGTGTTTCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAAc                                                                                                                                  >  1:1816134/1‑127 (MQ=255)
                                          aggcTCCAGCAACCTGGTCGATGTGTATCAGGACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAAc                                                                                                                                  <  2:1816134/127‑1 (MQ=255)
                                                cAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTc                                                                                                     <  2:2176957/150‑1 (MQ=255)
                                                 aGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCg                                                                                                    >  1:505278/1‑150 (MQ=255)
                                                         ggTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACgg                                                                                            <  2:1995710/150‑1 (MQ=255)
                                                          gTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGc                                                                                           <  1:2278565/150‑1 (MQ=255)
                                                           tCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCCACACCCGCGGCTCGAGCTACCGCt                                                                                          <  2:31153/150‑1 (MQ=255)
                                                                 tgtATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGccc                                                                                    <  2:182808/150‑1 (MQ=255)
                                                                  gtATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGcccc                                                                                   >  1:1881669/1‑150 (MQ=255)
                                                                   tATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCg                                                                                  <  2:685940/150‑1 (MQ=255)
                                                                   atctaGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCg                                                                                  <  1:553210/146‑1 (MQ=255)
                                                                     tCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGgc                                                                                >  1:1927621/1‑150 (MQ=255)
                                                                     tCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCCGCTCGGCATACGGCGGGGCCCCGgc                                                                                >  2:58512/1‑150 (MQ=255)
                                                                        gTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGcc                                                                             >  1:1946197/1‑150 (MQ=255)
                                                                        gTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGcc                                                                             >  1:2154546/1‑150 (MQ=255)
                                                                        gTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGAGGCAACCTGATCGCCTTCAACACCTGCCGCGGCTCGGCCTCCGGCTGTGCCCCggatat                                                                             >  2:553210/1‑146 (MQ=255)
                                                                         tACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCa                                                                            >  2:1902169/1‑150 (MQ=255)
                                                                           cgcgCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATc                                                                          <  1:1902169/150‑1 (MQ=255)
                                                                            gcgcAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCg                                                                         <  2:1825291/150‑1 (MQ=255)
                                                                             cgcAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGa                                                                        <  1:629986/150‑1 (MQ=255)
                                                                               cAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATc                                                                      <  2:2159233/150‑1 (MQ=255)
                                                                               cAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATc                                                                      <  2:2039447/150‑1 (MQ=255)
                                                                                    cgtgcgtgCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCt                                                                 <  2:2030239/150‑1 (MQ=255)
                                                                                     gtgcgtgcCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCg                                                                                                               <  2:1375302/103‑1 (MQ=255)
                                                                                     gtgcgtgcCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCg                                                                                                               >  1:1375302/1‑103 (MQ=255)
                                                                                     gtgcgtgcCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTg                                                                <  2:2217655/150‑1 (MQ=255)
                                                                                      tgcgtgcCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCACGGCGCCATCGATCAAGCTgg                                                               <  2:2195222/150‑1 (MQ=255)
                                                                                         gtgcCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGcaa                                                            <  2:1946197/150‑1 (MQ=255)
                                                                                            cccAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGcaacca                                                         <  2:2154546/150‑1 (MQ=255)
                                                                                              cAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAAc                                                                                                                                  >  1:1664431/1‑75 (MQ=255)
                                                                                              cAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAAc                                                                                                                                  <  2:1664431/75‑1 (MQ=255)
                                                                                              cAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGcaaccaac                                                       <  1:673733/150‑1 (MQ=255)
                                                                                                    cTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCgg                                                 >  2:2175263/1‑150 (MQ=255)
                                                                                                    cTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCgg                                                 <  2:1904877/150‑1 (MQ=255)
                                                                                                        tgtTTATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGgtgt                                             >  2:534127/1‑150 (MQ=255)
                                                                                                        tgtTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGgtgt                                             <  1:2273997/150‑1 (MQ=255)
                                                                                                          ttCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGa                                                                        >  1:680668/1‑121 (MQ=255)
                                                                                                          ttCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGa                                                                        <  2:680668/121‑1 (MQ=255)
                                                                                                          ttCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTc                                           <  1:1337871/150‑1 (MQ=255)
                                                                                                             aTGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAg                                        <  2:505278/150‑1 (MQ=255)
                                                                                                               ggACACCCCCGGCTACGACACGGTCTCGGCCACCGGTGAGGTGGCCGGCCGGGCCAACCTGAGCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCa                                      >  2:384805/1‑150 (MQ=255)
                                                                                                                  cacCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCa                                   <  1:12312/150‑1 (MQ=255)
                                                                                                                     ccccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTc                                           >  2:1783040/1‑139 (MQ=255)
                                                                                                                     ccccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTc                                           <  1:1783040/139‑1 (MQ=255)
                                                                                                                     ccccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCa                                >  1:1443925/1‑150 (MQ=255)
                                                                                                                       ccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGcc                                                                             >  2:571270/1‑103 (MQ=255)
                                                                                                                       ccGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGcc                                                                             <  1:571270/103‑1 (MQ=255)
                                                                                                                               gACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATg                      >  2:1538281/1‑150 (MQ=255)
                                                                                                                                aCCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATgg                     >  1:985696/1‑150 (MQ=255)
                                                                                                                                 cccGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGa                                         >  2:2308787/1‑129 (MQ=255)
                                                                                                                                 cccGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGa                                         <  1:2308787/129‑1 (MQ=255)
                                                                                                                                  ccGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCgg                                                 >  2:907118/1‑120 (MQ=255)
                                                                                                                                  ccGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCgg                                                 <  1:907118/120‑1 (MQ=255)
                                                                                                                                  ccGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGAc                   <  2:310975/150‑1 (MQ=255)
                                                                                                                                   ccccgccGGCCACCTGACCGGTGGCCGACGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACg                  >  1:488286/7‑150 (MQ=255)
                                                                                                                                      cgccGGCCACCTCACAGGCGGCCGACGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCa               >  2:1375685/4‑150 (MQ=255)
                                                                                                                                          ggCCACCGGTCAGGTGGCGGGGGGGGCCGACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTg           >  2:1960170/1‑150 (MQ=255)
                                                                                                                                          ggCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTg           >  1:1106762/1‑150 (MQ=255)
                                                                                                                                          ggCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAAACGGGTGTTCGAACACCAGCAAGAGGACATGGACGTCAAtgg           >  1:945139/1‑148 (MQ=255)
                                                                                                                                           gCCACCTGACCGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTgc          >  1:698528/1‑150 (MQ=255)
                                                                                                                                             cACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTgcgg        >  1:903288/1‑150 (MQ=255)
                                                                                                                                             cACCGGTCAGGTGGCCGGCGGGGCCAACCTGAGCGCCTTCACCACCGGCCGCAGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCCGCAAGAGGACATGGACGTCAATTgcgg        >  2:2096514/1‑150 (MQ=255)
                                                                                                                                              aCCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTgcggc       <  1:718621/150‑1 (MQ=255)
                                                                                                                                                   tCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGCGGCGGCAt  >  1:2281720/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACCACTATCCTCGAGAAATCCCTGGGCGCAGTGGCCAAGGCAGGCTCCAGCAACCTGGTCGATGTGTATCAGTACGCGCAAGCCGTGCGTGCCCAGGGCCTGGTGTTCATGGACACCCCCGGCTACGACCCGGTCTCGGCCACCGGTCAGGTGGCCGGCGGGGCCAACCTGATCGCCTTCACCACCGGCCGCGGCTCGGCCTACGGCTGTGCCCCGGCGCCATCGATCAAGCTGGCAACCAACAACCGGGTGTTCGAGCACCAGCAAGAGGACATGGACGTCAATTGCGGCGGCAT  >  NC_002947/3236046‑3236342

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