Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	612,076	A→G	12.8%	T130A (ACA→GCA)	sdrC →	Ser‑Asp rich fibrinogen/bone sialoprotein‑binding protein SdrC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	612,076	0	A	G	12.8%	78.9 / 3.0	33	T130A (ACA→GCA)	sdrC	Ser‑Asp rich fibrinogen/bone sialoprotein‑binding protein SdrC
Reads supporting (aligned to +/- strand):  ref base A (13/15);  new base G (2/2);  total (15/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 = 9.65e-01

ACGCAAACTGCAACTGCAGATCAGCCTAAAGTGACAATGAGTGATAGTGCAACAGTTAAAGAAACTAGTAGTAACATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAATCGCATGGCAGTGA  >  CP000730/611944‑612208                                                                                                                                     |                                                                                                                                     aCGCAAACTGCAACTGCAGATCAGCCTAAAGTGACAATGAGTGATAGTGCAACAGTTAAAGAAACTAGTAGTAACATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAAt                                                                                                                              >  2:251502/1‑141 (MQ=255)   gCAAACTGCAACTGCAGATCAGCCTAAAGTGACAATGAGTGATAGTGCAACAGTTAAAGAAACTAGTAGTAACATGCAATCACCACAAAACGCTCCAGTTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGa                                                                                                                            <  1:118804/141‑1 (MQ=255)    cAAACTGCAACTGCAGATCAGCCTAAAGTGACAATGAGTGATAGTGCAACAGTTAAAGAAACTAGTAGTAACATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGAAATGTAACAACAAATGaa                                                                                                                           >  2:128365/1‑140 (MQ=255)          gcaactgcaGATCAGCCTAAAGTGACAATGAGTGATAGTGCAACAGTTAAAGAAACTAGTAGTAACATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAatca                                                                                                                     <  2:248082/141‑1 (MQ=255)            aactgcaGATCAGCCTAAAGTGACAATGAGTGATAGTGCAACAGTTAAAGAAACTAGTAGTAACATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAatcata                                                                                                                   >  2:101027/1‑140 (MQ=255)                    aTCAGCCTAAAGTGACAATGAGTGATAGTGCAACAGTTAAAGAAACTAGTAGTAACATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAGCAACAAATGATAAATCATCAACTACat                                                                                                           >  1:148641/1‑141 (MQ=255)                    aTCAGCCTAAAGTGACAATGAGTGATAGTGCAACAGTTAAAGAAACTAGTAGTAACATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAGCAACAAATGAAAACTCATCAACTCCat                                                                                                           >  1:228599/1‑141 (MQ=255)                                    aaTGAGTGATAGTGCAACAGTTAAAGAAACTAGTAGTAACATGCAATAACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGa                                                                                           >  2:10270/1‑141 (MQ=255)                                         gtcatatggCAACAGTTAAATAAACTATTAGTAATATGCTTTCTCCACTTCCCCCTACATATAATCTATCTACTACAAACACTAGCAATTTATCTACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATaaaa                                                                                      <  1:237408/133‑1 (MQ=255)                                                         taaataaaCTAGTAGTAACATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGc                                                                      >  1:184021/1‑141 (MQ=255)                                                            aGAAACTAGTAGTAACATGGAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCaaa                                                                   <  2:21872/141‑1 (MQ=255)                                                             gAAACTAGTAGTAACATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAAATAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAAATGATAAAAGTAATTTAAAACaacaaata                                                                  >  2:127525/1‑135 (MQ=255)                                                                     tagtaACACGCAAGCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACGAAATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTc                                                          <  1:68469/141‑1 (MQ=255)                                                                          aCATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTa                                                     <  1:223966/141‑1 (MQ=255)                                                                               cAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCt                                                >  1:255715/1‑141 (MQ=255)                                                                               cAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCt                                                >  1:161697/1‑141 (MQ=255)                                                                                 accaccacCAAACGTTACAGCTAATCAATCTACTACAACAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTaa                                              <  1:268819/139‑1 (MQ=255)                                                                                        caAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACg                                       <  1:251502/141‑1 (MQ=255)                                                                                        caAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACg                                       <  1:39259/141‑1 (MQ=255)                                                                                        caAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAAATACATATAGTAATGAAACTGATAAAAGTAATTTAACAAAAGCAAAAGAAGTTTCAACTACACCTAAaacaaca                                       >  2:64669/1‑140 (MQ=255)                                                                                             cGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTAt                                  <  1:186026/141‑1 (MQ=255)                                                                                                tACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTaa                               <  1:77165/141‑1 (MQ=255)                                                                                                            atcGACAACAAAAACTGGCAATGTGACAACAAATGATAAATCAGCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACttt                   >  2:158214/1‑141 (MQ=255)                                                                                                             actactaCAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTa                  <  2:61516/140‑1 (MQ=255)                                                                                                                actacAAAAACTAGCAATTTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAAt               <  1:10270/141‑1 (MQ=255)                                                                                                                actacAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAAt               >  1:227075/1‑141 (MQ=255)                                                                                                                 ctacAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAATc              <  1:258594/141‑1 (MQ=255)                                                                                                                 ctacAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAATc              <  1:258570/141‑1 (MQ=255)                                                                                                                     aaaaaCTAGCAATGTAGCAACAAATGATAAATCATCAACTACATATAGTAATGTAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAATCGCAt          <  2:148641/141‑1 (MQ=255)                                                                                                                     aaaaaCTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAATCGCAt          >  1:157730/1‑141 (MQ=255)                                                                                                                     aaaaaCTAGCAATGGAGCAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATGTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAATCGCAt          <  2:228599/141‑1 (MQ=255)                                                                                                                            aGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAATCGCATGGCAGTg   <  1:187534/141‑1 (MQ=255)                                                                                                                             gCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAATCGCATGGCAGTGa  >  2:59506/1‑141 (MQ=255)                                                                                                                                     |                                                                                                                                     ACGCAAACTGCAACTGCAGATCAGCCTAAAGTGACAATGAGTGATAGTGCAACAGTTAAAGAAACTAGTAGTAACATGCAATCACCACAAAACGCTACAGCTAATCAATCTACTACAAAAACTAGCAATGTAACAACAAATGATAAATCATCAACTACATATAGTAATGAAACTGATAAAAGTAATTTAACACAAGCAAAAGATGTTTCAACTACACCTAAAACAACGACTATTAAACCAAGAACTTTAAATCGCATGGCAGTGA  >  CP000730/611944‑612208

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