Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,230,081 G→C 11.6% P61R (CCG→CGG) ‡ PP_4607 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,230,0810GC11.6% 121.1 / 8.6 43P61R (CCG→CGG) ‡PP_4607transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base C (2/3);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

CACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCA  >  NC_002947/5229935‑5230221
                                                                                                                                                  |                                                                                                                                            
caccaGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCgggg                                                                                                                                           <  2:994311/150‑1 (MQ=255)
 accaGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGggc                                                                                                                                          <  1:406249/150‑1 (MQ=255)
     gCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGgcggcg                                                                                                                                      >  2:387409/1‑150 (MQ=255)
       gACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCt                                                                                                                                    <  1:737353/150‑1 (MQ=255)
       gACGCCCACGACCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCCCGGCGgcgcg                                                                                                                                    <  2:926923/150‑2 (MQ=255)
        aCGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCCCGGCGgcncgg                                                                                                                                   <  2:700193/150‑3 (MQ=255)
        aCGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCCCCGgcggcgcggctggc                                                                                                                                   <  1:95497/150‑5 (MQ=255)
         cGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCttt                                                                                                                                  <  2:215970/150‑1 (MQ=255)
            ccACGGCCTGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCAGCCGGCGCCAGCACGCCGCCACTGTCGACATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCt                                                                                                                               <  1:684519/150‑1 (MQ=255)
            ccACGGCCAGCACACGCCAGGGTGCCAACGCCAGGCCTGCAACCAGAAGCCAGTAGCAGTTCTTCTGTGATCTCGCCCTGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCt                                                                                                                               <  1:684587/150‑1 (MQ=255)
             cACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTg                                                                                                                              >  1:858624/1‑150 (MQ=255)
                      cACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCAc                                                                                                                     <  2:692247/150‑1 (MQ=255)
                       aCCAGGCAGGCTGCCAACGCGAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCAcc                                                                                                                    >  1:563223/1‑150 (MQ=255)
                                    ccAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGCCGt                                                                                                       >  2:158219/1‑150 (MQ=255)
                                     cAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTa                                                                                                      <  1:157074/150‑1 (MQ=255)
                                       aCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAgg                                                                                                    <  2:48931/150‑1 (MQ=255)
                                       aCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAgg                                                                                                    <  2:23504/150‑1 (MQ=255)
                                       aCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAgg                                                                                                    <  1:114175/150‑1 (MQ=255)
                                                    cACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGc                                                                                       >  2:511113/1‑150 (MQ=255)
                                                     aCCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCCTGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTCCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGCGCCCGGTCGTAGGCTGCGCggcctgcc                                                                                      >  2:168237/1‑144 (MQ=255)
                                                     aCCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCg                                                                                      >  1:732174/1‑150 (MQ=255)
                                                     aCCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCg                                                                                      <  1:1032526/150‑1 (MQ=255)
                                                                   accacGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCaa                                                                        >  2:520162/1‑150 (MQ=255)
                                                                       cgcgCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGTGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAg                                                                    <  2:532530/150‑1 (MQ=255)
                                                                                 gtCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTCTAAGTTGTCGGCCAGCAACCAGCGCTTGAATg                                                          >  1:840799/1‑150 (MQ=255)
                                                                                  tCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgc                                                         >  1:682075/1‑150 (MQ=255)
                                                                                     ggCCATGGCATCGAGGTAGCGATGTAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgccgc                                                      <  2:563223/150‑1 (MQ=255)
                                                                                        cATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGACGTAGGCTGCGCGGTGTGCGAGGTCGGCCAGCAACCAGCGCATGAATGCCGcccc                                                   >  2:166971/1‑150 (MQ=255)
                                                                                            gCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCg                                               >  1:122404/1‑150 (MQ=255)
                                                                                              aTCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGc                                             >  2:1000886/1‑150 (MQ=255)
                                                                                              aTCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGc                                             >  1:9480/1‑150 (MQ=255)
                                                                                                         gcTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGtt                                  >  1:99617/1‑150 (MQ=255)
                                                                                                           tGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGc                                >  1:821367/1‑150 (MQ=255)
                                                                                                              aGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCg                             <  1:511113/150‑1 (MQ=255)
                                                                                                               gggCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGt                            >  1:780996/1‑150 (MQ=255)
                                                                                                                       tGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCGCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCgg                    >  2:273640/1‑150 (MQ=255)
                                                                                                                        gCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCtt                                                               >  1:902353/1‑106 (MQ=255)
                                                                                                                        gCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCtt                                                               <  2:902353/106‑1 (MQ=255)
                                                                                                                                   gCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGt                                   <  1:396009/123‑1 (MQ=255)
                                                                                                                                   gCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGt                                   >  2:396009/1‑123 (MQ=255)
                                                                                                                                    ccAGCCGCGCCGCCGGGGCGGCTCTGTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAGTGCCGCCGCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCa       <  2:556233/150‑1 (MQ=255)
                                                                                                                                        aagcgccgcCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGGGAGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGc   >  1:690255/3‑150 (MQ=255)
                                                                                                                                         agcgccgcCCCGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTACTCGGGGATCCGGTTGAACCATTGCAGGGCa  >  1:841133/2‑150 (MQ=255)
                                                                                                                                          gcgccgcCCCGGCGTCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgccgc                                                      <  1:911763/97‑1 (MQ=255)
                                                                                                                                          gcgccgcCCCGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgccgc                                                      >  2:911763/1‑97 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCA  >  NC_002947/5229935‑5230221

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