Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,085,849 G→A 8.7% intergenic (+515/‑80) PP_1864 → / → PP_1865 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,085,8490GA 8.7% 126.0 / 3.0 46intergenic (+515/‑80)PP_1864/PP_1865hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (18/24);  new base A (2/2);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

ACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAGG  >  NC_002947/2085700‑2085953
                                                                                                                                                     |                                                                                                        
aCTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCg                                                                                                          <  1:1717729/150‑1 (MQ=255)
aCTGGAGGTGAAAGCGCACAACCCAAAGGCCCAAAAACTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCg                                                                                                          >  2:1013198/1‑150 (MQ=255)
 cTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGc                                                                                                         <  1:673197/150‑1 (MQ=255)
           aaGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAg                                                                                               >  2:4019/1‑150 (MQ=255)
               gcACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAg                                                                                           <  2:720640/150‑1 (MQ=255)
               gcACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGCAGAATTCGCTGAAATCCAGTCAg                                                                                           <  2:949715/150‑1 (MQ=255)
                    aCCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGa                                                                                      >  2:523145/1‑150 (MQ=255)
                     cccAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAc                                                                                     <  2:1532461/150‑1 (MQ=255)
                       cAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAccc                                                                                   >  1:685460/1‑150 (MQ=255)
                       cAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAccc                                                                                   <  1:658670/150‑1 (MQ=255)
                        aaaGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCg                                                                                  >  2:809587/1‑150 (MQ=255)
                         aaGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGa                                                                                 <  1:1387777/150‑1 (MQ=255)
                         aaGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTCAACGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGa                                                                                 <  1:1040884/150‑1 (MQ=255)
                            gCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCg                                                                              >  1:484412/1‑150 (MQ=255)
                              ccAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGAt                                                                            >  1:1458130/1‑150 (MQ=255)
                              ccAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGAt                                                                            >  2:1369970/1‑150 (MQ=255)
                                   agctcaaaagctTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg                                                                       <  1:1642839/150‑1 (MQ=255)
                                        aaaagctTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCTAGGACGTCGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGAt                                                                  >  2:1215577/1‑150 (MQ=255)
                                        aaaagctTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGAt                                                                  <  1:4019/150‑1 (MQ=255)
                                              tttAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTa                                                            <  2:1365478/150‑1 (MQ=255)
                                                tAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAgc                                                          <  1:1013198/150‑1 (MQ=255)
                                                   cgcCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAgg                                                                                          <  2:67863/115‑1 (MQ=255)
                                                   cgcCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAgg                                                                                          >  1:67863/1‑115 (MQ=255)
                                                         ttCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCga                                                 <  1:762529/150‑1 (MQ=255)
                                                          tCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCgag                                                >  2:1609878/1‑150 (MQ=255)
                                                            tGATCGGGCGTACAGTCCGTTTGCGATGTCGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaa                                              >  2:996431/1‑150 (MQ=255)
                                                               tCGGGCGTACAGTCCGTTTGCGATGCCGCCGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGGCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaaaaa                                           >  1:1789967/1‑150 (MQ=255)
                                                                cGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaaaaaa                                          <  1:1369970/150‑1 (MQ=255)
                                                                                 tGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAAt                         <  2:158896/150‑1 (MQ=255)
                                                                                  gCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATg                        <  1:1422493/150‑1 (MQ=255)
                                                                                             gCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGc             <  2:281054/150‑1 (MQ=35)
                                                                                              ccAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGcc            <  1:1215577/150‑1 (MQ=35)
                                                                                                 aCTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTg         <  1:523145/150‑1 (MQ=32)
                                                                                                  cTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTgg        <  2:148922/150‑1 (MQ=32)
                                                                                                   tCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGt                                                                                              <  2:382917/63‑1 (MQ=37)
                                                                                                   tCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGt                                                                                              >  1:382917/1‑63 (MQ=35)
                                                                                                   tCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGa       >  2:302366/1‑150 (MQ=32)
                                                                                                   tCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGa       >  2:841146/1‑150 (MQ=32)
                                                                                                   tCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGa       >  1:539228/1‑150 (MQ=32)
                                                                                                     aCCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaaaaaa                                          <  2:984184/113‑1 (MQ=33)
                                                                                                     aCCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaaaaaa                                          >  1:984184/1‑113 (MQ=33)
                                                                                                      ccTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGAGCa    <  1:809587/150‑1 (MQ=32)
                                                                                                       cTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAg   >  2:204497/1‑150 (MQ=32)
                                                                                                       cTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCACTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaaaaaa                                          <  2:300913/111‑1 (MQ=17)
                                                                                                       cTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCACTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaaaaaa                                          >  1:300913/1‑111 (MQ=18)
                                                                                                       cTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCACTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaaaaaa                                          <  2:300168/111‑1 (MQ=18)
                                                                                                       cTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCACTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaaaaaa                                          >  1:300168/1‑111 (MQ=18)
                                                                                                        tttCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAgg  <  2:393803/150‑1 (MQ=32)
                                                                                                                                                     |                                                                                                        
ACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTTTGCGATGACGACGCCAACTCACCTTTCCGCCCTTACGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTGCCTTGGAGCAGG  >  NC_002947/2085700‑2085953

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: ‑