Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,624,187 G→C 14.3% intergenic (‑92/+81) PP_3194 ← / ← PP_3195 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,624,1870GC14.3% 92.5 / 8.7 35intergenic (‑92/+81)PP_3194/PP_3195hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/17);  new base C (2/3);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GGCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGC  >  NC_002947/3624047‑3624332
                                                                                                                                            |                                                                                                                                                 
ggCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATAc                                                                                                                                          <  2:7431/150‑1 (MQ=255)
        ccGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGct                                                                                                                                  >  2:499682/1‑150 (MQ=255)
                             tgcaggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAAc                                                                                                             <  1:174918/150‑1 (MQ=255)
                                       gggCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATc                                                                                                   >  2:232005/1‑150 (MQ=255)
                                         gCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACTGCTCCAATCgg                                                                                                 <  2:458403/150‑1 (MQ=255)
                                           aCTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATAc                                                                                                                                          >  1:354079/1‑107 (MQ=255)
                                           aCTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATAc                                                                                                                                          <  2:354079/107‑1 (MQ=255)
                                                 tGGTGATCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCGATTGGcgccgttc                                                                                         <  1:165686/150‑4 (MQ=255)
                                                    tGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAg                                                                                      >  2:479301/1‑150 (MQ=255)
                                                        ctctctCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAg                                                                                  <  1:356381/150‑1 (MQ=255)
                                                             tcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAg                                                                                  >  1:109230/1‑145 (MQ=255)
                                                             tcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAg                                                                                  <  2:109230/145‑1 (MQ=255)
                                                             tcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCg                                                                             >  1:603437/1‑150 (MQ=255)
                                                             tcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTAGTGATGGCCCTGTGTCGATGGCAGCGATACCTGGCGCTCTCAGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCg                                                                             <  2:50305/150‑1 (MQ=255)
                                                             tcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGACCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCAGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCg                                                                             <  2:50798/150‑1 (MQ=255)
                                                                tCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCt                                                                                                                                        <  2:322112/88‑1 (MQ=255)
                                                                tCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCt                                                                                                                                        >  1:322112/1‑88 (MQ=255)
                                                                  ttGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAgc                                                                        <  2:337849/150‑1 (MQ=255)
                                                                        gtggtgCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGgc                                                                                                >  1:56286/1‑120 (MQ=255)
                                                                        gtggtgCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGgc                                                                                                <  2:56286/120‑1 (MQ=255)
                                                                             gCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGc                                                             >  1:208278/1‑150 (MQ=255)
                                                                                      ggctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGgcc                                                    <  1:479301/150‑1 (MQ=255)
                                                                                             gCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAgg                                             <  2:400840/150‑1 (MQ=255)
                                                                                                            cagcagATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACga                                                                                                                 <  1:247884/67‑1 (MQ=255)
                                                                                                            cagcagATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACga                                                                                                                 >  2:247884/1‑67 (MQ=255)
                                                                                                                    gCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGAt                      <  1:232005/150‑1 (MQ=255)
                                                                                                                        ttGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAg                  <  2:208278/150‑1 (MQ=255)
                                                                                                                            tGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGcc              >  2:33268/1‑150 (MQ=255)
                                                                                                                              gCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAgg                                                              <  2:422286/100‑1 (MQ=255)
                                                                                                                              gCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAgg                                                              >  1:422286/1‑100 (MQ=255)
                                                                                                                                   gTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGACGGACTCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCAtt       <  1:499682/150‑1 (MQ=255)
                                                                                                                                        gATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGc  >  2:477534/1‑150 (MQ=255)
                                                                                                                                        gATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGc  >  2:544723/1‑150 (MQ=255)
                                                                                                                                        gATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCgg           <  1:524959/141‑1 (MQ=255)
                                                                                                                                        gATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCgg           >  2:524959/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GGCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGC  >  NC_002947/3624047‑3624332

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