Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,166,004 C→A 8.7% G33V (GGC→GTC)  PP_3666 ← MHS family metabolite MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,166,0040CA 8.7% 132.9 / 4.0 46G33V (GGC→GTC) PP_3666MHS family metabolite MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base A (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

CCAGCGGGCGCACCAGGAACGACAGGGCGAATACCGCAAAGGTGGCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTCTTTTAATTGGAGGCAGGGTCTGG  >  NC_002947/4165858‑4166143
                                                                                                                                                  |                                                                                                                                            
ccAGCGGGCGCACCAGGAACGACAGGGCGAATACCGCAAAGGTGGCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGcc                                                                                                                                             >  1:230728/1‑148 (MQ=255)
             cAGGAACGACAGGGCGAATACCGCAAAGGTGGCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGccc                                                                                                                               >  1:158684/1‑149 (MQ=255)
                 aaCGACAGGGCGAATACCGCAAAGGTGGCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGggc                                                                                                                           <  6:360708/149‑1 (MQ=255)
                      cAGGGCGAATACCGCAAAGGTGGCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCt                                                                                                                      <  3:138891/149‑1 (MQ=255)
                           cGAATACCGCAAACGTGGCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGACGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTgcgc                                                                                                                 <  2:274190/149‑1 (MQ=255)
                                     aaaGGTGGCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGc                                                                                                       <  5:233085/149‑1 (MQ=255)
                                         gTGGCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTc                                                                                                   >  2:259203/1‑149 (MQ=255)
                                            gCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCggg                                                                                                                            <  1:100134/121‑1 (MQ=255)
                                            gCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCggg                                                                                                                            >  2:100134/1‑121 (MQ=255)
                                              cAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGc                                                                                                                                 <  7:72973/114‑1 (MQ=255)
                                              cAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGc                                                                                                                                 >  8:72973/1‑114 (MQ=255)
                                               aGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGtt                                                                                             >  6:266403/1‑149 (MQ=255)
                                                   aaCCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTtgtgt                                                                                         <  2:76916/149‑1 (MQ=255)
                                                                gTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCgg                                                                               <  2:230728/146‑1 (MQ=255)
                                                                  cAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATg                                                                         <  4:283558/149‑1 (MQ=255)
                                                                          gggCAGAAGGTGGCGGCGATGAGGGGGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATgg                                                                 <  4:208359/149‑1 (MQ=255)
                                                                          gggAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATgg                                                                 >  7:270843/1‑149 (MQ=255)
                                                                           ggAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATg                                                                         <  5:365775/140‑1 (MQ=255)
                                                                           ggAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATggg                                                                >  7:107299/1‑149 (MQ=255)
                                                                           ggAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGGTCGGCGTTTTGTGCGGGCCCGGCG‑ATg                                                                         >  6:365775/1‑140 (MQ=255)
                                                                                aagGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGtt                                                                                                    >  3:88156/1‑109 (MQ=255)
                                                                                aagGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGtt                                                                                                    <  4:88156/109‑1 (MQ=255)
                                                                                      gcggcgATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGggtgg                                                      <  8:270843/148‑1 (MQ=255)
                                                                                         gcgATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCGCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCt                                                  >  4:309246/1‑149 (MQ=255)
                                                                                         gcgATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGACCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCgg                                                                               >  7:98210/1‑121 (MQ=255)
                                                                                         gcgATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGACCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCgg                                                                               >  7:93318/1‑121 (MQ=255)
                                                                                         gcgATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGACCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCgg                                                                               <  8:93318/121‑1 (MQ=255)
                                                                                         gcgATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGACCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCgg                                                                               <  8:98210/121‑1 (MQ=255)
                                                                                             tgatgGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGAc                                              <  3:3100/149‑1 (MQ=255)
                                                                                              gatgGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACGACTCGACGAAGTTGCCCATGAAGCTGGCCCGTGAGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATTGCGGATGGGTCTAGGGTGGTGCTGGAc                                              <  2:158684/148‑1 (MQ=255)
                                                                                              gatgGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACa                                             >  5:69849/1‑149 (MQ=255)
                                                                                                tggtggCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAg                                                           <  5:92269/133‑1 (MQ=255)
                                                                                                tggtggCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAg                                                           >  6:92269/1‑133 (MQ=255)
                                                                                                             gCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCctct                              >  4:224226/1‑149 (MQ=255)
                                                                                                               cGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGGTGGGTGCGGGGCCGGCGTATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTg                             <  1:259203/149‑1 (MQ=255)
                                                                                                                      ggcgTAGTCGAACCACTCGCCGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCAGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTCtttt                     <  7:52616/149‑1 (MQ=255)
                                                                                                                         gTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTCTTTTAAt                  >  5:36300/1‑149 (MQ=255)
                                                                                                                            gTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTc                         <  7:270257/139‑1 (MQ=255)
                                                                                                                            gTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTc                         >  8:270257/1‑139 (MQ=255)
                                                                                                                               gAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGtt                                                                                                    <  6:212224/62‑1 (MQ=255)
                                                                                                                               gAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGtt                                                                                                    >  5:212224/1‑62 (MQ=255)
                                                                                                                               gAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTCTTTTAATTGGAg             <  6:69849/148‑1 (MQ=255)
                                                                                                                                    aCTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTCTTTTAATTGGAGGCAggg       >  7:331724/1‑149 (MQ=255)
                                                                                                                                     cTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTAGCCTCTGGTTCTTTTAATTGGAGGCAGGGt      >  5:320797/1‑149 (MQ=255)
                                                                                                                                      tCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTCTTTTAATTGGAGGCAGGGTc     >  3:222252/1‑149 (MQ=255)
                                                                                                                                        gacgaAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTCTTTTAATTGGAGGCAGGGTCtg   >  6:15601/1‑149 (MQ=255)
                                                                                                                                         acgaAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTCTTTTAATTGGAGGCAGGGTCtgg  >  3:36193/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
CCAGCGGGCGCACCAGGAACGACAGGGCGAATACCGCAAAGGTGGCCAGCAAACCGGTGGTCTTGTCAGTTTGCGGGAAGAAGGTGGCGGCGATGATGGTGGCCAGATAGCCGTAGGCGGCGTAGTCGAACCACTCGACGAAGTTGCCCATGAAGCTGGCCCGGGCGGCCTTGCGCAGCGCCTGGCGTTCGGCGTTGTGTGCGGGGCCGGCG‑ATGGCGGATGGGTCTAGGGTGGTGCTGGACATGACGGGTACCCTCTGGTTCTTTTAATTGGAGGCAGGGTCTGG  >  NC_002947/4165858‑4166143

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