Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,071,213 C→G 12.5% intergenic (+165/+62) PP_1846 → / ← PP_1847 group II intron‑encoding maturase/TonB‑dependent ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,071,2130CG12.5% 92.5 / 3.6 38intergenic (+165/+62)PP_1846/PP_1847group II intron‑encoding maturase/TonB‑dependent ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (15/18);  new base G (3/2);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCCA  >  NC_002947/2071067‑2071354
                                                                                                                                                   |                                                                                                                                             
gTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCgt                                                                                                                                             <  4:133696/149‑1 (MQ=255)
    gCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGccc                                                                                                                                               <  8:151228/143‑1 (MQ=255)
    gCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGccc                                                                                                                                               >  7:151228/1‑143 (MQ=255)
        ggtggtGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGG                                                                                                                                     <  7:105898‑M1/149‑1 (MQ=255)
         gcggtgtgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGG                                                                                                                                    <  6:252325‑M1/147‑1 (MQ=255)
               tgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCC                                                                                                                              <  7:91752‑M1/149‑1 (MQ=255)
                   agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCT                                                                                                                          <  3:4538‑M1/149‑1 (MQ=255)
                    ggACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTC                                                                                                                         >  7:214501‑M1/1‑149 (MQ=255)
                        ggcggGTGTGAATCCGCCTCCTACTCGATTTTCTAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCA                                                                                                                      >  1:173088‑M1/1‑149 (MQ=255)
                                     ccGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAgg                                                                                                        <  5:257440/149‑1 (MQ=255)
                                        cctcctACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAg                                                                                                         <  6:203246/145‑1 (MQ=255)
                                        cctcctACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAg                                                                                                         >  5:203246/1‑145 (MQ=255)
                                               cTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATAGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGc                                                                                              >  2:26288/1‑149 (MQ=255)
                                                tCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAgg                                                                                                        >  1:187702/1‑138 (MQ=255)
                                                  gATTTTCAGGGCCATGGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCg                                                                                             <  7:252120/147‑1 (MQ=255)
                                                    ttttCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGccc                                                                                         >  5:118189/1‑149 (MQ=255)
                                                             gCCATTGGGGCTGCTGTGCAGCCCATCGCGCCACAAGCGGGCTTGCCCGGGAAGCCACCACCGTGGTGGATGGCACGGGCTTCGCCGGTGTGCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATc                                                                                 <  3:91895/149‑1 (MQ=37)
                                                                         gctgTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATc                                                                     >  6:103060/1‑149 (MQ=255)
                                                                          ctgTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCTCAGAGCGGCCCCGGCAATCTAGAAGCTGATc                                                                     <  5:103060/148‑1 (MQ=255)
                                                                               cAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagcc                                                               <  5:99549/149‑1 (MQ=255)
                                                                               cAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagc                                                                <  2:173088/148‑1 (MQ=255)
                                                                                         ccaccccaGCGGGCTTGCCCGCTAACCAAACCCCGGGGGGCTTGGCAGGGGCTGGGCCGGTGTGCGCGGCAACGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAg                                                     <  6:5847/142‑1 (MQ=255)
                                                                                         ccaaccaAGCGGCCTTGCCCCCCAAGCCCACACCGGGGTGGATGGCAGGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCCCAGGCCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAg                                                     <  4:175520/144‑1 (MQ=21)
                                                                                            cacaAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGc                                                                          >  6:252483/1‑125 (MQ=255)
                                                                                            cacaAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGc                                                                          <  5:252483/125‑1 (MQ=255)
                                                                                               aaGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTAc                                               >  5:239017/1‑149 (MQ=255)
                                                                                                 gCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTAcgc                                             >  5:67954/1‑149 (MQ=255)
                                                                                                      cTTGCCCGCGAAGCAAGCACTGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGgtaaggtgctgcggtttagactttgtggtacagctggctgccttcctggcggaaccgctcggcctgctcacgca                                        >  6:53407‑M1/1‑75 (MQ=255)
                                                                                                       ttGCCCGCGAAGCAAGCACTGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGgtgaggtgctgcggtttagactttgtggtacagctggctgccttcctggcggaaccgctcggcctgctcacgcat                                       >  8:171334‑M1/1‑74 (MQ=255)
                                                                                                          cccGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGc                                                                          <  4:115453/111‑1 (MQ=255)
                                                                                                          cccGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGc                                                                          >  3:115453/1‑111 (MQ=255)
                                                                                                             gcgAAGCAAGCACTGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGgtaaggtgctgcggtttagactttgtggtacagctggctgccttcctggcggaaccgctcggcctgctcacgcatgccttg                                 >  4:222853‑M1/1‑68 (MQ=255)
                                                                                                                 aGCCAACCCCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACa                             <  6:67954/149‑1 (MQ=255)
                                                                                                                   ccAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAgg                           >  6:54478/1‑149 (MQ=255)
                                                                                                                          cGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTAcgc                                             >  4:218280/1‑124 (MQ=255)
                                                                                                                          cGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTAcgc                                             <  3:218280/124‑1 (MQ=255)
                                                                                                                            tggtggATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagccgaaa                                                           >  4:100647/1‑108 (MQ=255)
                                                                                                                            tggtggATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagccgaaa                                                           <  3:100647/108‑1 (MQ=255)
                                                                                                                                       aCGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTACTTATCCCCCACGAAAta       >  4:4341/1‑148 (MQ=255)
                                                                                                                                            cTTCGCCCGTGTTCGGGGGCAAGCCCGCTCCCACAGGTCGCTAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCca  <  3:4341/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGTAATCCCCCA  >  NC_002947/2071067‑2071354

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