Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,536 C→A 14.4% intergenic (+52/+42) PP_4593 → / ← PP_4594 hypothetical protein/cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,212,5360CA14.4% 106.2 / 5.8 41intergenic (+52/+42)PP_4593/PP_4594hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base C (15/20);  new base A (2/4);  total (17/24)
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.92e-01

CGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCA  >  NC_002947/5212391‑5212684
                                                                                                                                                  |                                                                                                                                                    
cGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACcgtc                                                                                                                                                  <  1:176050/150‑1 (MQ=255)
    aCCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAg                                                                                                                                              <  1:312326/150‑1 (MQ=255)
         cgacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTTACTCCACCTGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCg                                                                                                                                         <  1:139055/150‑1 (MQ=255)
           acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGACGGT‑GTAGGAGCggg                                                                                                                                       <  2:686999/150‑1 (MQ=255)
                 ccGTCAACCGCCTGCGCAAAGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGT‑GTAGGAGCgggcttgcc                                                                                                                                 <  1:496974/150‑7 (MQ=37)
                   gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCcg                                                                                                                               >  2:881399/1‑150 (MQ=255)
                   gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCcg                                                                                                                               >  2:881075/1‑150 (MQ=255)
                        ccGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGaaa                                                                                                                          >  1:863008/1‑150 (MQ=255)
                              gcgcAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAACCCGCTCCTAC‑GACGGTGTAGGAGCGGGCTTGCCCGCGaagaggccg                                                                                                                    <  1:855940/150‑8 (MQ=25)
                              cagagcagAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGATAAAAC‑GCCGTCTGAGGCGCGGGTTTACCCGCGAAAGGGCCt                                                                                                                    >  2:820412/8‑150 (MQ=255)
                                    ggAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                                              <  1:829264/150‑1 (MQ=255)
                                            tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑CCCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacaca                                                                                                      <  2:993477/150‑1 (MQ=255)
                                            tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacaca                                                                                                      <  1:881399/150‑1 (MQ=255)
                                            tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacaca                                                                                                      <  1:881075/150‑1 (MQ=255)
                                                       gACAAGATCGAACGGATCATTCATCCGCCCAAATCGTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcgc                                                                                           >  1:366364/1‑150 (MQ=255)
                                                        aCAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                                          >  2:962402/1‑150 (MQ=255)
                                                         cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCtt                                                                                         <  2:464725/150‑1 (MQ=255)
                                                          aaGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCATCTTCGCGGGCAAACCCGCTCCTAC‑GACGGTGTAGGAGCGGGctt                                                                                                                                    >  1:779961/1‑106 (MQ=25)
                                                          aaGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCATCTTCGCGGGCAAACCCGCTCCTAC‑GACGGTGTAGGAGCGGGctt                                                                                                                                    <  2:779961/106‑1 (MQ=25)
                                                            gATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAg                                                                                              >  1:535561/1‑142 (MQ=255)
                                                            gATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAg                                                                                              <  2:535561/142‑1 (MQ=255)
                                                                  aCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTCTTGACACCTACGGCGTCTTCGCTGGCAAGCCCGAAACAAC‑TCCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGccc                                                                                <  1:520085/150‑1 (MQ=37)
                                                                       tcattcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGAGCTTTCAGCGCCACTTcc                                                                           >  2:793951/1‑150 (MQ=255)
                                                                          ttcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAgg                                                                        >  2:156973/1‑150 (MQ=255)
                                                                               ccGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                                   <  1:156973/150‑1 (MQ=255)
                                                                                     aaaTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAg                                                             <  2:863008/150‑1 (MQ=255)
                                                                                          cTGACAGGCATCTATTGAATCCACCGGCCTCTTCCAGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCAc                                                        <  2:366364/150‑1 (MQ=255)
                                                                                            gACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACTGGCGAAAGGGCCGGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg                                                      >  1:411485/1‑150 (MQ=255)
                                                                                                 gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAgg                                                                                                                        >  1:936034/1‑79 (MQ=255)
                                                                                                 gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAgg                                                                                                                        <  2:936034/79‑1 (MQ=255)
                                                                                                         tGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATcc                                         <  1:795017/150‑1 (MQ=255)
                                                                                                         tGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATcc                                         <  1:793951/150‑1 (MQ=255)
                                                                                                               cACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTGTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCg                                   >  2:739139/1‑150 (MQ=255)
                                                                                                                                    gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATGCTCTGTCCCTACCGACAACCGCAGCAGGTTCTCGc              >  2:88328/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCTACGACGGT‑GTAGGAGCGGGCTTACCCGTGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc                                                                                             <  1:1075478/69‑1 (MQ=17)
                                                                                                                                      ccGCTCCTACGACGGT‑GTAGGAGCGGGCTTACCCGTGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc                                                                                             >  2:1075478/1‑69 (MQ=17)
                                                                                                                                      ccGCTCCTAC‑GACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACg                                                       <  2:437359/107‑1 (MQ=38)
                                                                                                                                      ccGCTCCTAC‑GACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACg                                                       >  1:437359/1‑107 (MQ=38)
                                                                                                                                      ccGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCGCTAACCCAACCGACAACCGCAGCAGCTTCCCGCTg            >  2:456427/1‑150 (MQ=39)
                                                                                                                                        gCTCCTAC‑GACGGTGTAGGAGCGGGCTTGCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgcagca                       <  1:375864/137‑1 (MQ=38)
                                                                                                                                             tAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCa     <  1:962402/150‑1 (MQ=255)
                                                                                                                                                 acCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCa  >  1:227587/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCATCA  >  NC_002947/5212391‑5212684

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