Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,561,243 G→C 6.9% I102M (ATC→ATG)  PP_4048 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,561,2430GC 6.9% 185.8 / 2.9 58I102M (ATC→ATG) PP_4048hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base C (2/2);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

CACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGA  >  NC_002947/4561105‑4561385
                                                                                                                                          |                                                                                                                                              
cACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGa                                                                                                                                      >  8:221093/1‑149 (MQ=255)
       tGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGt                                                                                                                               >  6:325510/1‑149 (MQ=255)
            gATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGtcat                                                                                                                            >  3:266443/1‑147 (MQ=255)
            gATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGACGATGATCGGCGCGATCAGGGtcatca                                                                                                                          <  1:10021/149‑1 (MQ=255)
                    gACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGcc                                                                                                                  >  2:145951/1‑149 (MQ=255)
                       tttCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGccac                                                                                                                >  4:312849/1‑148 (MQ=255)
                       tttCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGccac                                                                                                                >  4:165920/1‑148 (MQ=255)
                       tttCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTACCTGCCCAACCACTGAGCCTGCCAACGAGCCCACGGGAGTGCCCACCGCCCCCAACACAGCGCCGATGATCCCCGCCACAAGGGTCATCAGGCAGGccac                                                                                                                >  7:109960/1‑148 (MQ=255)
                           aCGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCGGGGTCATCAGGCAGGCCACCAGc                                                                                                            >  7:32774/1‑148 (MQ=255)
                                tCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCGGGGTCATCAGGCAGGCCACCAGCAACAcc                                                                                                      <  8:32774/149‑1 (MQ=255)
                                tCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCATCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACAcc                                                                                                      <  7:138743/149‑1 (MQ=255)
                                 cAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACAcc                                                                                                      >  6:97127/1‑148 (MQ=255)
                                       caagcaCAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCActc                                                                                                 <  5:122163/147‑1 (MQ=255)
                                       caagcaCAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCActc                                                                                                 >  6:122163/1‑147 (MQ=255)
                                       caagcaCAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCtt                                                                                               <  3:165920/149‑1 (MQ=255)
                                       caagcaCAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCtt                                                                                               >  5:145842/1‑149 (MQ=255)
                                        aagcaCAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTg                                                                                              <  1:241259/149‑1 (MQ=255)
                                               aaCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTcc                                                                                       <  6:232635/149‑1 (MQ=255)
                                               aaCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATCATCGgcgctgtg                                                                                                                                    >  2:63180/1‑100 (MQ=255)
                                               aaCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATCATCGgcgctgtg                                                                                                                                    <  1:63180/104‑5 (MQ=255)
                                                              ctgctgctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCTACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGtt                                                                        <  2:164008/149‑1 (MQ=255)
                                                                 ctgctgCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGc                                                                      <  8:210922/148‑1 (MQ=255)
                                                                 ctgctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGc                                                                      <  7:221093/148‑1 (MQ=255)
                                                                    ctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCt                                                                    <  6:145842/147‑1 (MQ=255)
                                                                    ctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTgg                                                                  >  1:312856/1‑149 (MQ=255)
                                                                      tCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCt                                                                >  4:102038/1‑149 (MQ=255)
                                                                       ccTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGTCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCt                                                                <  5:325510/148‑1 (MQ=255)
                                                                          gCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCg                                                            >  2:152281/1‑149 (MQ=255)
                                                                          gCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGACTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCg                                                            >  1:50555/1‑149 (MQ=255)
                                                                             cAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTcc                                                                                       >  2:145875/1‑119 (MQ=255)
                                                                             cAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTcc                                                         >  2:228525/1‑149 (MQ=255)
                                                                                      cgcCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTa                                                 <  1:152281/148‑1 (MQ=255)
                                                                                      cgcCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTa                                                 <  2:295195/148‑1 (MQ=255)
                                                                                       gcCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcag                                                                                                                         >  8:184096/1‑75 (MQ=255)
                                                                                       gcCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcag                                                                                                                         <  7:184096/75‑1 (MQ=255)
                                                                                       gcCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAgc                                               >  7:241191/1‑149 (MQ=255)
                                                                                         cTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAgcgc                                             >  6:283069/1‑149 (MQ=255)
                                                                                            ccATCGCGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcggc                                                                                    <  8:221883/107‑1 (MQ=255)
                                                                                            ccATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCtt                                           <  5:283069/148‑1 (MQ=255)
                                                                                              aTCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGcc                                        >  2:263825/1‑149 (MQ=255)
                                                                                               tCGAGCCCACGGGAGTGCCCACAGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGcc                                        >  6:293852/1‑148 (MQ=255)
                                                                                                    cccACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGgcg                                  >  8:117609/1‑149 (MQ=255)
                                                                                                    cccACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGgcg                                  >  7:35390/1‑149 (MQ=255)
                                                                                                    cccACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGgcg                                  >  1:237071/1‑149 (MQ=255)
                                                                                                     ccACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGgcgc                                 >  3:18379/1‑149 (MQ=255)
                                                                                                       aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCtt                                           <  1:142211/137‑1 (MQ=255)
                                                                                                       aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCtt                                           >  2:142211/1‑137 (MQ=255)
                                                                                                       aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAg                               <  2:312856/149‑1 (MQ=255)
                                                                                                         gggAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGc                             <  8:241191/149‑1 (MQ=255)
                                                                                                             gTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCacac                         <  7:117609/149‑1 (MQ=255)
                                                                                                                ccccccGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGc                      <  4:18379/149‑1 (MQ=255)
                                                                                                                 ccACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCa                     >  3:260574/1‑149 (MQ=255)
                                                                                                                   aCCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGc                   >  3:139514/1‑149 (MQ=255)
                                                                                                                            ctgatcgcgcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACg          >  6:235743/7‑149 (MQ=255)
                                                                                                                             ccacAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACgg         <  5:97127/149‑1 (MQ=255)
                                                                                                                               acaGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGCCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACggg        >  4:188715/1‑148 (MQ=255)
                                                                                                                                    gcCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGa  <  4:139514/149‑1 (MQ=255)
                                                                                                                                    gcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGa  <  5:235743/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                              
CACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGA  >  NC_002947/4561105‑4561385

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