Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,858,702 C→G 8.3% intergenic (‑208/+17) PP_mr50 ← / ← ligA Bacteria_large_SRP ncRNA/DNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,858,7020CG 8.3% 150.9 / 3.1 48intergenic (‑208/+17)PP_mr50/ligABacteria_large_SRP ncRNA/DNA ligase
Reads supporting (aligned to +/- strand):  ref base C (25/19);  new base G (2/2);  total (27/21)
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.73e-01

ATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATGGGTCTTGC  >  NC_002947/4858557‑4858849
                                                                                                                                                 |                                                                                                                                                   
aTACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAgg                                                                                                                                                  >  7:231258/1‑149 (MQ=255)
 tACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCATGTTTTGAGAAGACAGTTGCCCCCCACAggg                                                                                                                                                 <  1:27739/149‑1 (MQ=255)
                         catcatTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTAt                                                                                                                         >  4:176991/1‑149 (MQ=255)
                             atTTCTCAATTACTTGTCAGGGTCGAGATGGCCTGGACGTGCCTCTTCGCGTGTAAACCCGCTCCTACTGTCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGcc                                                                                                                      <  5:200577/148‑1 (MQ=39)
                                 ctcAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAGCCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAgg                                                                                                                                                  >  6:297968/1‑116 (MQ=255)
                                    aaTGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCt                                                                                                                                   <  4:39142/128‑1 (MQ=255)
                                    aaTGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCt                                                                                                                                   >  3:39142/1‑128 (MQ=255)
                                    aaTGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCgg                                                                                                              >  8:271746/1‑149 (MQ=255)
                                         gttGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAg                                                                                                          <  6:202784/147‑1 (MQ=255)
                                                   tCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCAACAGGGGGGCACGGATGCCTAGAACGCTATGCCCTGTTCGGCCAGGAATGTGACa                                                                                                >  1:120365/1‑148 (MQ=255)
                                                            ccTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCt                                                                                       >  6:104814/1‑148 (MQ=255)
                                                             cTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAggggg                                                                                                                                               >  6:76957/1‑91 (MQ=255)
                                                             cTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTc                                                                                     >  7:111711/1‑149 (MQ=255)
                                                             cTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTc                                                                                     <  8:231258/149‑1 (MQ=255)
                                                                gACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTc                                                                                  >  3:69371/1‑149 (MQ=255)
                                                                            cgcgGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACaaa                                                                                              <  2:324156/125‑1 (MQ=255)
                                                                            cgcgGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACaaa                                                                                              >  1:324156/1‑125 (MQ=255)
                                                                            cgcgGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTca                                                                       <  3:176991/148‑1 (MQ=255)
                                                                             gcgGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTcaca                                                                     >  5:6872/1‑149 (MQ=255)
                                                                                   aaaCCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGcc                                                                                                                      >  3:232685/1‑94 (MQ=255)
                                                                                   aaaCCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGcc                                                                                                                      <  4:232685/94‑1 (MQ=255)
                                                                                      cccGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCg                                                            <  8:111711/149‑1 (MQ=255)
                                                                                       ccGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCg                                                                                    <  7:18626/124‑1 (MQ=255)
                                                                                       ccGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCg                                                                                    >  8:18626/1‑124 (MQ=255)
                                                                                       ccGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCgg                                                           >  2:147017/1‑149 (MQ=255)
                                                                                       ccGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCgg                                                           >  4:8008/1‑149 (MQ=255)
                                                                                       ccGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCgg                                                           >  1:221012/1‑149 (MQ=255)
                                                                                             cTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCt                                                     <  7:271746/149‑1 (MQ=255)
                                                                                                cGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTgg                                                  <  1:93212/149‑1 (MQ=255)
                                                                                                cGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTgg                                                  >  2:93212/1‑149 (MQ=255)
                                                                                                 gggCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGc                                                 <  6:6872/149‑1 (MQ=255)
                                                                                                        cgcgCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTc                                          <  2:221012/149‑1 (MQ=255)
                                                                                                                 ttttCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGcc                                 >  8:55714/1‑149 (MQ=255)
                                                                                                                     cAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAg                             >  6:98142/1‑149 (MQ=255)
                                                                                                                      aaGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                            >  8:35887/1‑149 (MQ=255)
                                                                                                                         ttttGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                            <  5:142377/146‑1 (MQ=255)
                                                                                                                         ttttGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                            >  6:142377/1‑146 (MQ=255)
                                                                                                                           ttGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCa                        <  7:55714/148‑1 (MQ=255)
                                                                                                                               gaagaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGccac                   >  3:130994/1‑149 (MQ=255)
                                                                                                                                aagaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                            <  7:184751/139‑1 (MQ=255)
                                                                                                                                aagaCAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                            >  8:184751/1‑139 (MQ=255)
                                                                                                                                     ccgtGCCCCCCTGTGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                            >  7:206281/4‑134 (MQ=255)
                                                                                                                                     ccgtGCCCCCCTGTGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAgg                            <  8:206281/131‑1 (MQ=255)
                                                                                                                                      cgtGCCCCCCTGTGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTgg                                                  <  2:216443/109‑1 (MQ=255)
                                                                                                                                      cgtGCCCCCCTGTGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTgg                                                  >  1:216443/3‑111 (MQ=255)
                                                                                                                                         gCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATgg         >  1:217780/1‑149 (MQ=255)
                                                                                                                                          ccccccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATgg         <  3:8008/148‑1 (MQ=255)
                                                                                                                                             cccACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGgccgcc                                                      <  2:29720/100‑1 (MQ=255)
                                                                                                                                                acaGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATGGGTCTTGc  >  2:196609/1‑149 (MQ=255)
                                                                                                                                                acaGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATGGGTCTTGc  <  4:69371/149‑1 (MQ=255)
                                                                                                                                                acaGGGGGGCACGGATGCCTAGAACGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATGGGTCTTGc  <  2:120365/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATACAACCACTGCGGCAGAAAAAGTCATCATTTCTCAATGACTTGTCAGGGTCGAGATGGCCTGGACGGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGACCGCGCAGGCTTTTCAAGTTTTGAGAAGACAGTTGCCCCCCACAGGGGGGCACGGATGCCTAGACCGCTATGCCCTGTTCGGCCAGGAATGTGACAAAGGCATCTTCGTCCAGCACCTTCACACCCAACTCGGCCGCCTTGGCCAGCTTCGAGCCTGCCCCAGGCCCAGCCACCACGCAATGGGTCTTGC  >  NC_002947/4858557‑4858849

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: