Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,294,502 T→A 10.5% E155V (GAA→GTA)  PP_2894 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,294,5020TA10.5% 102.6 / 3.9 38E155V (GAA→GTA) PP_2894hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/18);  new base A (2/2);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

GACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTGCGGCCTTGCTCTCGCTCGCG  >  NC_002947/3294365‑3294648
                                                                                                                                         |                                                                                                                                                  
gACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCt                                                                                                                                                 >  1:26768/1‑141 (MQ=255)
gACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCt                                                                                                                                                 <  2:26768/141‑1 (MQ=255)
           cAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGAGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAg                                                                                                                                 <  4:122763/146‑1 (MQ=255)
                 cgACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTgg                                                                                                                        <  3:168673/149‑1 (MQ=255)
                  gACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGc                                                                                                                       >  1:137534/1‑149 (MQ=255)
                   aCTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCg                                                                                                                      <  2:127046/149‑1 (MQ=255)
                               gatgCGTGGCGGCCATCTAGGGACGGTGGGCGGCTAGATGAACTGGTTGGTGTCGGGTGAGGGCGGTTGGGTGGTGGCGCGTACCGGGACGATGAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCt                                                                                                          <  2:295377/146‑1 (MQ=255)
                                   cGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGaa                                                                                                      <  7:45095/149‑1 (MQ=255)
                                      ggcggcCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGtt                                                                                                                   <  8:1987/133‑1 (MQ=255)
                                      ggcggcCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGtt                                                                                                                   >  7:1987/1‑133 (MQ=255)
                                         ggcCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGtt                                                                                                                   >  2:318814/1‑130 (MQ=255)
                                         ggcCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGtt                                                                                                                   <  1:318814/130‑1 (MQ=255)
                                                aaCGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGGGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGAt                                                                                          <  8:269108/148‑1 (MQ=255)
                                                     acggTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATACCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCAc                                                                                    >  3:267451/1‑149 (MQ=255)
                                                                    aTGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATACCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTcgcg                                                                     >  3:38254/1‑149 (MQ=255)
                                                                         cTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTcgcg                                                                     <  6:349494/144‑1 (MQ=255)
                                                                         cTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTcgcg                                                                     >  5:349494/1‑144 (MQ=255)
                                                                         cTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGtttt                                                                >  5:101804/1‑149 (MQ=255)
                                                                          tCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTc                                                                        >  1:151155/1‑140 (MQ=255)
                                                                          tCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTc                                                                        <  2:151155/140‑1 (MQ=255)
                                                                            gTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGc                                                             <  6:101804/149‑1 (MQ=255)
                                                                             ttGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATACCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGcc                                                            <  4:38254/149‑1 (MQ=255)
                                                                                 tgtCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATATGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGa                                                        >  1:340994/1‑149 (MQ=255)
                                                                                      ggTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGATGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGcc                                                   >  1:54584/1‑149 (MQ=255)
                                                                                            ggCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTc                                                               <  7:180429/131‑1 (MQ=255)
                                                                                            ggCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTc                                                               >  8:180429/1‑131 (MQ=255)
                                                                                                   gggtggTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATg                                      >  3:86311/1‑149 (MQ=255)
                                                                                                              cgTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTg                             >  3:264361/1‑147 (MQ=255)
                                                                                                                   cGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTg                      >  2:86948/1‑149 (MQ=255)
                                                                                                                   cGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTg                      <  3:3892/149‑1 (MQ=255)
                                                                                                                     ggACGATCAGGTTGTTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTGCg                    <  2:54584/149‑1 (MQ=255)
                                                                                                                         gATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTGCGGCCt                <  4:264361/149‑1 (MQ=255)
                                                                                                                          aTCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTGCGGCCtt               >  4:208547/1‑149 (MQ=255)
                                                                                                                           tCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTGCGGCCTTg              <  1:86948/149‑1 (MQ=255)
                                                                                                                                   gtgtATACCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTGCGGCCTTGCTctcgc       <  4:267451/148‑1 (MQ=255)
                                                                                                                                     gtATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTGCGGCCTTGc             >  7:176315/1‑140 (MQ=255)
                                                                                                                                     gtATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTGCGGCCTTGc             <  8:176315/140‑1 (MQ=255)
                                                                                                                                       aTTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTGCGGCCTTGCTCTCGCTcgcg  >  2:59436/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
GACATCTTCTACAAGCGCGACTTCTAGGCCCGCAGCGTGGCGGCCATCAACGGACGGTGGCCGCCTCGATGAACTCGTTGGTGTCGGGTGAGGGCGGTGGGGTGGTGGCGCGTACCGGGACGATCAGGTTGGTGTATTCCTCGATCAACCGGTTCAGGGCCTCTGGCGGTTCGCTCGCCTGGAACTGCATGCGGATCTGCACTTCGTCCGGGTCGCGGTTTTCGCCCTGACGGCCTTGGCCACCAATGGTGACCTTGAAGCTTGCGGCCTTGCTCTCGCTCGCG  >  NC_002947/3294365‑3294648

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