Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,172,489 A→G 10.5% T240A (ACG→GCG)  PP_2784 → short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,172,4890AG10.5% 97.6 / 2.2 38T240A (ACG→GCG) PP_2784short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (21/13);  new base G (2/2);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GGCGGCTCAAGCGCGCATCGCCGGGATCACCCATGAGCAGATCCTGCGCAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAGCGGCGATCGCCGAGTGGTTGCCGGAGCATCATGAAGCGG  >  NC_002947/3172344‑3172638
                                                                                                                                                 |                                                                                                                                                     
ggcggcTCAAGCGCGCATCGCCGGGATCACCCATGAGCAGATCCTGCGCAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGgc                                                                                                                                                   <  1:737289/150‑1 (MQ=255)
   ggcTCAAGCGCGCATCGCCGGGATCACCCATGAGCAGATCCTGCGCAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCa                                                                                                                                                >  1:1130538/1‑150 (MQ=255)
            gcgcATCGCCGGGATCACCCATGAGCAGATCCTGCGCAATGCCAGGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCtgat                                                                                                                                       <  2:870113/150‑1 (MQ=255)
            gcgcATCGCCGGGATCACCCATGAGCAGATCCTGCGCAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCtgat                                                                                                                                       <  1:273671/150‑1 (MQ=255)
                aTCGCCGGGATCACCCATGAGCAGATCCTGCGCAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATc                                                                                                                                   >  2:1059851/1‑150 (MQ=255)
                             cccATGAGCAGATCCTGCGCAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCg                                                                                                                      <  2:1130538/150‑1 (MQ=255)
                              ccATGAGCAGATCCTGCGCAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCgg                                                                                                                     >  2:389564/1‑150 (MQ=255)
                                aTGAGCAGATCCTGCGCAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGt                                                                                                                   >  2:780080/1‑150 (MQ=255)
                                  gAGCAGATCCTGCGCAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTg                                                                                                                 <  1:127540/150‑1 (MQ=255)
                                               gcAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCAt                                                                                                    >  1:119028/1‑150 (MQ=255)
                                                           agaagTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTggg                                                                                        <  1:457932/150‑1 (MQ=255)
                                                                 tGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGac                                                                                  <  1:278902/150‑2 (MQ=255)
                                                                      cTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCGCGTCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGAtt                                                                                     <  1:25512/142‑1 (MQ=255)
                                                                      cTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCGCGTCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGAtt                                                                                     >  2:25512/1‑142 (MQ=255)
                                                                               ttcattcaGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGGTGGGGGATGGACTGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTgg                                                                    >  1:209211/1‑150 (MQ=255)
                                                                               ttcattcaGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTgg                                                                    >  2:29818/1‑150 (MQ=255)
                                                                                 cattcaGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCAt                                                                                                    <  2:849165/116‑1 (MQ=255)
                                                                                 cattcaGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCAt                                                                                                    >  1:849165/1‑116 (MQ=255)
                                                                                    tcaGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCcgcg                                                               >  1:820687/1‑150 (MQ=255)
                                                                                        ggCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAAGCGAGTGGTTTCACGGCGCAAAGCCTGATGATCGACGGCGATGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGCGTGAACCTCGGCCTTTTGGCCGCGCTgg                                                           >  2:168037/1‑150 (MQ=255)
                                                                                        ggCGACGAAGTGGCGGCCATGGCTGGTTATCTGGCGAGCACGAAGGCGAGGGGTGTCAGGGCGCAAAGCCTGATGGTCTACGGCGGTGCGGGGTTGGGAATGTACCGATGGCTCGAGTGGGATTGAACCGCGGACTTTTGGCTGCGCTgg                                                           >  1:579124/1‑150 (MQ=255)
                                                                                                         cATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCa                                          >  1:23130/1‑150 (MQ=255)
                                                                                                             gTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAgcgg                                      <  1:780080/150‑1 (MQ=255)
                                                                                                             gTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAgcgg                                      <  1:430280/150‑1 (MQ=255)
                                                                                                             gTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGCGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAgcgg                                      >  1:974379/1‑150 (MQ=255)
                                                                                                              ttCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAgcggc                                     >  2:1035650/1‑150 (MQ=255)
                                                                                                                    aTCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAGCGGCGATCGc                               >  1:739051/1‑150 (MQ=255)
                                                                                                                        ggCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAAt                                                                                                            <  1:1004677/69‑1 (MQ=255)
                                                                                                                        ggCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAAt                                                                                                            >  2:1004677/1‑69 (MQ=255)
                                                                                                                               ccGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAAcc                                                                                >  2:163983/1‑90 (MQ=255)
                                                                                                                               ccGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAAcc                                                                                <  1:163983/90‑1 (MQ=255)
                                                                                                                                cGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAGCGGCGATCGCCGAGTGGTTGcc                   >  2:595998/1‑150 (MQ=255)
                                                                                                                                cGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAGCGGCGATCGCCGAGTGGTTGcc                   >  2:759977/1‑150 (MQ=255)
                                                                                                                                            gtgtCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCa                                          <  1:539914/115‑1 (MQ=255)
                                                                                                                                            gtgtCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCa                                          >  2:539914/1‑115 (MQ=255)
                                                                                                                                            gtgtCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAGCGGCGATCGCCGAGTGGTTGCCGGAGCATCATGa       >  1:719609/1‑150 (MQ=255)
                                                                                                                                             tgtCGCGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGcc                                                                  <  2:395535/90‑1 (MQ=255)
                                                                                                                                             tgtCGCGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGcc                                                                  >  1:395535/1‑90 (MQ=255)
                                                                                                                                                 aCGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCATTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAGCGGCGATCGCCGAGTGGTTGCCGGAGCATCATGAAGCgg  <  1:389564/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GGCGGCTCAAGCGCGCATCGCCGGGATCACCCATGAGCAGATCCTGCGCAATGCCAAGAAGAAGTTGCCCCTGCGGCGCTTCATTCAGGGCGACGAAGTGGCGGCCATGGTTCGTTATCTGGCGAGCCCGCAGGCGAGTGGTGTCACGGCGCAAAGCCTGATGATCGACGGCGGTGCCGGGTTGGGAATGTAGCCATGGCTCGACTGGGATTGAACCGCGGCCTTTTGGCCGCGCTGGTCAGCCTCGGCCCGCCAGCGGCGATCGCCGAGTGGTTGCCGGAGCATCATGAAGCGG  >  NC_002947/3172344‑3172638

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