Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,375,442 C→A 13.5% A32S (GCG→TCG)  PP_5637 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,375,4420CA13.5% 97.1 / 8.0 37A32S (GCG→TCG) PP_5637hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/14);  new base A (2/3);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

CACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTG  >  NC_002947/4375294‑4375586
                                                                                                                                                    |                                                                                                                                                
caccacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTcgc                                                                                                                                                  <  3:5563/149‑1 (MQ=255)
caccacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTcgc                                                                                                                                                  <  4:107607/149‑1 (MQ=255)
caccacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTcgc                                                                                                                                                  >  5:218123/1‑149 (MQ=255)
   cacGACTGCGCATATAGCACGCCTTCGACCTCTTCAACCCCGTTGATGATGATGCCGAGCTGGGCCATGCCGGTGCCCTTCGCATCGTGCAGCCGGGGGATACTGTCGGGCCCAGGCGTAGGGGTCAACACCCAGGCATGGGTCGCGAc                                                                                                                                               >  2:248095/1‑149 (MQ=255)
   aacGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCAGTTGACCTTAGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAAGACACAGGACTGTGTCGCGa                                                                                                                                                >  4:31292/2‑148 (MQ=255)
             cATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTGGATCCCGCGCTGGGCCATGCCGTTCACCTCCGCAGCGTGCAGCCGGGGTATAATGGCGGGCCCCGGCGTAGGGTTAAACCCCCCGGCCTGGGTCGCGACTCGACCCaga                                                                                                                                     <  6:218123/149‑1 (MQ=255)
             cATAAAGCACGCGTTCCACCTCTTCTACCCGTTTGATTTTGTGGCCGCGCGGGGCCAGGCGTTTCCCCTTCGCTTGGTCCAGCGGGGGTATAATGTGGGGCCCCGGCTTTGGTTTAAACACCCCGGCCTGGGTCGCGACTCGACCCaga                                                                                                                                     <  2:164028/149‑1 (MQ=255)
              atatAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCagag                                                                                                                                    >  4:293102/1‑149 (MQ=255)
                 taGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGAGTCGCGACCCaggcctg                                                                                                                                 <  5:229595/149‑6 (MQ=255)
                   gCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGAGTCGCGACCCaggcctgg                                                                                                                                <  1:57838/148‑7 (MQ=255)
                        ccTTCGACCACTTCTCCCCCGTTGAGGTGGATGCCGCGCTGGGCCAGGCCGTTGACCTCCGCATCGGGCAGCCGGGGGATAATGTCGGGCCCTGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCtgtg                                                                                                                          <  4:11687/149‑1 (MQ=255)
                        ccTTCGACAGCTGATACCTCGTTGATGTTGTTGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGTGTATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCTGGCCTGGGGCGCGACTCGACCCAGAGGCTCGCtgtg                                                                                                                          <  6:201261/149‑1 (MQ=255)
                                        ccccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCagag                                                                                                                                    <  1:90169/123‑1 (MQ=255)
                                        ccccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCagag                                                                                                                                    >  2:90169/1‑123 (MQ=255)
                                         cccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGc                                                                                                                  <  1:118017/140‑1 (MQ=255)
                                         cccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGc                                                                                                                  >  2:118017/1‑140 (MQ=255)
                                                      atgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAg                                                                                            <  3:97242/149‑1 (MQ=255)
                                                       tgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAg                                                                                            <  3:293102/148‑1 (MQ=255)
                                                        gCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAggg                                                                                          >  1:134283/1‑149 (MQ=255)
                                                                     aTGCGGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCg                                                                             <  7:248338/149‑1 (MQ=255)
                                                                        ccGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCTTAGGGTTAAACACCAAGGCCTGTGTCGCGACTAGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAg                                                                          <  3:31292/149‑1 (MQ=255)
                                                                           ttGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCg                                                                       >  2:25810/1‑149 (MQ=255)
                                                                                      aTCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGc                                                             >  8:236766/1‑148 (MQ=255)
                                                                                                  ggggATAAGGTCGGGCCCAGGGGTAGGGTTAAACACCCAGGCCGGGGTGGCGAGTCGCCCCAGTGGCTCGCGGTGATGGTCGCCAATGTGGACGTCTGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCAt                                                <  1:248095/149‑1 (MQ=255)
                                                                                                                            ggTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCa                                                      >  3:276572/1‑117 (MQ=255)
                                                                                                                            ggTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCa                                                      <  4:276572/117‑1 (MQ=255)
                                                                                                                             gTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaaggaag                      >  6:129747/1‑148 (MQ=255)
                                                                                                                              ttAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGt                                                                            >  4:125508/1‑93 (MQ=255)
                                                                                                                              ttAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGt                                                                            <  3:125508/93‑1 (MQ=255)
                                                                                                                              ttAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTa                    >  6:148276/1‑149 (MQ=255)
                                                                                                                              ttAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTACGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTa                    >  2:63289/1‑149 (MQ=255)
                                                                                                                                    acCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCg              >  8:251584/1‑149 (MQ=255)
                                                                                                                                     cccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGAACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCgg             >  7:35676/1‑149 (MQ=255)
                                                                                                                                        aGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCAc          >  6:110310/1‑149 (MQ=255)
                                                                                                                                        aGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCAc          >  6:329542/1‑149 (MQ=255)
                                                                                                                                          ctcTGGGTCGAGTCGCGACCCAGGGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAg                                                                <  1:150955/91‑1 (MQ=38)
                                                                                                                                          ctcTGGGTCGAGTCGCGACCCAGGGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAg                                                                >  2:150955/3‑93 (MQ=38)
                                                                                                                                          ctcTGGGTCGAGTCGCGACCCAGGGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCAcc         >  7:304843/3‑148 (MQ=255)
                                                                                                                                               ggTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGAACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATa     <  8:35676/147‑1 (MQ=255)
                                                                                                                                                 tCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTg  >  6:306113/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CACCACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAGGGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTG  >  NC_002947/4375294‑4375586

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