Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,989,188 G→C 10.7% G132A (GGC→GCC)  PP_1780 → mannosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,989,1880GC10.7% 105.9 / 3.3 37G132A (GGC→GCC) PP_1780mannosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (18/15);  new base C (2/2);  total (20/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 = 9.74e-01

TGCCGTTGAAAATGCCGACGCCGTTGCCAAGCTTAAATCGTTGGCCGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCTGCTACTGAGCTTTTTGGCGGAATACGAAGGCTTTGGTCGTAGC  >  NC_002947/1989039‑1989330
                                                                                                                                                     |                                                                                                                                               
tGCCGTTGAAAATGCCGACGCCGTTGCCAAGCTTAAATCGTTGGCCGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCcgg                                                                                                                                                 >  2:486219/1‑150 (MQ=255)
         aaaTGCCGACGCCGTTGCCAAGCTTAAATCGTTGGCCGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGgctgct                                                                                                                                        <  1:324655/150‑1 (MQ=255)
             gCCGAGGTCGTTTCCAAGCTTAAGTCGTTGGCCGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTACGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATa                                                                                                                                    <  2:331811/150‑1 (MQ=255)
                gACGCCGTTGCCAAGCTTAAATCGTTGGCCGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGgc                                                                                                                                            <  1:782339/139‑1 (MQ=255)
                gACGCCGTTGCCAAGCTTAAATCGTTGGCCGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGgc                                                                                                                                            >  2:782339/1‑139 (MQ=255)
                    ccGTTGCCAAGCTTAAATCGTTGGCCGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGaaa                                                                                                                            <  2:94608/150‑1 (MQ=255)
                            aaGCTTAAATCGTTGGCCGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTCGAGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAAc                                                                                                                    >  1:607411/1‑150 (MQ=255)
                            aaGCTTAAATCGTTGGCCGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAAc                                                                                                                    >  2:330970/1‑150 (MQ=255)
                                             cGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCAtggt                                                                                                   >  2:351418/1‑150 (MQ=255)
                                                      ccagGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGCCGGgcttga                                                                                                                                        >  2:708082/1‑102 (MQ=255)
                                                      ccagGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGCCGGgcttga                                                                                                                                        <  1:708082/105‑4 (MQ=255)
                                                              tGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCtttt                                                                                  >  2:581398/1‑150 (MQ=255)
                                                                         ttCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGgc                                                                                                                                            <  2:200033/82‑1 (MQ=255)
                                                                         ttCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGgc                                                                                                                                            >  1:200033/1‑82 (MQ=255)
                                                                             acaTTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACACCAAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGcc                                                                   >  1:451578/1‑150 (MQ=255)
                                                                                   ccGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGaaa                                                             <  1:287183/150‑1 (MQ=255)
                                                                                   ccGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTCTCCGCCACATGGCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTAACTCGACCCCAGCCATGaaa                                                             >  1:251048/1‑150 (MQ=255)
                                                                                     gAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAg                                                           >  1:68881/1‑150 (MQ=255)
                                                                                         aTGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTa                                                       <  1:330970/150‑1 (MQ=255)
                                                                                              ccGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGaaaa                                                            >  2:244768/1‑140 (MQ=255)
                                                                                              ccGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGaaaa                                                            <  1:244768/140‑1 (MQ=255)
                                                                                               cGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAAccc                                                                                                                         <  2:427392/78‑1 (MQ=255)
                                                                                               cGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAAccc                                                                                                                         >  1:427392/1‑78 (MQ=255)
                                                                                                gTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTCGAAACCCCCAACGCCCAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCaa                                                >  2:88537/1‑150 (MQ=255)
                                                                                                   cccGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCt                                             >  2:728502/1‑150 (MQ=255)
                                                                                                     cGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCTGc                                           <  2:716148/150‑1 (MQ=255)
                                                                                                        ttttGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGcc                                                         >  1:745287/1‑133 (MQ=255)
                                                                                                        ttttGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGcc                                                         <  2:745287/133‑1 (MQ=255)
                                                                                                                 aGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCTGCTACTGAGCtttt                               >  1:738211/1‑150 (MQ=255)
                                                                                                                     ggTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCTGCTACTGAGCTTTTTGGc                           >  2:195597/1‑150 (MQ=255)
                                                                                                                     ggTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCTGCTACTGAGCTTTTTGGc                           <  1:351418/150‑1 (MQ=255)
                                                                                                                         gggCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCTGCTACTGAGCTTTTTGGCGGaa                       <  1:728502/150‑1 (MQ=255)
                                                                                                                               gcgcTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCAATGAAAAGCCTATCCACAACCTGCTAGTGGGCTTTTTGGCGGAATACGaa                 <  2:607411/150‑1 (MQ=255)
                                                                                                                                      gagTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGc                                                                    <  2:510994/92‑1 (MQ=255)
                                                                                                                                      gagTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGc                                                                    >  1:510994/1‑92 (MQ=255)
                                                                                                                                          tcagcaGCCCGCCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCTGCTACTGAGCttttt                              >  1:28555/6‑126 (MQ=255)
                                                                                                                                          tcagcaGCCCGCCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCTGCTACTGAGCttttt                              <  2:28555/121‑1 (MQ=255)
                                                                                                                                              aaGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCTGCTACTGAGCTTTTTGGCGGAATACGAAGGCTTTGGTCGTAGc  >  1:633609/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TGCCGTTGAAAATGCCGACGCCGTTGCCAAGCTTAAATCGTTGGCCGACGCCAGCCAGGTGCTGGTGTCCGGGTTCCACATTGCCGAACATGTCCCGTTCCCGGTTTTGCAACAGCTGGTAGGGCAGGCGCTTAGAGTACTCAAGCCCGGCGGGCTGCTGATACTGGAAACCCCCAACGCCGAAAACCTCATGGTGGGGACCAACAGCTTTTACCTCGACCCCAGCCATGAAAAGCCTATCCCCAACCTGCTACTGAGCTTTTTGGCGGAATACGAAGGCTTTGGTCGTAGC  >  NC_002947/1989039‑1989330

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: