Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,683,630 C→T *663* (TAG→TAA)  yhjK ← cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,683,6300CT97.4% 119.5 / ‑6.5 39*663* (TAG→TAA) yhjKcyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/18);  minor base G (0/1);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCA  >  NC_000913/3683517‑3683762
                                                                                                                 |                                                                                                                                    
aGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCtt                                                                                                             >  1:268787/1‑139 (MQ=255)
aGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCtt                                                                                                             >  1:444743/1‑139 (MQ=255)
    aacaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                         <  1:747844/139‑1 (MQ=255)
      caTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                       >  1:106513/1‑139 (MQ=255)
              aaaTGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATa                                                                                               >  1:406296/1‑139 (MQ=255)
                aTGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                            >  2:362865/1‑108 (MQ=255)
                aTGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                            <  1:362865/108‑1 (MQ=255)
                aTGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAgg                                                                                             >  1:547160/1‑139 (MQ=255)
                 tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                            >  2:798639/1‑107 (MQ=255)
                 tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                            <  1:798639/107‑1 (MQ=255)
                 tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAgg                                                                                                                       >  1:62627/1‑112 (MQ=255)
                 tGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAgg                                                                                                                       <  2:62627/112‑1 (MQ=255)
                   tACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAg                                                                                          >  1:429094/1‑139 (MQ=255)
                        ggATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                            <  1:10831/100‑1 (MQ=255)
                        ggATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTttctt                                                                                                                            >  2:10831/1‑100 (MQ=255)
                            cATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGc                                                                                 <  2:627773/139‑1 (MQ=255)
                                   gAAACAAAAAGGTTAAAACGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGGTACTTTTCTTCCAGGTAAc                                                                                                                   <  2:555258/98‑1 (MQ=255)
                                   gAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                                          <  2:444743/139‑1 (MQ=255)
                                   gAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                                          <  2:268787/139‑1 (MQ=255)
                                   gAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAAACAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAAc                                                                                                                   >  1:555258/1‑98 (MQ=255)
                                    aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTcc                                                                                                                          >  2:647506/1‑90 (MQ=255)
                                    aaaCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTcc                                                                                                                          <  1:647506/90‑1 (MQ=255)
                                      aCAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAg                                                                       >  2:706981/1‑139 (MQ=255)
                                        aaaaaGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGcc                                                                     <  2:106513/139‑1 (MQ=255)
                                             ggTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGTCCTggg                                                                <  2:429094/139‑1 (MQ=255)
                                                                     cacaAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTgcca                                        >  2:71138/1‑139 (MQ=255)
                                                                          aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGccc                                                                    >  2:678622/1‑106 (MQ=255)
                                                                          aCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGccc                                                                    <  1:678622/106‑1 (MQ=255)
                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTggg                                                                >  2:207925/1‑100 (MQ=255)
                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTggg                                                                <  1:207925/100‑1 (MQ=255)
                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTcgcg                              >  2:104370/1‑134 (MQ=255)
                                                                                    tttAAAGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTcgcg                              <  1:104370/134‑1 (MQ=255)
                                                                                       aaaGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCtt                      <  2:141142/139‑1 (MQ=255)
                                                                                         aGTTTTGTAATCAGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCa                    <  2:547160/139‑1 (MQ=255)
                                                                                                     aGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCa             >  2:648076/1‑134 (MQ=255)
                                                                                                     aGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCa             <  1:648076/134‑1 (MQ=255)
                                                                                                     aGTTTGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCt        >  1:72226/1‑139 (MQ=255)
                                                                                                       tttGGGGTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTc      >  1:305773/1‑139 (MQ=255)
                                                                                                           gggTAGTTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCa  <  1:71138/139‑1 (MQ=255)
                                                                                                                 |                                                                                                                                    
AGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCA  >  NC_000913/3683517‑3683762

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