Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	981,030	T→G	28.8%	N155T (AAC→ACC)	U712_04790 ←	Alkaline phosphatase 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	981,030	0	T	G	28.8%	46.8 / 5.6	40	N155T (AAC→ACC)	U712_04790	Alkaline phosphatase 4
Reads supporting (aligned to +/- strand):  ref base T (5/23);  new base G (10/2);  total (15/25)
Fisher's exact test for biased strand distribution p-value = 1.67e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAAGATTTTCCGCCGCCGAGCAGCACGTCTATTTTATGTTTGCCTTTTATCTTGTCATCCATATAGCTGTTGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCATTGTTATATGTCTTAACGCCTGTCGC  >  CP006881/980893‑981157                                                                                                                                          |                                                                                                                                aaTAAGATTTTCCGCCGCCGAGCAGCACGTCTATTTTATGTTTGCCTTTTATCTTGTCATCCATATAGCTGTTGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGtt                                                                                                                                <  7:180850/139‑1 (MQ=255)       aTTTTCCGCCGCCGAGCAGCACGTCTATTTTATGTTTGCCTTTTATCTTGTCATCCATATAGCTGTTGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTTGTTTAttt                                                                                                                           >  8:238842/1‑138 (MQ=255)                                  ttttATGTTTGCCTTTTATCTTGTCATCCATATAGCTGTTGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTTTTTTGCGACGTGGCGACCAGCCCGGTTTActt                                                                                               >  6:9543/1‑139 (MQ=255)                                            gCCTTTTATCTTGTCATCCATATAGCTGTTGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGGTTATTTAAGACGTGGCGACAATCCCTGTTGGCGTGCCTTGCTGt                                                                                     >  5:164955/1‑139 (MQ=255)                                                   aTCTTGTCATCCTTATTGCTGTTGGCGATTTGGTCCATGTTTATCCGTGATTCAGTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGcc                                                                               <  2:147634/138‑1 (MQ=255)                                                      agacagtCCATATAGCTGTTGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCtt                                                                           <  5:202942/133‑1 (MQ=255)                                                       tGTCATCCATATAGCTGTTGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGGTTATTTTTAACGGGGCGCCCAGCCCCGGTGGCGTGGCTGGGTGTTTgggcgcgt                                                                           >  6:266968/1‑132 (MQ=255)                                                            tCCATATAGCTGTTGGCTATTTGGTCCTGGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTa                                                                     <  2:78174/139‑1 (MQ=255)                                                                   aGCTGTTGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGGTTATTTCATACTGTTCGAAAAGCCCTGTTTGCGTTGCTTGGTGTGTTGCCGCTTGAAGTAAAGCttt                                                              >  7:326939/1‑139 (MQ=255)                                                                        ttGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGt                                                                                                                                 >  4:280115/1‑67 (MQ=255)                                                                        ttGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGt                                                                                                                                 <  3:280115/67‑1 (MQ=255)                                                                        ttGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGGTTATTTAAGACGTGGGGACAAGCCCAGGTTGCTTGACTTTGTGGTTGGGCGGTTGAAGTAAAAATTTAAcat                                                         >  6:332202/1‑137 (MQ=255)                                                                        ttGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGGTTATTTAAGACGTGGCGCCAAGCCCGGTTGAGATGCCTTGGGGGTTGGGCTGTTGAAGTTCAGAATTGACtt                                                         >  1:228023/1‑139 (MQ=255)                                                                        ttGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGGTTATTTAAGACGTGGCGCCAAGCCCAGTTGGCTTTCCTTGGTGTTTGGCCTCGTGAAGTTAATAATTCACtt                                                         >  4:152915/1‑139 (MQ=255)                                                                          ggCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGGTTATTTGATACGAGGGGGCAACCCCAGTTGACATGGCGTGGTGTTTGGGGGCTTAAAGGAAAGAATTGACtttt                                                       >  3:257578/1‑139 (MQ=255)                                                                          ggCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGGTTAATTAAGACGTGGCGACAAGCCCGGGTGACTTGGCTTGCTGTTTGGGCGAGTAAAGGACAGGTTTCaagttt                                                       >  8:322228/1‑134 (MQ=255)                                                                          ggCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGGGGTTTATTTTAGACGTGGCGACAAGCCCTGTTGACTTGCCGTGCTGTTTGGGCGCTTCAAGTACAGATTTGACtttt                                                       >  5:20312/1‑139 (MQ=255)                                                                          ggCGATTTGGTCCATGTTTTTCCGTGATTCATTCTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTTAATTAACACCGGGCCACACCCCCCGTTCACTTTCCTTGCTGTTTGGGCTCATGAAGGTAAGATTTCACgttt                                                       >  3:155345/1‑139 (MQ=255)                                                                             gATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACttttttt                                                    <  5:92948/139‑1 (MQ=255)                                                                                        tGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATc                                         <  7:322228/139‑1 (MQ=255)                                                                                                cGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCaa                                 <  2:189533/139‑1 (MQ=255)                                                                                                    aTTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTg                              <  6:164955/138‑1 (MQ=255)                                                                                                     ttCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCa                            <  5:332202/139‑1 (MQ=255)                                                                                                     ttCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCa                            <  7:181167/139‑1 (MQ=255)                                                                                                     ttCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCa                            <  1:180636/139‑1 (MQ=255)                                                                                                     ttCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCa                            <  5:266968/139‑1 (MQ=255)                                                                                                     ttCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCa                            <  5:247751/139‑1 (MQ=255)                                                                                                     ttCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCa                            <  1:173691/139‑1 (MQ=255)                                                                                                     ttCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCa                            <  3:152915/139‑1 (MQ=255)                                                                                                     ttCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCa                            <  7:84828/139‑1 (MQ=255)                                                                                                     ttCATTGTGGGCGCCATATGCGGCGGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCa                            <  5:112990/139‑1 (MQ=255)                                                                                                     acagttGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCa                            <  7:216899/135‑1 (MQ=255)                                                                                                       cATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCAt                           <  5:65098/138‑1 (MQ=255)                                                                                                            tgGGCGCCATATGCGGCTGGAGTGGCGGGGGTTATTTTAGACCTGGCGCCAAGCCCAGGTTACATGCCTTTCAGTTTTGCCAGTTAAAGTAAAGCTTTAAAATTTTTTACGTTTTTTTAGAAGCCAATTTCATTTTta                      >  6:247751/1‑138 (MQ=255)                                                                                                             gggCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCATTGTtata                    <  1:166182/139‑1 (MQ=255)                                                                                                             gggCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCATTGTtata                    <  3:124350/139‑1 (MQ=255)                                                                                                             ggcggcCATATGCGCCGGTAGGGCCGGGGTTAATTTCAGCCGGGGCTACATGCCTGGTGTACTTGCCTTGCGGTTTGGCCTCTCCACGTCCTGATTTCACTTTTTTTCCGTTTTAATCGCGGCCAATGGCATGGTtata                    <  2:228023/135‑1 (MQ=255)                                                                                                               gcgcCATATGCGGCTGGAGTGGCGTGGTTTATTTAAGACGGGGCGGCCCGCCCGGGTTGCTTGCCGTGGGGGTTTGGCGCGTCAAGTTAAGCTTTCACATTTTTTGCGGTTTTATGGACGCCAATTGGATTGTTttatt                  >  2:180636/1‑138 (MQ=255)                                                                                                                    atatGCGGCTGGAGTGGCGTGGGTAATTTGAGACGATGCGCCAAGCCCAGTTGACATTGCATGGTGTTTGGCCACTTCAACTACAGAATTAACGTTTTTTTCGTTTTTATAGACCCCAATTGCATTTTTTTATGTcag              >  8:216899/1‑136 (MQ=255)                                                                                                                         cGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCATTGTTATATGTCTTAACGc         <  4:257578/138‑1 (MQ=255)                                                                                                                          ggCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCATTGTTATATGTCTTAACGCCt       <  1:144808/139‑1 (MQ=255)                                                                                                                          ggCTGGAGTGGCGTGGGTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCATTGTTATATGTCTTAACGCCt       <  4:146106/139‑1 (MQ=255)                                                                                                                               gAGGGGCGTGGGTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCATTGTTATATGTCTTAACGCCTGTCGc  <  8:211006/139‑1 (MQ=255)                                                                                                                                          |                                                                                                                                AATAAGATTTTCCGCCGCCGAGCAGCACGTCTATTTTATGTTTGCCTTTTATCTTGTCATCCATATAGCTGTTGGCGATTTGGTCCATGTTTTTCCGTGATTCATTGTGGGCGCCATATGCGGCTGGAGTGGCGTGGTTAATTTCAGACGTGGCGACAAGCCCTGTTGACTTGCCTTGCTGTTTGGCCTCTTCAAGTACAGATTTCACTTTTTTTCCGTTTTTATCGACGCCAATTGCATTGTTATATGTCTTAACGCCTGTCGC  >  CP006881/980893‑981157

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: ‑