Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,054,732 C→G 13.7% intergenic (‑207/+101) PP_2666 ← / ← PP_2667 hypothetical protein/ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,054,7320CG13.7% 109.1 / 10.4 44intergenic (‑207/+101)PP_2666/PP_2667hypothetical protein/ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (13/25);  new base G (3/3);  total (16/28)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

ACGATAGACAGTAAATGCAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGT  >  NC_002947/3054586‑3054874
                                                                                                                                                   |                                                                                                                                              
acgaTAGACAGTAAATGCAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAg                                                                                                                                              >  2:276436/1‑149 (MQ=255)
 cgaTAGACAGTAAATGCAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGc                                                                                                                                             <  7:186293/149‑1 (MQ=255)
   aTAGACAGTAAATGCAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGCAAGGGCCGCGCCGCgg                                                                                                                                           >  7:274306/1‑149 (MQ=255)
        cAGTAAATGCAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGGCGCGAAAAGGCCCCCCAACGGCCCCa                                                                                                                                      >  7:257733/1‑149 (MQ=37)
          gTAAATGCAGTACTAAAGTAGTAGATGGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGGGCGGCCTTGCGTGGGGAAAGGGGCCCCCAGCGGcccccg                                                                                                                                     >  1:221026/1‑147 (MQ=37)
          gTAAATGCAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAgg                                                                                                                                    >  7:62799/1‑149 (MQ=255)
          gTAAATGCAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAg                                                                                                                                     >  4:68832/1‑148 (MQ=255)
           tAAATGCAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCTGCGCGGccctttcg                                                                                                                                   <  3:287122/149‑6 (MQ=255)
           tAAATGCAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGa                                                                                                                                   >  3:28975/1‑149 (MQ=255)
             aaTGCAGTACTAAAGTACTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGAt                                                                                                                                  <  3:232638/148‑1 (MQ=255)
                 cAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCg                                                                                                                              <  3:55441/148‑1 (MQ=255)
                  aGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGc                                                                                                                            <  4:227727/149‑1 (MQ=255)
                         aaGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGcatcat                                                                                                                       <  4:28975/147‑1 (MQ=255)
                            tagtagAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATa                                                                                                                  <  7:229060/149‑1 (MQ=255)
                               tagAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGAt                                                                                                               <  4:226971/149‑1 (MQ=255)
                                       ccGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGGCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTggg                                                                                                        <  1:276436/149‑1 (MQ=255)
                                        cGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTgggg                                                                                                       <  3:68832/149‑1 (MQ=255)
                                                   aGGCATGGAGACGCCCTAGGGAGCAACTGCCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCgg                                                                                            <  8:155946/149‑1 (MQ=255)
                                                       aTGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGCAAGGGCCGCGCCGCGGGCCCCGGGGCGCGGGCTGAAGCATAGAGTGCTGGGGCCGCGTTGGGGcgc                                                                                         >  3:280194/1‑146 (MQ=21)
                                                                     gggAGCAACTGTCTTGCTCAATTTCTGATAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAgg                                                                          <  8:62799/149‑1 (MQ=255)
                                                                                     cTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGt                                                          <  2:132543/149‑1 (MQ=255)
                                                                                        attttCTGAAAGCTGGCTCGTTCGCTGTGTGAGCGGCCTTGTTTCGCGAAAGGGCCGCGCTGCGGACGCTGGATCTCGTCATCATGCTTAGATTGCTGGTGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTac                                                       <  4:335200/147‑1 (MQ=25)
                                                                                                aacgCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCa                                                 <  7:163118/144‑1 (MQ=255)
                                                                                                aaaGCTGGCGCGATCCCTGTGGGCGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAg                                                <  2:221026/148‑1 (MQ=255)
                                                                                                  aGCTGGCGCGATGCCTGTGGGAGCGGCCTGGCGTCGCGAATGGGCCGCGCTGCGGCCCCAGGAGCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGtt                                             <  8:257733/149‑1 (MQ=37)
                                                                                                       gcgcGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTc                                          <  3:111212/147‑1 (MQ=255)
                                                                                                                cTGGGGGAGCGGCCTTGCCTCGCGAACGGGCCGCGCAGCGGCCCCAGGATCTCGCCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTg                               <  7:179823/149‑1 (MQ=255)
                                                                                                                         cGGCATTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAgcgc                       <  4:234569/148‑1 (MQ=255)
                                                                                                                              ttGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATg                 >  3:203900/1‑149 (MQ=255)
                                                                                                                                gCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCgg                                   <  8:148779/129‑1 (MQ=255)
                                                                                                                                gCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCgg                                   >  7:148779/1‑129 (MQ=255)
                                                                                                                                  gTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTgg             <  3:44584/149‑1 (MQ=255)
                                                                                                                                     gcgAAAGGGCCGCGCAGCGGCCCCAGTATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTgggg           <  8:137463/148‑1 (MQ=255)
                                                                                                                                      cgAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGt         <  4:280194/149‑1 (MQ=255)
                                                                                                                                        ctggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGtt        >  3:265888/4‑148 (MQ=255)
                                                                                                                                         tggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCCCACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCt               >  6:14779/3‑140 (MQ=38)
                                                                                                                                         tggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCt               <  5:14779/138‑1 (MQ=255)
                                                                                                                                         tggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTg       >  6:105139/3‑148 (MQ=255)
                                                                                                                                           gggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGaa     <  4:265888/148‑1 (MQ=255)
                                                                                                                                             gCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCa                   <  2:265975/132‑1 (MQ=255)
                                                                                                                                             gCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCa                   >  1:265975/1‑132 (MQ=255)
                                                                                                                                               cgcgCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGt  >  6:381758/1‑147 (MQ=255)
                                                                                                                                               cgcgCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCAACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGa      <  6:5588/143‑1 (MQ=255)
                                                                                                                                               cgcgCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCAACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGa      >  5:5588/1‑143 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
ACGATAGACAGTAAATGCAGTACTAAAGTAGTAGAT‑GGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGT  >  NC_002947/3054586‑3054874

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