Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,516,520 C→T 7.4% G438D (GGT→GAT)  yfeA ← putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,516,5200CT 7.4% 163.6 / 3.3 54G438D (GGT→GAT) yfeAputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (24/26);  new base T (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

TAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTT  >  NC_000913/2516374‑2516665
                                                                                                                                                  |                                                                                                                                                 
tAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCAc                                                                                                                                                   >  7:352179/1‑147 (MQ=255)
   tGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATa                                                                                                                                               <  7:161683/148‑1 (MQ=255)
      tCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCATCATAGcc                                                                                                                                            >  8:7626/1‑148 (MQ=255)
      tCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCATCATAGcc                                                                                                                                            >  8:14012/1‑148 (MQ=255)
         gCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGATGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCa                                                                                                                                          >  1:21893/1‑147 (MQ=255)
                      ctcctcTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTa                                                                                                                            >  6:143783/1‑148 (MQ=255)
                             ctGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGccc                                                                                                                                           <  6:165453/126‑1 (MQ=255)
                             ctGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGccc                                                                                                                                           >  5:165453/1‑126 (MQ=255)
                                cTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAg                                                                                                                   <  2:375467/147‑1 (MQ=255)
                                cTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAg                                                                                                                   <  1:208531/147‑1 (MQ=255)
                                            aGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGAc                                                                                                      >  8:158798/1‑148 (MQ=255)
                                              cacaTGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTa                                                                                                    >  2:408798/1‑148 (MQ=255)
                                               acaTGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTAt                                                                                                   <  4:340874/148‑1 (MQ=255)
                                                caTGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCCGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATACCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGACGACTAtt                                                                                                  >  6:411171/1‑148 (MQ=255)
                                                caTGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTAtt                                                                                                  <  7:179566/148‑1 (MQ=255)
                                                caTGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTAtt                                                                                                  >  6:344575/1‑148 (MQ=255)
                                                        tgtgCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCAGGTATTGTTACAGTGAATTTGCCGACTATTCAGGATGt                                                                                          >  6:326094/1‑148 (MQ=255)
                                                            cgcAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCATCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTaa                                                                                       <  7:7626/147‑1 (MQ=255)
                                                            cgcAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCATCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTaa                                                                                       <  7:14012/147‑1 (MQ=255)
                                                             gcAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAAc                                                                                      >  2:408993/1‑147 (MQ=255)
                                                                ggATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCata                                                                                  >  4:388471/1‑148 (MQ=255)
                                                                 gATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCatat                                                                                 <  2:21893/148‑1 (MQ=255)
                                                                      cTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTgg                                                                            <  5:143783/148‑1 (MQ=255)
                                                                      cTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTg                                                                             <  7:204386/147‑1 (MQ=255)
                                                                        ccGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAg                                                                          >  8:45453/1‑148 (MQ=255)
                                                                         cGCCTGCCAGCTTATGTGCCCCAACTAGGGTTGCAGGGTTTCCTGTTTTCCATCATAACGCCGGCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAg                                                                          <  8:249667/147‑1 (MQ=255)
                                                                           ccagccagCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGGTCCAGTGAATTTGCCGACTAGTCAGGATGTTAACCATATGCTGGAGTCg                                                                       >  1:306348/1‑148 (MQ=255)
                                                                                cagcTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCtt                                                                  >  7:258708/1‑148 (MQ=255)
                                                                                 agcTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTc                                                                 <  7:45453/148‑1 (MQ=255)
                                                                                    ttAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTgg                                                                            <  8:139475/134‑1 (MQ=255)
                                                                                    ttAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTgg                                                                            >  7:139475/1‑134 (MQ=255)
                                                                                       aCTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTc                                                           <  3:368007/148‑1 (MQ=255)
                                                                                       aCTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTc                                                           <  1:408993/148‑1 (MQ=255)
                                                                                              caacaaGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCg                                                    >  5:243828/1‑148 (MQ=255)
                                                                                               aacaaGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGc                                                   >  3:300271/1‑148 (MQ=255)
                                                                                                   aGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAg                                               <  7:158798/148‑1 (MQ=255)
                                                                                                         gCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTa                                          <  3:388471/147‑1 (MQ=255)
                                                                                                           aGGGTTTCCTGATATCCATCAAACCGCCCCCAGGCAGCACCATAGCCCATATCCCGCCCGGTATTGTTCCAGTCAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAAca                                       <  5:411171/148‑1 (MQ=255)
                                                                                                              gTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTgg                                                                            <  7:168167/108‑1 (MQ=255)
                                                                                                              gTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTgg                                                                            >  8:168167/1‑108 (MQ=255)
                                                                                                                ttCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAgg                                                                                              >  5:296120/1‑88 (MQ=255)
                                                                                                                ttCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAgg                                                                                              <  6:296120/88‑1 (MQ=255)
                                                                                                                 tCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGTTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGt                                  <  8:352179/147‑1 (MQ=255)
                                                                                                                   cTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTc                                <  6:243828/147‑1 (MQ=255)
                                                                                                                     gATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTACAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCAc                              >  4:77249/1‑147 (MQ=255)
                                                                                                                      aTTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAAttt                                                                                                            >  6:428422/1‑68 (MQ=255)
                                                                                                                      aTTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAAttt                                                                                                            <  5:428422/68‑1 (MQ=255)
                                                                                                                         tCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGATGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTAccc                         >  8:136612/1‑148 (MQ=255)
                                                                                                                             tCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGc                      >  2:64327/1‑147 (MQ=255)
                                                                                                                              cAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCat                                                                                   <  4:386450/85‑1 (MQ=255)
                                                                                                                              cAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCat                                                                                   >  3:386450/1‑85 (MQ=255)
                                                                                                                                            gcagcaCCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACtttt      <  3:77249/148‑1 (MQ=255)
                                                                                                                                             gagcaCCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAAc                                                                                      <  4:241457/66‑1 (MQ=255)
                                                                                                                                             cagcaCCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGATAAc                                                                                      >  3:241457/1‑67 (MQ=255)
                                                                                                                                                caCCATAGCCCATATCCAGCCCTGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACATCAGTTCACTACCCGGCAATTGATACAACTTTTCGtt  <  5:326094/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCAGGATGTTAACCATATGCTGGAGTCGCCCTTCCGTTTCCGGCCCGCTCAGCACTAACAGCAGTTCACTACCCGGCAATTGATACAACTTTTCGTT  >  NC_000913/2516374‑2516665

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