Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,474,639	C→G	A9129P (GCA→CCA)	ebh ←	extracellular matrix binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,474,639	0	C	G	89.7%	87.1 / 3.9	29	A9129P (GCA→CCA)	ebh	extracellular matrix binding protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base G (13/13);  total (15/14)
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.80e-01

GCGACCTTTAAATCCATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGCTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAAC  >  CP000730/1474502‑1474759                                                                                                                                          |                                                                                                                         gCGACCTTTAAATCCATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTg                                                                                                                       >  1:53748/1‑141 (MQ=255)        ttAAATCCATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACt                                                                                                                >  2:214786/1‑141 (MQ=255)               cATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAATTGCTACACGTTCATCACCATTTAATGCTTGTGTTGCATTTAAAATATTTTGAAGTGCTTGATCCACTTGTTGCt                                                                                                         >  1:181181/1‑141 (MQ=255)                  cTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTAt                                                                                                      <  2:293411/141‑1 (MQ=255)                  cTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGCTTGATCCACTTGTTGCTTAt                                                                                                      <  1:225569/141‑1 (MQ=255)                       tgAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCa                                                                                                 <  2:197613/141‑1 (MQ=255)                       tgAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCa                                                                                                 <  2:118370/141‑1 (MQ=255)                         aGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCAtt                                                                                               >  2:70322/1‑141 (MQ=255)                         aGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCAtt                                                                                               <  1:250904/141‑1 (MQ=255)                               gTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGttt                                                                                         >  1:101555/1‑141 (MQ=255)                                              ttGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTtcagc                                                                          >  1:218213/1‑141 (MQ=255)                                                ggTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGTAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTtcagcat                                                                        >  2:222672/1‑141 (MQ=255)                                                 ttCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTtcagcatc                                                                       <  1:104650/140‑1 (MQ=255)                                                   cGATGTCATGTTTACCCTTTGTTTTAGCAAGTGCTACACGCTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTtcagcatcag                                                                     <  2:269381/141‑1 (MQ=255)                                                    gATGTCATGTTTACCATTTTTTTTAGCAAGTGCTACACGTTCATCACCATTTAATTCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTtcagcatcagc                                                                    <  1:195922/141‑1 (MQ=255)                                                     cattaatatttCCCATTTGTTTTAATCAGTACTACACGTTCATCACCATTAAATGCTTGTTTTGTATTTACAATATTTTCAAGTCGATGATCCACTTGTTGCTTATTCTCATTTTGTTTCAATACATTTtcagcatcagca                                                                   <  2:282456/133‑1 (MQ=255)                                                     aTGTCATCTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATATTGTTTCAATACATTTtcagcatcagca                                                                   >  1:269381/1‑141 (MQ=255)                                                       gTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAAt                                                                 <  2:245877/141‑1 (MQ=255)                                                          aTGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGc                                                              <  2:16301/141‑1 (MQ=255)                                                             tttACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCttt                                                           >  2:185086/1‑140 (MQ=255)                                                             tttACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGCTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTg                                                           >  1:227011/1‑141 (MQ=255)                                                                                     tACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGa                                   >  2:179006/1‑141 (MQ=255)                                                                                           ttCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGc                             <  1:36602/141‑1 (MQ=255)                                                                                               tcaCCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTtgt                         >  1:206919/1‑141 (MQ=255)                                                                                                caCCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTtgtg                        <  1:69879/141‑1 (MQ=255)                                                                                                    aTTTAAAGCTTGATTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTGGTTGCTTATTCGCATTTTGTTTCAATAAATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAAtcttctaatttgta                    >  2:177447/1‑133 (MQ=255)                                                                                                           gCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACa             >  1:86617/1‑141 (MQ=255)                                                                                                                 tttGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTTCTTGTGTGTAATCTACAGTTTGc       <  2:118573/141‑1 (MQ=255)                                                                                                                      aTTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAAc  >  2:223028/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                         GCGACCTTTAAATCCATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGCTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAAC  >  CP000730/1474502‑1474759

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