Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,088,091 G→C 15.0% intergenic (+34/+16) yiiG → / ← frvR DUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,088,0910GC15.0% 113.0 / 14.6 40intergenic (+34/+16)yiiG/frvRDUF3829 family lipoprotein/putative frv operon regulator; contains a PTS EIIA domain
Reads supporting (aligned to +/- strand):  ref base G (19/15);  new base C (2/4);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

GATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCG  >  NC_000913/4087947‑4088239
                                                                                                                                                |                                                                                                                                                    
gATAAAGAGCAACTTCAGGACGCAAACTAAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCa                                                                                                                                                  >  4:702938/1‑149 (MQ=255)
     agagCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGa                                                                                                                                             <  3:480110/149‑1 (MQ=255)
     agagCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGa                                                                                                                                             <  4:169674/149‑1 (MQ=255)
                aGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCAcc                                                                                                                                  >  4:759327/1‑149 (MQ=255)
                                      ggTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCtt                                                                                                                              <  4:1012347/131‑1 (MQ=255)
                                      ggTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCtt                                                                                                                              >  3:1012347/1‑131 (MQ=255)
                                          gAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                                                                        <  3:759327/149‑1 (MQ=255)
                                          gAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGcc                                                                                                        >  3:691036/1‑149 (MQ=255)
                                                tctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATTATAAATGCAGTGTCGGATGCGCCGCATCCGaca                                                                                                                                           <  4:694347/108‑1 (MQ=255)
                                                tctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATTATAAATGCAGTGTCGGATGCGCCGCATCCGaca                                                                                                                                           >  3:694347/1‑108 (MQ=255)
                                                tctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggc                                                                                                                                                     <  2:770892/98‑1 (MQ=255)
                                                tctTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATgcggc                                                                                                                                                     >  1:770892/1‑98 (MQ=255)
                                                   ttCGCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGa                                                                                                                                             <  4:966855/103‑1 (MQ=255)
                                                   ttCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGCCGCATCCGa                                                                                                                                             >  3:966855/1‑103 (MQ=255)
                                                       cGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                                                                           >  4:447265/1‑149 (MQ=255)
                                                              aTTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTc                                                                                                                  <  1:957620/119‑1 (MQ=255)
                                                              aTTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTc                                                                                                                  >  2:957620/1‑119 (MQ=255)
                                                                     gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGGTGCGGCGCATCCGAcaca                                                                                                                                         <  4:980379/89‑1 (MQ=255)
                                                                     gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGATTTTGATGATAAATGCAGTGTTGGGCGCGGCGCATCCGAcaca                                                                                                                                         >  3:980379/1‑89 (MQ=255)
                                                                            aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCTACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgccg                                                                      >  1:939275/1‑149 (MQ=255)
                                                                                          gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat                                                        <  3:19063/149‑1 (MQ=255)
                                                                                             tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa                                                              <  1:985100/140‑1 (MQ=255)
                                                                                             tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa                                                              >  2:985100/1‑140 (MQ=255)
                                                                                             tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg                                                     >  1:624767/1‑149 (MQ=255)
                                                                                                aataGCCTGCGTTAAACGTTTTAGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTAAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt                                                  <  3:702938/149‑1 (MQ=255)
                                                                                                     ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAg                                             >  3:898605/1‑149 (MQ=255)
                                                                                                      cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGa                                            >  4:331101/1‑149 (MQ=255)
                                                                                                       tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGAt                                           <  4:898605/149‑1 (MQ=255)
                                                                                                           ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCt                                       >  3:9647/1‑149 (MQ=255)
                                                                                                              aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTaaa                                    >  4:423466/1‑149 (MQ=255)
                                                                                                                  tttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCt                                       >  3:435729/1‑142 (MQ=255)
                                                                                                                  tttttGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCt                                       <  4:435729/142‑1 (MQ=255)
                                                                                                                               aTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                                                                         >  4:747137/1‑47 (MQ=255)
                                                                                                                               aTGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGaa                                                                                                                         <  3:747137/47‑1 (MQ=255)
                                                                                                                                 gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt                                                  >  2:953101/1‑116 (MQ=255)
                                                                                                                                 gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt                                                  <  1:953101/116‑1 (MQ=255)
                                                                                                                                 gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGAcgc                 <  2:624767/149‑1 (MQ=255)
                                                                                                                                 gCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGAcgc                 <  3:447265/149‑1 (MQ=255)
                                                                                                                                         ggATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAtttt                              <  1:399042/128‑1 (MQ=255)
                                                                                                                                         ggATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATAtttt                              >  2:399042/1‑128 (MQ=255)
                                                                                                                                           aTGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                                                                    <  2:569195/88‑1 (MQ=255)
                                                                                                                                           aTGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                                                                    >  1:569195/1‑88 (MQ=255)
                                                                                                                                                gcgcATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGgcggcggcg  >  3:152724/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GATAAAGAGCAACTTCAGGACGCAAACTCAAGCTGGATGGTCGAAGACTCTTTCCCGAGAGCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTAAATATTTTTAGCTCATGACGCGCATCGGCGGCGGCG  >  NC_000913/4087947‑4088239

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: