Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,991,165 C→G 8.9% A805A (GCC→GCG) ‡ PP_1780 → mannosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,991,1650CG 8.9% 131.1 / 4.5 45A805A (GCC→GCG) ‡PP_1780mannosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base G (2/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTTCG  >  NC_002947/1991018‑1991308
                                                                                                                                                   |                                                                                                                                               
gCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGcc                                                                                                                                                 >  6:186856/1‑148 (MQ=255)
gCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGcc                                                                                                                                                 >  8:254948/1‑148 (MQ=255)
gCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGAGCAGGCGGTTTTAGCCAGTGcc                                                                                                                                                 >  4:136199/1‑148 (MQ=255)
 cATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCaa                                                                                                                                               <  3:23497/149‑1 (MQ=255)
      cTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGc                                                                                                                                          >  4:144059/1‑149 (MQ=255)
           cccTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCaaa                                                                                                                                     <  6:124497/149‑1 (MQ=255)
              tCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAacac                                                                                                                                  >  7:195681/1‑149 (MQ=255)
              tCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAacac                                                                                                                                  >  6:377963/1‑149 (MQ=255)
                  gATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGaaa                                                                                                                              <  7:254948/149‑1 (MQ=255)
                      ggCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGcc                                                                                                                          <  5:186856/149‑1 (MQ=255)
                       gCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCcg                                                                                                                         >  1:106668/1‑149 (MQ=255)
                               cgcgGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGgcc                                                                                                                   <  5:157698/147‑1 (MQ=255)
                                gcgGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGgccgcc                                                                                                                <  2:331210/149‑1 (MQ=255)
                                            tcgAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTGAGCCAGTTGGCACTGGCTAAAaccgcct                                                                                                                               <  3:174375/122‑6 (MQ=38)
                                            tcgAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTGAGCCAGTTGGCACTGGCTAAAaccgcct                                                                                                                               >  4:174375/1‑117 (MQ=38)
                                             cgAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACATGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGAcc                                                                                                    >  1:272568/1‑148 (MQ=255)
                                              gAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCa                                                                                                   >  6:21664/1‑148 (MQ=255)
                                                          gCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTggtggtggt                                                                                      >  6:103063/1‑149 (MQ=255)
                                                          gCACTCGAGCAGCCCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGTGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTggtggtggt                                                                                      >  3:373306/1‑149 (MQ=255)
                                                             cTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAAc                                                                                                >  6:147300/1‑136 (MQ=255)
                                                             cTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAAc                                                                                                <  5:147300/136‑1 (MQ=255)
                                                                gAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTtctc                                                                                 <  2:106668/148‑1 (MQ=255)
                                                                   cagcGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTg                                                                             <  2:131971/149‑1 (MQ=255)
                                                                   cagcGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTg                                                                             >  5:282191/1‑149 (MQ=255)
                                                                    agcGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTgg                                                                            <  8:195681/149‑1 (MQ=255)
                                                                               cAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATg                                                                 >  2:317540/1‑149 (MQ=255)
                                                                                aaGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGc                                                                >  7:317206/1‑149 (MQ=255)
                                                                                     aTTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCAt                                                           >  2:333904/1‑149 (MQ=255)
                                                                                         aaCTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCt                                                       <  6:282191/149‑1 (MQ=255)
                                                                                           cTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTgg                                                                            <  6:151100/126‑1 (MQ=255)
                                                                                           cTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGa                                                     <  1:333904/149‑1 (MQ=255)
                                                                                                   cagcaTCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTc                                             <  1:342328/149‑1 (MQ=255)
                                                                                                     gcaTCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAg                                           <  8:317206/149‑1 (MQ=255)
                                                                                                     gcaTCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAg                                           <  5:21664/149‑1 (MQ=255)
                                                                                                        tCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCa                                        >  3:183525/1‑149 (MQ=255)
                                                                                                           acacTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTg                                     >  5:137505/1‑149 (MQ=255)
                                                                                                               ttGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAAAACGAAATGCAGGAGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACAATCAGGCATTGCt                                   >  1:46334/1‑147 (MQ=255)
                                                                                                                   gAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGa                             <  5:103063/149‑1 (MQ=255)
                                                                                                                    aGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAg                            >  1:7532/1‑149 (MQ=255)
                                                                                                                         agTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCaa                       <  1:317540/149‑1 (MQ=255)
                                                                                                                                aGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGAc                <  3:144059/149‑1 (MQ=255)
                                                                                                                                 ggCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCa                                 <  5:254593/131‑1 (MQ=255)
                                                                                                                                 ggCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCa                                 >  6:254593/1‑131 (MQ=255)
                                                                                                                                    tgTTTGAGCCAGTTGGCACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATaca                                                <  4:354136/112‑1 (MQ=255)
                                                                                                                                    tgTTTGAGCCAGTTGGCACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATaca                                                >  3:354136/2‑113 (MQ=255)
                                                                                                                                    ggTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGGTCAGTGAGTTCAAGCAGGACGAt             >  2:118406/1‑148 (MQ=255)
                                                                                                                                          aGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTacgacg      >  1:157925/1‑149 (MQ=255)
                                                                                                                                               gTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTTCg  >  1:378738/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCGCCGGCGCCAAGGCATTGAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCGGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTTCG  >  NC_002947/1991018‑1991308

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