Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,037,644 C→T 10.8% C254C (TGC→TGT)  PP_0896 → carbon‑nitrogen hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,037,6440CT10.8% 92.8 / 3.6 37C254C (TGC→TGT) PP_0896carbon‑nitrogen hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base C (15/18);  new base T (2/2);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

CAGAGCCCGACGGTGGGGGTGGCGCCCTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCC  >  NC_002947/1037504‑1037792
                                                                                                                                             |                                                                                                                                                    
cAGAGCCCGACGGTGGGGGTGGCGCCCTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGt                                                                                                                                               <  4:92348/148‑1 (MQ=255)
   agCCCGACGGTGGGGGTGGCGCCCTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGGAATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGc                                                                                                                                            <  1:66865/149‑1 (MQ=255)
    gCCCGACGGTGGGGGTGGCGCCCTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGAGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTg                                                                                                                                          <  2:232147/149‑1 (MQ=255)
       cGACGGTGGGGGTGGCGCCCTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGgtg                                                                                                                                       >  4:68529/1‑149 (MQ=255)
          cGGTGGGGGTGGCGCCCTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGgtgtgt                                                                                                                                    <  3:241286/149‑1 (MQ=255)
             tGGGGGTGGCGCCCTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATAACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAg                                                                                                                                 >  1:93865/1‑149 (MQ=255)
                  gTGGCGCCCTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGgtgt                                                                                                                                      >  4:18494/1‑139 (MQ=255)
                  gTGGCGCCCTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGgtgt                                                                                                                                      <  3:18494/139‑1 (MQ=255)
                        cccTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTgga                                                                                                                      <  3:68529/149‑1 (MQ=255)
                                    gCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTTGTGTGTGAGGTAGAATTGCAGGAAGTgcggct                                                                                                          >  6:187941/1‑148 (MQ=255)
                                             gagaATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑AGGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCg                                                                                                 >  6:173014/1‑149 (MQ=255)
                                                    tCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAgg                                                                                           <  3:147637/148‑1 (MQ=255)
                                                      ggCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGGAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGc                                                                                         <  5:187941/148‑1 (MQ=255)
                                                       gCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAg                                                                                       <  8:94745/149‑1 (MQ=255)
                                                                  ttGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTg                                                                            >  5:80236/1‑149 (MQ=255)
                                                                    gTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGc                                                                          >  5:184604/1‑149 (MQ=255)
                                                                     ttCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCg                                                                         <  5:220058/149‑1 (MQ=255)
                                                                               cgGATAACGGGAGGCCGGGGG‑AGGGGGTGGTGGCGTTGAGTGTGCTGCTGAACCGGGAGTGCTGCGTGTGGCTGGTGTGTGAGGTAGATTTGGTGGAAGTGCGGCGGGTGAGGCGGGATGGGGAGGTGTTTACCTGGCGGGACTGGCCg                                                               <  2:93865/149‑1 (MQ=255)
                                                                                  aTCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGGTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGGGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAg                                                            <  6:79683/149‑1 (MQ=255)
                                                                                   tCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCgg                                                                        >  2:9273/1‑136 (MQ=255)
                                                                                   tCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCgg                                                                        <  1:9273/136‑1 (MQ=255)
                                                                                      cGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGt                                                                                                                >  3:221559/1‑93 (MQ=255)
                                                                                      cGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGt                                                                                                                <  4:221559/93‑1 (MQ=255)
                                                                                      cGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGt                                                        >  1:131593/1‑149 (MQ=255)
                                                                                             ccGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCga                                                    >  3:138555/1‑146 (MQ=255)
                                                                                                     aTGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCtg                                          <  2:131593/149‑1 (MQ=255)
                                                                                                            ggtggCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGAc                                   >  6:99662/1‑149 (MQ=255)
                                                                                                                  gTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTggc                             >  2:239310/1‑149 (MQ=255)
                                                                                                                      agtgagCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTggcggc                          <  1:239310/148‑1 (MQ=255)
                                                                                                                        tgagCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCctc                       >  7:24778/1‑149 (MQ=255)
                                                                                                                        tgagCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCct                        <  6:184604/148‑1 (MQ=255)
                                                                                                                                tGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCgg                >  7:99136/1‑148 (MQ=255)
                                                                                                                                cgagctcGGACTGTCGCCAGTGGCTGGTGTGGGAGGTAGATTTGGCGGACGGGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCggg               <  4:49143/143‑1 (MQ=255)
                                                                                                                                  aaCCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCgg                <  8:99136/146‑1 (MQ=255)
                                                                                                                                       ggACTGTAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGCGTCCTGTGATGACGGTGGCGGCCTCTTCGCGGGCagctcg        >  3:49143/1‑144 (MQ=255)
                                                                                                                                       ggACTGTAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCggg               <  4:45541/142‑1 (MQ=255)
                                                                                                                                       ggACTGTAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCggg               >  3:45541/1‑142 (MQ=255)
                                                                                                                                             cAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTccc  <  3:230119/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CAGAGCCCGACGGTGGGGGTGGCGCCCTGGTCGCCGGCATTGGACGAGAATATCGGCCGGGCAGGCTTGTTCGTACCGCCGGATCACGGGATGCCGGGGG‑ATGGGGTGGTGGCGTTGAGTGAGCAGCTGAACCCGGACTGCAGCCAGTGGCTGGTGTGTGAGGTAGATTTGGAGGAAGTGCGGCGGGTGAGGCGGGAAGGGCAGGTGTTTACCTGGCGGGACTGGCCGGAGCAGTTCGAGAGAGTCCTGTGATGACTGTGGCGGCCTCTTCGCGGGCAAGCCCGCTCCC  >  NC_002947/1037504‑1037792

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