Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,955,253 T→A 20.0% intergenic (+40/‑153) PP_1752 → / → PP_1753 hypothetical protein/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,955,2530TA20.0% 50.3 / 9.3 25intergenic (+40/‑153)PP_1752/PP_1753hypothetical protein/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base T (8/12);  new base A (2/3);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

CTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAG  >  NC_002947/1955110‑1955399
                                                                                                                                               |                                                                                                                                                  
cTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCcgc                                                                                                                                              >  1:265886/1‑150 (MQ=255)
   ggACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTAAACCCGCGaa                                                                                                                                           <  1:2052/150‑1 (MQ=255)
      cGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                        <  1:265378/150‑1 (MQ=255)
           gtggAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAg                                                                                                                                   >  2:336385/1‑150 (MQ=255)
                   cAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGt                                                                                                                                      >  2:370224/1‑139 (MQ=255)
                   cAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGt                                                                                                                                      <  1:370224/139‑1 (MQ=255)
                         cAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCt                                                                                                                     <  2:107926/150‑1 (MQ=255)
                          aaGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTc                                                                                                                    <  1:336385/150‑1 (MQ=255)
                          aaGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTGCTGACCCTTTCGCGGGTAAACCCGCGaagcggcc                                                                                                                                     >  2:88401/1‑127 (MQ=25)
                          aaGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTGCTGACCCTTTCGCGGGTAAACCCGCGaagcggcc                                                                                                                                     <  1:88401/133‑7 (MQ=37)
                                        ggTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcc                                                                                                      <  2:179290/150‑1 (MQ=255)
                                                    gcgcACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGa                                                                                          <  1:158149/150‑1 (MQ=255)
                                                      gcACAATCAACCAGGCGCTTGAGCCGGCGTTTGCCAAGCACAAGAAATGGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGAcc                                                                                        <  2:28122/150‑1 (MQ=255)
                                                                      gcTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGCATGACCTGATGAATCCTTTtct                                                                        <  2:340738/150‑1 (MQ=255)
                                                                                        cAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcc                                                                                                      >  2:254580/1‑102 (MQ=255)
                                                                                        cAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcc                                                                                                      <  1:254580/102‑1 (MQ=255)
                                                                                         aaGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtctc                                                     >  2:268206/1‑150 (MQ=255)
                                                                                          aGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTctct                                                    <  2:265886/150‑1 (MQ=255)
                                                                                                                 aGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTca                             >  1:445499/1‑150 (MQ=255)
                                                                                                                                   gCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGc           >  2:379936/1‑150 (MQ=255)
                                                                                                                                    tttCGCGGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTa                                    >  1:393336/2‑124 (MQ=255)
                                                                                                                                    tttCGCGGGTAAACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTa                                    <  2:393336/123‑1 (MQ=255)
                                                                                                                                       cgcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGa       <  1:268206/150‑1 (MQ=255)
                                                                                                                                           ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa   <  2:196077/150‑1 (MQ=255)
                                                                                                                                            gTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAg  >  2:353310/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CTTGGACGCTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAG  >  NC_002947/1955110‑1955399

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: