Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	ACB122	3,508,387	A→G	E1333G (GAA→GGA)	rhsA →	putative deoxyribonuclease RhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,508,387	0	A	G	96.5%	88.4 / NA	29	E1333G (GAA→GGA)	rhsA	putative deoxyribonuclease RhsA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/17);  minor base T (0/1);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCAAAGCACCTACTCGTGCCTACTCACCCACAAACCTTGTGCGTATGGCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGAACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAATTTGCTCGCGCTTACTCCGTGGGAGCATGAGGCTGTCGACCAGTTTAGGCATGTAG  >  ACB122/3508256‑3508503                                                                                                                                    |                                                                                                                     ggCCAAAGCACCTACTCGTGCCTACTCACCCACAAACCTTGTGCGTATGGCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGAc                                                                                                                    <  1:206712/134‑1 (MQ=255)             tACTCGTGCCTACTCACCCACAAACCTTGTGCGTATGGCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATg                                                                                                        >  4:168610/1‑134 (MQ=255)             tACTCGTGCCTACTCACCCACAAACCTTGTGCGTATGGCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGa                                                                                                       <  3:14585/135‑1 (MQ=255)               cTCGTGCCTACTCACCCACAAACCTTGTGCGTATGTCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGata                                                                                                     >  7:138048/1‑135 (MQ=255)                          tCACCCACCAACCTTGTGCGGATGGCTGAGCGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGc                                                                                            <  6:57577/133‑1 (MQ=255)                            aCCCACAAACCTTGTGCGTATGGCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGCCTTCATCATAATGATTTTCCTCAGCGAgt                                                                                        >  1:208865/1‑135 (MQ=255)                                      cTTGTGCGTATGGCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCg                                                                              >  8:204231/1‑135 (MQ=255)                                           gCGTATGGCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTtcatca                                                                                                             <  8:60837/99‑1 (MQ=255)                                           gCGTATGGCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTtcatca                                                                                                             >  7:60837/1‑99 (MQ=255)                                           gCGTATGGCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGgtgt                                                                         >  7:197100/1‑135 (MQ=255)                                                   cTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGAtt                                                                 >  8:66811/1‑135 (MQ=255)                                                         ggAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCaa                                                           >  7:77715/1‑135 (MQ=255)                                                             aGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAAttt                                                       <  7:204231/135‑1 (MQ=255)                                                             aGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAAttt                                                       <  7:66811/135‑1 (MQ=255)                                                             aGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAAttt                                                       <  5:194121/135‑1 (MQ=255)                                                             aGCGCCCAAAATGACAGTTGAGGTGATTTCGAGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAAttt                                                       <  7:93664/135‑1 (MQ=255)                                                                gcCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAATTTGCt                                                    >  6:51402/1‑135 (MQ=255)                                                                   cAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAATTTGCTcgc                                                 <  5:51402/135‑1 (MQ=255)                                                                             gTTGAGGTGATTTCGCGCAAGACCTATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAAttt                                                       <  5:55318/119‑1 (MQ=255)                                                                             gTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAAttt                                                       >  6:55318/1‑119 (MQ=255)                                                                             gTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAATTTGCTCGCGCTTACTCCg                                       <  3:39028/135‑1 (MQ=255)                                                                                       attagcgcAAGACCGATAAAATATAAATTAGAGAGCATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAATTTGCTCGCGCTTACTCCGGGGGAGCATg                             <  7:77801/131‑1 (MQ=255)                                                                                               aaGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGTGTGTACATGATTCGTCAAATTTGCTCGCGCTTACTCCGTGGGAGCATGAGGCTGTc                     <  2:37733/135‑1 (MQ=255)                                                                                               aaGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAATTTGCTCGCGCTTACTCCGTGGGAGCATGAGGCTGTc                     <  3:168610/135‑1 (MQ=255)                                                                                               aaGACCGATAAAATATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAATTTGCTCGCGCTTACTCCGTGGGAGCATGAGGCTGTc                     <  8:77715/135‑1 (MQ=255)                                                                                               aaGACCGATAAAATATCAATTAGAGAGTATGAATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTAAATGATTCGTCAAATTTGCTCGCGCTTAATCCGTGGGAGCATGAGGCTGTc                     <  4:168008/135‑1 (MQ=255)                                                                                               aaGACATATAAAATATCAATTAGAGAGTATGCATTATTACTTCATCCTAATGATATTCATCAGCGAGTGGGTGCGGCGGGTGTACATGATTCGTCAAATTTGCTCGCGCTTACTCCGTGGTAGAATGAGGCTGTc                     <  3:111815/135‑1 (MQ=255)                                                                                                         aaaTATCAATTAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAATTTGCTCGCGCTTACTCCGTGGGAGCATGAGGCTGTCGACCAGTTTa           >  7:128960/1‑135 (MQ=255)                                                                                                                  ttAGAGAGTATGCATTAGGACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAATTTGCTCGCGCTTACTCCGTGGGAGCATGAGGCTGTCGACCAGTTTAGGCATGTAg  <  5:99551/135‑1 (MQ=255)                                                                                                                                    |                                                                                                                     GGCCAAAGCACCTACTCGTGCCTACTCACCCACAAACCTTGTGCGTATGGCTGAGGGGAAAGCGCCCAAAATGACAGTTGAGGTGATTTCGCGCAAGACCGATAAAATATCAATTAGAGAGTATGCATTAGAACTTCATCATAATGATATTCCTCAGCGAGTGGGTGGGGCGGGTGTACATGATTCGTCAAATTTGCTCGCGCTTACTCCGTGGGAGCATGAGGCTGTCGACCAGTTTAGGCATGTAG  >  ACB122/3508256‑3508503

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑