Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,714,844 C→A 11.4% S101S (TCG→TCT)  metH ← cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,714,8440CA11.4% 101.4 / 2.9 35S101S (TCG→TCT) metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base C (16/15);  new base A (2/2);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGAGCAGGTTCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTAGGCCTTCTCGATGGCGGCGATCACGTCCGGGCGGCTGAGCAGCAAAAGGTCGTTGTTACCCT  >  NC_002947/2714713‑2714989
                                                                                                                                   |                                                                                                                                                 
aTCGAGCAGGTTCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTACATGCCGTAGTCg                                                                                                                                 >  1:251994/1‑150 (MQ=255)
   gAGCAGGTTCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGAGACTCCATGCCGTAGTCGGcc                                                                                                                              >  1:170207/1‑150 (MQ=255)
     gCAGGTTCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTg                                                                                                                            >  2:1126594/1‑150 (MQ=255)
     gCAGGTTCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGAGACTCCATGCCGTAGTCGGCCTg                                                                                                                            >  2:1106455/1‑150 (MQ=255)
       aGGTTCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTggg                                                                                                                          >  2:964511/1‑150 (MQ=255)
       aGGTTCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTggg                                                                                                                          >  1:1007270/1‑150 (MQ=255)
         gTTCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTc                                                                                                                                  <  2:247241/140‑1 (MQ=255)
         gTTCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTc                                                                                                                                  >  1:247241/1‑140 (MQ=255)
          ttCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGaaa                                                                                                                       <  1:943923/150‑1 (MQ=255)
               cTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGc                                                                                                                  <  1:1161204/150‑1 (MQ=255)
               cTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGAGACTCCATGCCGTAGTCGGCCTGGGAAATCTGc                                                                                                                  <  1:1106455/150‑1 (MQ=255)
                          gAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGt                                                                                                                                   <  2:604732/122‑1 (MQ=255)
                          gAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGt                                                                                                                                   >  1:604732/1‑122 (MQ=255)
                          gAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGaa                                                                                                       <  2:984898/150‑1 (MQ=255)
                                        cGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCgg                                                                                         >  1:1208056/1‑149 (MQ=255)
                                        cGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCga                                                                                         <  2:251994/150‑1 (MQ=255)
                                         gAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCgag                                                                                        >  2:705986/1‑150 (MQ=255)
                                            gcgcGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAat                                                                                     >  2:888064/1‑150 (MQ=255)
                                                            gTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATc                                                                                  <  2:1096728/137‑1 (MQ=255)
                                                            gTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATc                                                                                  >  1:1096728/1‑137 (MQ=255)
                                                            gTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCg                                                                                 <  1:334955/138‑1 (MQ=255)
                                                            gTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCg                                                                                 >  2:334955/1‑138 (MQ=255)
                                                                 cAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAgg                                                                                                     <  1:1152149/113‑1 (MQ=255)
                                                                 cAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAgg                                                                                                     >  2:1152149/1‑113 (MQ=255)
                                                                 cAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTa                                                                >  2:1124767/1‑150 (MQ=255)
                                                                          ggCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTAGGCCTtctc                                                       <  2:1007270/150‑1 (MQ=255)
                                                                                  cAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTAGGCCTTCTCGATggcgg                                               <  1:1126594/150‑1 (MQ=255)
                                                                                      cTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTAGGCCTTCTCGATGGCGGCGAt                                           <  1:1035406/150‑1 (MQ=255)
                                                                                       tGGCGAGCGATGCGCGCGCCTTAGACGGGCAACTCGTACCCCAGAGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGCCGTTGAAGGTGTTGCTTTCCAGGATGTCGGCAACGGCATCCAGGTAGGCCTTCTCGATGGCGGCGATc                                          <  2:170207/150‑1 (MQ=255)
                                                                                                               aCGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTAGGCCTTCTCGa                                                     >  2:145094/1‑115 (MQ=255)
                                                                                                               aCGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTAGGCCTTCTCGa                                                     <  1:145094/115‑1 (MQ=255)
                                                                                                                 gTTCAACTCGTAAACCAGCGACTCAATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTATCGAGAATATCGGCACCGGCATCCAGGTAGGCCTTCTCGATGGCGGCGATCACGTCCGGGCGGCTGAGCAGCAAAAg                >  1:886576/1‑150 (MQ=255)
                                                                                                                  ttCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTAGGCCTTCTCGATGGCGGCGATCACGTCCGGGCGGCTGAGCAGCAAAAgg               <  2:1208056/150‑1 (MQ=255)
                                                                                                                       cTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTAGGCCTTCTCGATGGCGGCGATCAAGTCCTGGCGGCTGAGCAGCAAAAGGTCgtt          <  1:1124767/150‑1 (MQ=255)
                                                                                                                               ccAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTAGGCCTTCTCGATGGCGGCGATCACGTCCGGGCGGCTGAGCAGCAAAAGGTCGTTGTTACCCt  >  2:115972/1‑150 (MQ=255)
                                                                                                                                   |                                                                                                                                                 
ATCGAGCAGGTTCGGCTGGTAGGGCCGAGCACGCCGGCAACGAAGCGCGGCTTGTCCGGCGTCTCCAGGGTCTTGGCATCGGCAACCTGGCGAGCGATGCGCGCGCCTTCGACGTTCAACTCGTAAACCAGCGACTCCATGCCGTAGTCGGCCTGGGAAATCTGCGTGGCGTTGAAGGTGTTGGTTTCGAGAATATCGGCACCGGCATCCAGGTAGGCCTTCTCGATGGCGGCGATCACGTCCGGGCGGCTGAGCAGCAAAAGGTCGTTGTTACCCT  >  NC_002947/2714713‑2714989

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