Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,400,023 C→T 9.7% V379I (GTT→ATT)  PP_3879 ← HK97 family portal protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,400,0230CT 9.7% 111.4 / 2.5 41V379I (GTT→ATT) PP_3879HK97 family portal protein
Reads supporting (aligned to +/- strand):  ref base C (18/19);  new base T (2/2);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

ATGTCGTTTCATGCGGCAAGCATCCTCAGGATGTCGTCGCGCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCAGAGAGGTTGAACTCGATGAAGTAATCCCGCCGCTCC  >  NC_002947/4399878‑4400168
                                                                                                                                                 |                                                                                                                                                 
aTGTCGTTTCATGCGGCAAGCATCCTCAGGATGTCGTCGCGCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAATggg                                                                                                                                                >  3:197484/1‑149 (MQ=255)
                         tCAGGATGTCGTCGCGCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATg                                                                                                                       >  4:305462/1‑149 (MQ=255)
                             gATGTCGTCGCGCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGt                                                                                                                   <  4:229601/149‑1 (MQ=255)
                              aTGTCGTCGCTCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACTGGCGGCATGTTTTCCAGGCGCCGGATGTCGtt                                                                                                                  <  8:199941/149‑1 (MQ=255)
                              aTGTCGTCGCGCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGtt                                                                                                                  <  5:312218/149‑1 (MQ=255)
                              aTGTCGTCGCGCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGtt                                                                                                                  <  3:174776/149‑1 (MQ=255)
                                            ctgTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCAcc                                                                                                    <  5:233355/149‑1 (MQ=255)
                                            ctgTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCAcc                                                                                                    >  6:233355/1‑149 (MQ=255)
                                              gTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCAcccc                                                                                                  >  2:280098/1‑149 (MQ=255)
                                                   aGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTg                                                                                             >  4:274191/1‑149 (MQ=255)
                                                    gCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTAGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGc                                                                                            >  6:12128/1‑149 (MQ=255)
                                                        ggCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGc                                                                                        <  1:280098/149‑1 (MQ=255)
                                                               aCGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGc                                                                                 <  5:123430/149‑1 (MQ=255)
                                                                        ttAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTAGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGcgtc                                                                        <  5:12128/149‑1 (MQ=255)
                                                                           gTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGcgtcg                                                                       >  1:126898/1‑147 (MQ=255)
                                                                           gTCATGGCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAATGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGcgtcgt                                                                      <  4:197484/148‑1 (MQ=255)
                                                                              aTGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACGGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCg                                                                  >  7:54280/1‑149 (MQ=255)
                                                                              aTGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCg                                                                  >  7:265339/1‑149 (MQ=255)
                                                                              aTGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCg                                                                  >  6:199521/1‑149 (MQ=255)
                                                                               tGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAATGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCATCGTAGCgg                                                                 >  3:6076/1‑149 (MQ=255)
                                                                                 tCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCt                                                               >  6:323307/1‑149 (MQ=255)
                                                                                   gggCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCTTCGTAGCGGCttt                                                             >  6:73156/1‑149 (MQ=255)
                                                                                      cAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCt                                                          >  3:279279/1‑149 (MQ=255)
                                                                                       aGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTg                                                         <  2:126898/149‑1 (MQ=255)
                                                                                        ggTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTgg                                                        >  1:155998/1‑149 (MQ=255)
                                                                                         gTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGt                                                       <  8:265339/149‑1 (MQ=255)
                                                                                            aTTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTTCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTcgc                                                    <  5:199521/149‑1 (MQ=255)
                                                                                            aTTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTcgc                                                    >  2:96884/1‑149 (MQ=255)
                                                                                            aTTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTcgc                                                    <  1:96884/149‑1 (MQ=255)
                                                                                            aTTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTcg                                                     <  5:323307/148‑1 (MQ=255)
                                                                                                  tACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTCGCcgcg                                                >  7:106270/1‑147 (MQ=255)
                                                                                                          tCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCa                                      <  3:274191/149‑1 (MQ=255)
                                                                                                            aTGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCaga                                    <  2:155998/149‑1 (MQ=255)
                                                                                                               ttGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCTTCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCagagag                                 <  5:73156/149‑1 (MQ=255)
                                                                                                                  aGCGGCTGCAGATAACTGTCGCCGTTGGCAATGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCATCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCAGAGAGGtt                              <  4:6076/149‑1 (MQ=255)
                                                                                                                       cTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCAGAGAGGTTGAACt                         <  4:279279/149‑1 (MQ=255)
                                                                                                                              tAACTGTCGCCGTTGGCAACGGGCGGCCTGATTGCCAGGCGCCGGATGTCGTTGATCGACAGCCCCCCACCCTGCCGGCCCATTGCATAGGCGGCGTAGCGGCTTTTCGGGTCGCCGCCCAGCCGGCCAGAGCGGTGGACCTCGagg                    >  2:239752/1‑145 (MQ=255)
                                                                                                                               aaCTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCAGAGAGGTTGAACTCGATGAAg                 <  6:42989/149‑1 (MQ=255)
                                                                                                                                               aaCGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCAGAGAGGTTGAACTCGATGAAGTAATCCCGCCGCTcc  >  4:76639/1‑148 (MQ=255)
                                                                                                                                               aaCGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCAGAGAGGTTGAACTCGATGAAGTAATCCCGCCGCTc   <  7:233625/147‑1 (MQ=255)
                                                                                                                                               aaCGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCAGAGAGGTTGAACTCGATGAAGTAATCCCGCCGCTc   >  8:233625/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ATGTCGTTTCATGCGGCAAGCATCCTCAGGATGTCGTCGCGCTGCTGTTCCAGCTGGGCGCGGACGTTGGGGTTAGTCATGTCGGGCAGGTCATTGGCTACGTCGGTCATGTTGAGCGGCTGCAGATAACTGTCGCCGTTGGCAACGGGCGGCATGTTTTCCAGGCGCCGGATGTCGTTGATCGACAGCCACCCCCACTGCCGGCCAATTGCATAGGCGTCGTAGCGGCTTTTCTGGTCGCCGCGCAGCAGGCCAGAGAGGTTGAACTCGATGAAGTAATCCCGCCGCTCC  >  NC_002947/4399878‑4400168

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