Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,020,315 G→C 12.9% W270C (TGG→TGC PP_2638 → cellulose synthase operon protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,020,3150GC12.9% 78.5 / 6.2 31W270C (TGG→TGCPP_2638cellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base C (2/2);  total (14/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 = 5.42e-01

TCCGGGCAAGCCGTGAGTGCCACCAGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCAGTACCCGGGCGATGCCAGCCTGCATGGGGCCCTGGGTTATGCGTTGATGCGCCAGGGCC  >  NC_002947/3020181‑3020457
                                                                                                                                      |                                                                                                                                              
tCCGGGCAAGCCGTGAGTGCCACCAGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACgg                                                                                                                                 >  1:53299/1‑150 (MQ=255)
 ccGGGCAAGCCGTTACGGCCACCAGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCCGCATGGCAGGCCGGTCAACcgg                                                                                                                                <  1:212339/150‑1 (MQ=255)
     gCAAGCCGTGAGTGCCACCAGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACCTGCCGGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTaa                                                                                                                            >  1:412606/1‑150 (MQ=255)
       aaGCCGTGAGTGCCACCAGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAgg                                                                                                                          <  1:121375/150‑1 (MQ=255)
           cGTGATTGCCACCAGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCAtt                                                                                                                      <  1:501370/150‑1 (MQ=255)
              gAGTGCCACCAGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCt                                                                                                                   <  1:371894/150‑1 (MQ=255)
                     accaGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAg                                                                                                            >  1:351984/1‑150 (MQ=255)
                       caGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACggg                                                                                                                                <  2:41237/128‑1 (MQ=255)
                       caGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACggg                                                                                                                                >  1:41237/1‑128 (MQ=255)
                           gcgGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCCTGCCATGCCgggtc                                                                                                                                     >  2:434743/1‑116 (MQ=255)
                           gcgGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCCTGCCATGCCgggtc                                                                                                                                     >  2:434712/1‑116 (MQ=255)
                           gcgGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCCTGCCATGCCgggtc                                                                                                                                     <  1:434712/119‑4 (MQ=255)
                           gcgGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCCTGCCATGCCgggtc                                                                                                                                     <  1:434743/119‑4 (MQ=255)
                                        aaGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGa                                                                                              <  2:47240/145‑1 (MQ=255)
                                        aaGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGa                                                                                              >  1:47240/1‑145 (MQ=255)
                                              ttcttcAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGc                                                                                   >  2:48747/1‑150 (MQ=255)
                                               tcttcAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCa                                                                                  <  1:134757/150‑1 (MQ=255)
                                                    aGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTgc                                                                             >  1:249132/1‑150 (MQ=255)
                                                      cGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTgcgc                                                                           >  1:519342/1‑150 (MQ=255)
                                                            ccGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCa                                                                     >  1:69875/1‑150 (MQ=255)
                                                                    gCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCa                                                             <  1:309217/150‑1 (MQ=255)
                                                                    gCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCa                                                             >  1:289355/1‑150 (MQ=255)
                                                                            tGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCAGTACCCgg                                                     <  1:239766/150‑1 (MQ=255)
                                                                                 caGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCAGTACCCGGGCGAt                                                <  2:289355/150‑1 (MQ=255)
                                                                                            gaCCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAAAGGGGTAAGGCATTGCAAGCAAAGGGCTTTTATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCAGTACCCGGGCGATGCCAGCCTGCa                                     >  1:380989/1‑150 (MQ=255)
                                                                                               cTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCAGTACCCGGGCGATGCCAGCCTGCATgg                                  <  2:519342/150‑1 (MQ=255)
                                                                                                   cAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCAGTACCCGGGCGATGCCAGCCTGCATGGGGcc                              >  1:219764/1‑150 (MQ=255)
                                                                                                       agcagcGCAAGCTGCTCGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAACCGGAAACGCAGCTGCGCCGCGCACTGCGCCAGTACCCGGGCGATGCCAGCCTGCATGGGGCCCTgg                          <  1:230801/150‑1 (MQ=255)
                                                                                                                 gctgctGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAATCGCAGCTGCGCCGCGCACTGCGCCAGTACCCGGGCGATGCCAGCCTGCATGGGGCCCTGGGTTATGCGtt                <  2:292445/150‑1 (MQ=255)
                                                                                                                        gCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCAGTACCCGGGCGATGCCAGCCTGCATGGGGCCCTGGGTTATGCGTTGATgcgc         <  2:249132/150‑1 (MQ=255)
                                                                                                                               cGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCAGTACCCGGGCGATGCCAGCCTGCATGGGGCCCTGGGTTATGCGTTGATGCGCCAGGgcc  <  1:203039/150‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                              
TCCGGGCAAGCCGTGAGTGCCACCAGCGCGGCAGCCTGGCAAGCCTTTCTTCAGCGTTACCCGGCTTCGCCGTTGCTGGCACAGGCCAGCGAGACCTTGCAACAGCAGCGCAAGCTGCTGGCCGACCCGGCATGGCAGGCCGGTCAACGGGGTAAGGCATTGCTCGACAAGGGCCTTAATGCCGAAGCGGAAACGCAGCTGCGCCGCGCACTGCGCCAGTACCCGGGCGATGCCAGCCTGCATGGGGCCCTGGGTTATGCGTTGATGCGCCAGGGCC  >  NC_002947/3020181‑3020457

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: