Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,120 A→C 16.7% intergenic (‑130/+24) PP_4560 ← / ← csbD ribonuclease BN/stress response protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,1200AC16.7% 65.6 / 11.7 36intergenic (‑130/+24)PP_4560/csbDribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base A (19/11);  new base C (3/3);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

CGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGTTT  >  NC_002947/5181993‑5182266
                                                                                                                               |                                                                                                                                                  
cGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaa                                                                                                                              <  1:592318/150‑1 (MQ=255)
cGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCAGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtcaa                                                                                                                              >  2:592318/1‑150 (MQ=255)
      cTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTGCACCCGAGAAGAGGCCGGCATCGCCtcaggtcaag                                                                                                                        >  1:552404/1‑150 (MQ=255)
              gcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGcc                                                                                                                                           <  2:726829/123‑1 (MQ=255)
              gcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGcc                                                                                                                                           >  1:726829/1‑123 (MQ=255)
               cTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGc                                                                                                                                            <  2:459368/121‑1 (MQ=255)
               cTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGc                                                                                                                                            >  1:459368/1‑121 (MQ=255)
                                  aGTCTAATAGACCGTGCCCGAGCTGGGCTGCGCCCTGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTcc                                                                                            <  1:196273/150‑1 (MQ=255)
                                        aTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGAc                                                                                      >  1:19279/1‑150 (MQ=255)
                                           gACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                              <  2:80272/91‑1 (MQ=37)
                                           gACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGCGAagag                                                                                                                                              >  1:80272/1‑91 (MQ=37)
                                                          gtgCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCttt                                                                    <  2:19279/150‑1 (MQ=255)
                                                             ctgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAg                                                                 >  1:665110/1‑150 (MQ=255)
                                                              tgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaaga                                                                                                                                               >  2:768648/1‑71 (MQ=21)
                                                              tgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaaga                                                                                                                                               <  1:768648/71‑1 (MQ=21)
                                                              tgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGAATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGc                                                                >  2:181650/1‑150 (MQ=255)
                                                                             ttGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTTTTGACCGCTTCCTTCACTTTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCgccttgcct                                                 >  1:93652/1‑150 (MQ=255)
                                                                             ttGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGGAAACGCTGCTTCACCGCATCATTAAAACTGCCCTTGACCGGCTGCGCCTCGCCTTTGAGCTCCCGCaccttgcct                                                 >  1:178613/1‑150 (MQ=255)
                                                                               gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTca                                                                                                   <  1:301421/98‑1 (MQ=255)
                                                                               gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTca                                                                                                   >  2:301421/1‑98 (MQ=255)
                                                                               gaagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCAACCCTCACATTTCCCTTGACCTGCTGCGCCCCGCCTTTGAACTCCTGCGCCTTGCCTTc                                               >  2:164138/1‑150 (MQ=255)
                                                                                aagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCg                                              >  2:439974/1‑150 (MQ=255)
                                                                                     tttGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCg                                         >  1:248551/1‑150 (MQ=255)
                                                                                          ttACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCtt                                                                                         <  2:287616/97‑1 (MQ=255)
                                                                                          ttACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCtt                                                                                         >  1:287616/1‑97 (MQ=255)
                                                                                              gAATAAGTGCCGGCTTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                                <  2:55678/150‑1 (MQ=255)
                                                                                              gAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                                >  2:385879/1‑150 (MQ=255)
                                                                                                aTAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCATCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGt                              <  2:552404/150‑1 (MQ=255)
                                                                                                        cGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTc                      >  2:51672/1‑150 (MQ=255)
                                                                                                         ggCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCa                     >  1:429449/1‑150 (MQ=255)
                                                                                                           ccTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCAcc                   >  2:438816/1‑150 (MQ=255)
                                                                                                           ccTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTGTACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCAcc                   >  1:376192/1‑150 (MQ=255)
                                                                                                            ctctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTc                                               >  2:582172/1‑121 (MQ=38)
                                                                                                            ctctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTc                                               <  1:582172/121‑1 (MQ=38)
                                                                                                               ttCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTc                          >  2:277866/1‑139 (MQ=255)
                                                                                                               ttCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTc                          <  1:277866/139‑1 (MQ=255)
                                                                                                                           cccGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTATTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGtt   <  1:439974/150‑1 (MQ=255)
                                                                                                                            ccGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAACTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGttt  <  1:305120/150‑1 (MQ=255)
                                                                                                                               |                                                                                                                                                  
CGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGTTT  >  NC_002947/5181993‑5182266

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: ‑