Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,777,525 C→G 8.9% intergenic (+259/+33) PP_5476 → / ← PP_1585 hypothetical protein/antidote protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,777,5250CG 8.9% 130.0 / 4.5 45intergenic (+259/+33)PP_5476/PP_1585hypothetical protein/antidote protein
Reads supporting (aligned to +/- strand):  ref base C (13/28);  new base G (2/2);  total (15/30)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

TTTAGTCGGTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCAG  >  NC_002947/1777378‑1777671
                                                                                                                                                   |                                                                                                                                                    
tttAGTCGGTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGccc                                                                                                                                                     >  1:189016/1‑149 (MQ=255)
    gTCGGTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTgg                                                                                                                                                 >  2:144527/1‑149 (MQ=255)
     tCGGTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGgc                                                                                                                                                >  6:96806/1‑149 (MQ=255)
       ggTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGgcg                                                                                                                                               <  5:107091/148‑1 (MQ=255)
       ggTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGgcg                                                                                                                                               <  1:193277/148‑1 (MQ=255)
           tAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGGGCGTGGCGccgc                                                                                                                                           <  7:184140/148‑3 (MQ=255)
             aCTTGTCGCCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGGTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGGCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGGGCGTGGCGccgcga                                                                                                                                         <  5:22930/148‑5 (MQ=255)
             aCTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCTGCGCCACGCCCGTGGCGCCCCTg                                                                                                                                         <  4:75202/148‑1 (MQ=255)
                  gCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTAtt                                                                                                                                   <  2:28010/149‑1 (MQ=255)
                  gCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTGGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTAt                                                                                                                                    <  2:151332/148‑1 (MQ=255)
                       ccAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGac                                                                                                                                >  6:47850/1‑146 (MQ=255)
                        cAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAg                                                                                                                             <  3:112348/149‑1 (MQ=255)
                        cAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAg                                                                                                                             <  8:222295/149‑1 (MQ=255)
                         aGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAg                                                                                                                             >  4:81076/1‑148 (MQ=255)
                           gCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGAt                                                                                                                          <  6:6773/149‑1 (MQ=255)
                                      gTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCg                                                                                                               <  3:81076/149‑1 (MQ=255)
                                         gagCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGAc                                                                                                             <  3:110797/148‑1 (MQ=255)
                                               gtAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGAc                                                                                                             <  3:225941/142‑1 (MQ=255)
                                               gtAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGAc                                                                                                             >  4:225941/1‑142 (MQ=255)
                                                  ggTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTg                                                                                                    <  6:98648/148‑1 (MQ=255)
                                                   gTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGAt                                                                                                                          <  1:242726/125‑1 (MQ=255)
                                                   gTCATTGTGTACATGGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATCTTTGTCGCGGGTTTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGga                                                                                                  <  1:222250/149‑1 (MQ=255)
                                                    taaTTGTGTACATTGGGGTTGCAATTAGCCCCAATCTTTCCGTTTTTCGATCATTGTCGCGGTTTTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGga                                                                                                  <  1:222252/146‑1 (MQ=255)
                                                      aTTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTg                                                                                                                                                  >  7:69843/1‑98 (MQ=255)
                                                      aTTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTg                                                                                                                                                  <  8:69843/98‑1 (MQ=255)
                                                            tACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGt                                                                                                                                        <  2:57012/102‑1 (MQ=255)
                                                            tACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGt                                                                                                                                        >  1:57012/1‑102 (MQ=255)
                                                                   gggTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGAtct                                                                                 <  2:189016/148‑1 (MQ=255)
                                                                    ggTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGTGACCAGACGCTGGACCGACCCGATGATCTCa                                                                               <  1:24756/149‑1 (MQ=255)
                                                                      ttGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATc                                                                             >  5:58863/1‑149 (MQ=255)
                                                                            gTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCa                                                                                                           >  8:159989/1‑115 (MQ=255)
                                                                            gTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCa                                                                                                           <  7:159989/115‑1 (MQ=255)
                                                                                 cccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTc                                                                                                                    >  8:1431/1‑101 (MQ=255)
                                                                                 cccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTc                                                                                                                    <  7:1431/101‑1 (MQ=255)
                                                                                             gTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCtgt                                                        >  1:174119/1‑147 (MQ=255)
                                                                                             gTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCtgt                                                        <  2:174119/147‑1 (MQ=255)
                                                                                                        tCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAg                                           <  1:144527/149‑1 (MQ=255)
                                                                                                        tCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGACCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGaa                                            <  5:47850/148‑1 (MQ=255)
                                                                                                                       tGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGtt                             <  3:60055/148‑1 (MQ=255)
                                                                                                                       tGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTg                            <  6:58863/149‑1 (MQ=255)
                                                                                                                             aaCGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGc                         >  7:169931/1‑146 (MQ=255)
                                                                                                                             aaCGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGc                         <  8:169931/146‑1 (MQ=255)
                                                                                                                                      aggggcgCCACGGGCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAggg             >  5:40806/4‑149 (MQ=255)
                                                                                                                                      aggggcgCCACGGGCGTGGCGCCCCTGGTGATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGCCCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGGGCGCAGATCGACGGCCGTCTGCCGGTTGAGCAAGAACTCAggg             >  7:93991/4‑149 (MQ=255)
                                                                                                                                              cACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGAGATCATCgggtc                                                                        >  6:65482/1‑78 (MQ=255)
                                                                                                                                              cACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGAGATCATCgggtc                                                                        <  5:65482/82‑5 (MQ=255)
                                                                                                                                                   ccGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGAATGCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCAg  >  1:80494/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTAGTCGGTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCAG  >  NC_002947/1777378‑1777671

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: