Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,796,838 A→T 9.4% intergenic (+22/‑110) PP_t43 → / → PP_2449 tRNA‑Val/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,796,8380AT 9.4% 145.6 / 5.8 53intergenic (+22/‑110)PP_t43/PP_2449tRNA‑Val/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/21);  new base T (2/3);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

ACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGC  >  NC_002947/2796700‑2796982
                                                                                                                                          |                                                                                                                                                
aCAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGG‑TAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGccc                                                                                                                                        <  6:182051/148‑1 (MQ=17)
 cAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGccc                                                                                                                                        <  8:255654/148‑1 (MQ=255)
     tACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTTGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCATTTCTGATCGCCCtttt                                                                                                                                    <  3:205441/148‑1 (MQ=17)
      aCGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTg                                                                                                                                  >  4:8122/1‑149 (MQ=255)
       cGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGc                                                                                                                                 <  5:37956/149‑1 (MQ=255)
                ccACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGcc                                                                                                                        >  1:68677/1‑149 (MQ=255)
                ccACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGCCAAGAAGAAATGATCGCCCTTTTTGATTTTgcgcg                                                                                                                        >  8:237570/1‑148 (MQ=17)
                 cACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGccc                                                                                                                       >  5:64240/1‑149 (MQ=255)
                  aCTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGCTCAGAAATGATCGCCATTTTTTCATTAGCGccc                                                                                                                       >  4:189076/1‑148 (MQ=255)
                    tAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGaa                                                                                                                    <  2:326146/149‑1 (MQ=255)
                      atatGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGCAATGATCGCCCTTTTTGCATTTGCGCCCGAag                                                                                                                   >  8:307982/1‑148 (MQ=255)
                      atatGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAag                                                                                                                   >  8:3952/1‑148 (MQ=255)
                         tGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAagcagc                                                                                                               <  7:3952/149‑1 (MQ=255)
                                caccCGACAGGCAGATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAACTCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAAc                                                                                                        <  3:189076/149‑1 (MQ=255)
                                cacaCGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAAc                                                                                                        <  5:23324/149‑1 (MQ=255)
                                     gACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTg                                                                                                   >  5:80870/1‑149 (MQ=255)
                                                     aGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCa                                                                                    >  7:232663/1‑148 (MQ=255)
                                                      gTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATa                                                                                  >  2:245490/1‑149 (MQ=255)
                                                      gTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATa                                                                                  >  8:301113/1‑149 (MQ=255)
                                                           ttAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCCAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCg                                                                             >  3:135829/1‑149 (MQ=255)
                                                                gCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGt                                                                         <  4:135829/148‑1 (MQ=255)
                                                                 ccccaccaTGTCAAGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGc                                                                       >  4:152063/3‑149 (MQ=255)
                                                                 caccacCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGc                                                                       >  1:205829/1‑149 (MQ=255)
                                                                   ccaccTTGACATGGTGGGGGTCGTTGGTTCTAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGc                                                                       >  6:203859/1‑147 (MQ=255)
                                                                   ccaccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                     <  3:8122/149‑1 (MQ=255)
                                                                   ccaccTTGACATGGTGGGGGGCCGTTGTTCGCGTCCAAATGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                     <  7:237570/149‑1 (MQ=38)
                                                                      ccTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGaa                                                                  <  6:64240/149‑1 (MQ=255)
                                                                           aCATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGc                                                             <  7:301113/149‑1 (MQ=255)
                                                                            cATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGc                                                             <  8:232663/148‑1 (MQ=255)
                                                                              tggtggGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATcc                                                               >  3:313432/1‑144 (MQ=255)
                                                                                         gTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTgg                                               <  7:224818/149‑1 (MQ=255)
                                                                                          ttGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGAAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTg                                                >  3:159504/1‑147 (MQ=255)
                                                                                          ttGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTCGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGc                                                                           >  7:301339/1‑120 (MQ=255)
                                                                                          ttGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTCGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGc                                                                           <  8:301339/120‑1 (MQ=255)
                                                                                            ggTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAgg                                            >  4:271793/1‑149 (MQ=255)
                                                                                            ggTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAg                                             <  3:92014/148‑1 (MQ=255)
                                                                                                              ccTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACt                          >  1:142016/1‑149 (MQ=255)
                                                                                                                gcgcAAATGCAAAAAGGGCGATCATTTCTGATCGCCCTTTTTCGATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTgg                        >  8:200444/4‑149 (MQ=37)
                                                                                                                 aCCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGc                       >  6:113082/1‑149 (MQ=255)
                                                                                                                   cAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                     <  2:29072/101‑1 (MQ=255)
                                                                                                                   cAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTa                                                                     >  1:29072/1‑101 (MQ=255)
                                                                                                                   cAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCtg                     >  5:38583/1‑149 (MQ=255)
                                                                                                                       tCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCgg                                                                                            <  3:150704/74‑1 (MQ=255)
                                                                                                                       tCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCgg                                                                                            >  4:150704/1‑74 (MQ=255)
                                                                                                                          aaaaaGGGCGATCATTTCTGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAAt              >  4:199334/1‑149 (MQ=255)
                                                                                                                          aaaaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAAt              >  3:151745/1‑149 (MQ=255)
                                                                                                                            aaaGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGAAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATgg            <  4:159504/149‑1 (MQ=255)
                                                                                                                              aGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCa          <  6:80870/149‑1 (MQ=255)
                                                                                                                                 gCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACAcgc                                     <  1:9666/119‑1 (MQ=255)
                                                                                                                                 gCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACAcgc                                     >  2:9666/1‑119 (MQ=255)
                                                                                                                                     tCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTc    >  2:263145/1‑148 (MQ=255)
                                                                                                                                     tCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTc    <  1:263145/148‑1 (MQ=255)
                                                                                                                                       aGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGc  <  1:245490/148‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                
ACAGGTACGCATCGGGCCACTAATATGCGCACCACACGACAGGCACATAGCTCAGTTGGTTAGAGCACCACCTTGACATGGTGGGGGTCGTTGGTTCGAGTCCAATTGTGCCTACCAAATCGAAAAAGGGCGATCAGAAATGATCGCCCTTTTTGCATTTGCGCCCGAAGCAGCGTGTAACGTTTGGTCGCGGTCCGGGCATACCCCGGCGTGCTAGAATCCGCCCCTTCGAATCTGGAGGTACACGCGTGCGCAAACTGGCTGTGGTAATGGCAGTGCTCGC  >  NC_002947/2796700‑2796982

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