Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,054,732 C→G 13.1% 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.1% 97.3 / 4.2 37intergenic (‑207/+101)PP_2666/PP_2667hypothetical protein/ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (12/20);  new base G (2/3);  total (14/23)
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.88e-01

ACGATAGACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGA  >  NC_002947/3054586‑3054870
                                                                                                                                                  |                                                                                                                                          
acgaTAGACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAg                                                                                                                                          >  2:104493/1‑149 (MQ=255)
  gaTAGACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCg                                                                                                                                        <  5:167396/149‑1 (MQ=255)
      gACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCATGCGTGGCGAAAAGGCCGCGCAGCCGccc                                                                                                                                    >  7:206977/1‑149 (MQ=255)
         aGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAg                                                                                                                                 >  8:365296/1‑149 (MQ=255)
          gTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCTGCGCGGccctttc                                                                                                                                <  8:169622/149‑5 (MQ=255)
          gTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCTGCGCGGccctttc                                                                                                                                <  8:163855/149‑5 (MQ=255)
          gTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCTGCGCGGccctttc                                                                                                                                <  3:183191/149‑5 (MQ=255)
                gCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCg                                                                                                                          >  5:126520/1‑149 (MQ=255)
                 cAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCg                                                                                                                          <  1:104493/148‑1 (MQ=255)
                  aGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCGACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTGGCGAAAGGGCCGCGCAGCCGCCCCAGGAtgtggg                                                                                                                         >  5:368240/1‑143 (MQ=255)
                            tagtagATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATa                                                                                                              <  4:181016/149‑1 (MQ=255)
                                                 aaGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGc                                                                                          <  3:88889/148‑1 (MQ=255)
                                                      aTGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGGCGTGCGTGGCGAAAGGGCCGCGGCCCGGCCCCCGGAGCTGGGGCGCAAGCAATGGTTGCGGGGGGCGCGTTGGGGgcc                                                                                     >  2:338102/1‑148 (MQ=12)
                                                      aTGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGCCGCAAAGCGGcccca                                                                                   <  7:365296/149‑3 (MQ=37)
                                                         gagaCGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAAt                                                                                  <  6:126520/148‑1 (MQ=255)
                                                              gCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCgg                                                                                        >  1:301574/1‑137 (MQ=255)
                                                              gCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCgg                                                                                        <  2:301574/137‑1 (MQ=255)
                                                                               tCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCcac                                                            <  4:243278/148‑1 (MQ=255)
                                                                                   tcTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACgg                                                        <  8:218450/147‑1 (MQ=255)
                                                                                   gCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACgg                                                        <  1:338102/148‑1 (MQ=255)
                                                                                             ttaaaGCTGGCGCTATCCCTGTGGGAGCGTCATTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCa                                             <  7:304315/147‑1 (MQ=255)
                                                                                                      gcgcGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTc                                      <  8:206977/147‑1 (MQ=255)
                                                                                                         cgATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAgcgc                                 <  6:170590/149‑1 (MQ=255)
                                                                                                                          ctcTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGc                <  1:374396/147‑1 (MQ=255)
                                                                                                                           ccTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCa               <  4:118940/149‑1 (MQ=255)
                                                                                                                            cTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCAt              <  2:230271/149‑1 (MQ=255)
                                                                                                                               gCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCt           >  8:357590/1‑149 (MQ=255)
                                                                                                                                 gTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAg                      <  5:383138/136‑1 (MQ=255)
                                                                                                                                 gTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAg                      >  6:383138/1‑136 (MQ=255)
                                                                                                                                   cgcgAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTggcg       >  4:5370/1‑147 (MQ=255)
                                                                                                                                   cgcgAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCAACATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTgggg       >  6:253747/1‑149 (MQ=255)
                                                                                                                                      cctggggCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGtttt                                       >  3:6143/5‑114 (MQ=255)
                                                                                                                                      cctggggCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGtttt                                       <  4:6143/110‑1 (MQ=255)
                                                                                                                                        tggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCAAGGTTTTCACCGCGGCTTGCGAAGCGCTGCATGCAGGCGGATg   >  4:21715/3‑148 (MQ=37)
                                                                                                                                        aaGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCACGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGa  <  6:368240/149‑1 (MQ=255)
                                                                                                                                             ccGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCg                                                                   <  1:360170/79‑1 (MQ=255)
                                                                                                                                             ccGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCg                                                                   >  2:360170/1‑79 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
ACGATAGACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGA  >  NC_002947/3054586‑3054870

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: