Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,442,214 T→C 9.6% H106R (CAT→CGT) ‡ nicP‑II ← porin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,442,2140TC 9.6% 136.1 / 5.3 52H106R (CAT→CGT) ‡nicP‑IIporin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (27/20);  new base C (3/2);  total (30/22)
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.03e-01

TGAACAGCAGCACCGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTGGGTGCTGGCGCCCT  >  NC_002947/4442068‑4442358
                                                                                                                                                  |                                                                                                                                                
tGAACTGCAGCACGGGTAGGTTTGGCATGAACTCGCCCTCTTTCAGCTCGGTTTGCGACACCTTGATCTTTTTCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCGGGACCGGGCAt                                                                                                                                                  <  8:183403/147‑1 (MQ=255)
tGAACAGCAGCACCGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGt                                                                                                                                                >  6:184052/1‑149 (MQ=255)
             cGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAGCCCGGAACATGCCCGGTCCGGGttgctgc                                                                                                                                   <  2:69455/149‑2 (MQ=38)
             cGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAGCCCGGAACATGCCCGGTCCGGGttgctgc                                                                                                                                   <  6:429345/149‑2 (MQ=38)
             cGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGt                                                                                                                                   <  8:411535/149‑1 (MQ=255)
               gCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTc                                                                                                                                  >  4:317817/1‑148 (MQ=255)
                   gTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCt                                                                                                                             >  7:164594/1‑149 (MQ=255)
                         cATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAAc                                                                                                                       >  6:285748/1‑149 (MQ=255)
                             aaCTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCcag                                                                                                                   <  8:88765/149‑1 (MQ=255)
                               cTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCcag                                                                                                                   >  7:89132/1‑147 (MQ=255)
                                    ccACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCg                                                                                                            >  2:79039/1‑149 (MQ=255)
                                    ccACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCg                                                                                                            <  6:7676/149‑1 (MQ=255)
                                        tttCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTaa                                                                                                        >  1:41939/1‑149 (MQ=255)
                                             gCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATc                                                                                                    >  1:168505/1‑148 (MQ=255)
                                                 ggTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCcac                                                                                                 <  8:89132/147‑1 (MQ=255)
                                                   tttGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACAccc                                                                                             >  8:275562/1‑149 (MQ=255)
                                                       cGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCacac                                                                                               <  6:250563/143‑1 (MQ=255)
                                                       cGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCacac                                                                                               >  5:250563/1‑143 (MQ=255)
                                                        gACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTAGAAACCCAGCAGGCCGGTAACATCCACACCCACGc                                                                                         >  5:238844/1‑148 (MQ=255)
                                                            ccTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAgg                                                                                    <  2:168505/149‑1 (MQ=255)
                                                            ccTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAgg                                                                                    <  8:202729/149‑1 (MQ=255)
                                                            ccTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAgg                                                                                    <  3:317817/149‑1 (MQ=255)
                                                             cTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAAcc                                                                                                                      <  8:212902/114‑1 (MQ=255)
                                                             cTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAAcc                                                                                                                      >  7:212902/1‑114 (MQ=255)
                                                                gAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGt                                                                                  <  7:77302/147‑1 (MQ=255)
                                                                   cTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGCCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCt                                                                             >  4:10633/1‑149 (MQ=255)
                                                                      gATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTgg                                                                           <  5:285748/148‑1 (MQ=255)
                                                                        tCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAATCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGta                                                                        >  1:49282/1‑149 (MQ=255)
                                                                           ccgcACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATaa                                                                     >  2:135917/1‑149 (MQ=255)
                                                                              cACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCg                                                                  <  6:267839/149‑1 (MQ=255)
                                                                                cccAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGa                                                                <  8:164594/149‑1 (MQ=255)
                                                                                 ccAGGCGGCTGTATTCGTCCGCAGACTTGCCGCTGTCCAGGCTGGGAAGCAACCCGGACCGGGCATCTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAg                                                               >  2:162598/1‑149 (MQ=255)
                                                                                       ggcTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGACTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCa                                                         >  5:204063/1‑149 (MQ=255)
                                                                                            tATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCCCGCCAAGGGTGACCTGGGTATAACCGGAGCTGGCATGGAg                                                    >  4:291702/1‑149 (MQ=255)
                                                                                              ttCGTCCGCCGACTTGCCGCTGTCCAGGCTTGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAgg                                                   <  8:173443/148‑1 (MQ=255)
                                                                                              ttCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAATCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAgg                                                   <  2:49282/148‑1 (MQ=255)
                                                                                              ttCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAgg                                                   <  7:275562/148‑1 (MQ=255)
                                                                                              ttCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACACACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGa                                                  >  5:408172/1‑149 (MQ=255)
                                                                                                    cgccgACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAAcc                                            >  8:208633/1‑149 (MQ=255)
                                                                                                     gccgACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAAccc                                           >  5:29855/1‑149 (MQ=255)
                                                                                                         aCTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTgc                                       <  6:408172/149‑1 (MQ=255)
                                                                                                               cGCCGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCAt                                 >  4:339817/1‑149 (MQ=255)
                                                                                                                 cTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTc                               >  8:408767/1‑149 (MQ=255)
                                                                                                                  tGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGGCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATTCCCGGAGCTGGCAGTGAGGATAAACCCCTGCGCCCCTTcc                              >  4:195263/1‑149 (MQ=255)
                                                                                                                      cAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTgg                           <  7:287477/148‑1 (MQ=255)
                                                                                                                      cAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGCCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTgg                           <  3:10633/148‑1 (MQ=255)
                                                                                                                         gctcgacagcagCCCGGAACATGCCCGGTCCGGGTTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTgc                       >  5:162961/7‑149 (MQ=37)
                                                                                                                            cgacagcagCCCGGAACATGCCCGGTCCGGGTTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGcgcg                    >  8:343464/4‑149 (MQ=38)
                                                                                                                               cagcagCCCGGAACATGCCCGGTCCGGGTTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACt                 >  8:188386/1‑149 (MQ=255)
                                                                                                                                 gcaACCCGGACCGTGCATGTTCCTGGCTGCTGGCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTgg               <  6:238844/149‑1 (MQ=255)
                                                                                                                                      ccGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTGGGTGCt          >  5:220714/1‑149 (MQ=255)
                                                                                                                                         gACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTGGGTGCTGgc       >  1:293644/1‑149 (MQ=255)
                                                                                                                                              ggCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTGGGTGCTGGCGCCCt  >  5:249648/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGAACAGCAGCACCGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTGGGTGCTGGCGCCCT  >  NC_002947/4442068‑4442358

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