Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,623,117	A→G	I245T (ATC→ACC)	dgcZ ←	diguanylate cyclase DgcZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,623,117	0	A	G	97.1%	112.7 / ‑4.4	35	I245T (ATC→ACC)	dgcZ	diguanylate cyclase DgcZ
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (15/19);  total (16/19)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTTGTTCGTCAATAAACATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTT  >  NC_000913/1622980‑1623251                                                                                                                                          |                                                                                                                                       tttGTTCGTCAATAAACATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGgta                                                                                                                                      >  2:390270/1‑140 (MQ=255) tttGTTCGTCAATAAACATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGgta                                                                                                                                      >  2:239115/1‑140 (MQ=255)     ttCGTCAATAAACATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGc                                                                                                                                   <  2:597550/139‑1 (MQ=255)     ttCGTCAATAAACATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGc                                                                                                                                   >  1:597550/1‑139 (MQ=255)       cGTCAATAAACATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCt                                                                                                                                <  2:21812/140‑1 (MQ=255)             tAAACATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCtt                                                                                                                               >  2:783829/1‑135 (MQ=255)             tAAACATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCtt                                                                                                                               <  1:783829/135‑1 (MQ=255)             tAAACATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGaa                                                                                                                          <  2:55041/140‑1 (MQ=255)              aaaCATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAAt                                                                                                                         <  1:184410/140‑1 (MQ=255)                 cATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAAtgtg                                                                                                                      <  1:390270/140‑1 (MQ=255)                         gATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCa                                                                                                              <  1:701687/140‑1 (MQ=255)                             tcttcCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGt                                                                                                          <  2:173224/140‑1 (MQ=255)                                            aCCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGc                                                                                           <  2:438273/140‑1 (MQ=255)                                                   tACATTGCTCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCt                                                                                    >  2:239997/1‑140 (MQ=255)                                                        tGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGacac                                                                               >  1:692181/1‑140 (MQ=255)                                                          cccGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCt                                                                             >  1:575671/1‑140 (MQ=255)                                                           ccGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTg                                                                            <  2:484879/140‑1 (MQ=255)                                                            cGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGacac                                                                               >  1:367114/1‑136 (MQ=255)                                                            cGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGacac                                                                               <  2:367114/136‑1 (MQ=255)                                                            cGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGc                                                                           <  2:692181/140‑1 (MQ=255)                                                            cGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGc                                                                           >  2:628582/1‑140 (MQ=255)                                                                   ctctTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGTCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACAt                                                                    <  2:276261/140‑1 (MQ=255)                                                                                cacaTCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTcatcat                                                       >  2:760380/1‑140 (MQ=255)                                                                                cacaTCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCAT‑GTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCAtt                                                      >  2:430076/1‑140 (MQ=255)                                                                                  caTCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTa                                                     >  2:345089/1‑140 (MQ=255)                                                                                    tCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGc                                                   <  1:239997/140‑1 (MQ=255)                                                                                             tctTGAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGa                                           <  2:58905/139‑1 (MQ=255)                                                                                             tctTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGa                                           >  1:58905/1‑139 (MQ=255)                                                                                               ttCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACaa                                        <  1:239115/140‑1 (MQ=255)                                                                                                         gcgcGACTCACACCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATaaa                              >  2:763176/1‑140 (MQ=255)                                                                                                                  acacCTGCTGTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGcc                     <  2:575671/140‑1 (MQ=255)                                                                                                                         ctgTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAg              >  2:244618/1‑140 (MQ=255)                                                                                                                         ctgTCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAg              >  1:44587/1‑140 (MQ=255)                                                                                                                            tCACGGTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGa           >  1:311427/1‑140 (MQ=255)                                                                                                                                 gTAATGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAAc      >  1:163759/1‑140 (MQ=255)                                                                                                                                     tGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGtt  <  2:44587/140‑1 (MQ=255)                                                                                                                                     tGTTGGTATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGtt  <  1:760380/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       TTTGTTCGTCAATAAACATGCAGCGATTTCTTCCGGTTTGCTTACCCTCATACATTGCCCGGTCCGCTCTTCCAATGACCACATCCAGAGGCTCTTCAGGAAATGCGCGACTCACACCTGCTGTCACGGTAATGTTGATATGCCCTTCAGAATGTGTGATGGCATGGTTATCGACTAACTGGCAAATTCTGACACCTGCACGACATGCTTCTTCATCATTAGCCGCTTTGACAATAATGATAAATTCTTCGCCCCCGTAGCGATAAACCGTT  >  NC_000913/1622980‑1623251

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