Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	1,419,549	+C	intergenic (‑38/+98)	PP_1242 ← / ← PP_1244	conserved protein of unknown function/conserved protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,419,548	1	.	C	94.1%	52.4 / ‑4.7	17	intergenic (‑37/+99)	PP_1242/PP_1244	conserved protein of unknown function/conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (7/9);  total (7/10)
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

CCCAACTCGGCATCATCCAGCAGCGCCTTGGCCATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GG‑CTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCACGGATTGGCGGGCCACGGACAGGGCGCAGAGTACGTACTCT  >  AE015451/1419255‑1419833                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                              cccAACTCGGCATCATCCAGCAGCGCCTTGGCCATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTc                                                                                                                                                                                                                                                                                         <  1:198563/301‑1 (MQ=255)                cccAGCAGCGCCTTGGCCATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCt                                                                                                                                                                                                                                                                          <  1:463390/300‑1 (MQ=255)                                         ccAGGTAATGCCAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGa                                                                                                                                                                                                                                                 <  2:476463/300‑1 (MQ=255)                                          cacgTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGAc                                                                                                                                                                                                                                                <  2:299351/298‑1 (MQ=255)                                                 tGCGAGGCCCAGATCAGTGTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgtg                                                                                                                                                                                                                                         <  1:237453/301‑1 (MQ=255)                                                                 gTTTCTCGTCGCCGTCCATCAGGCCGGTGATGAGCCGGTCGCGCACTTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGTGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGc                                                                                                                                                                                                                         <  1:137211/301‑1 (MQ=255)                                                                     gggggcgcggTCCATAAGGCCGGTGATGAGCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTGGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGc                                                                                                                                                                                                                     <  2:754221/292‑1 (MQ=255)                                                                                 cATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGCCGTGCCCGCGAag                                                                                                                                                                                                         >  1:484081/1‑301 (MQ=255)                                                                                                                                   cctgatggagggggTCCTGGCAGGGGATGCCTGGTTCGATGCGGAACAGCGGGTGAGTGGTGTTTTGGCCTTGGGAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCAGTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCGCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGAcc                                                                                                                                                       <  1:697527/288‑1 (MQ=255)                                                                                                                                    cTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGGCCCTGTGGGAGCGAGCGTGCCAGGGAAGGGGCGTGACCTAGTAATGCCAggcgggggggtatctggggaagtgcgctcc                                                                                                                                                      >  2:463030/1‑273 (MQ=255)                                                                                                                                                           ggATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGCGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGGGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCg                                                                                                                               >  1:119957/1‑301 (MQ=255)                                                                                                                                                                                                gttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACCGGTACTGTGCATTGCGTGCCGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCGCGCGAAGAGGCCTGCCCTGCTACTGCAAGGTGCGCGGCTCTGTGGGTCAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGAt                                                                                          >  1:433448/1‑301 (MQ=255)                                                                                                                                                                                                   ttcgctgtggtcggggtgggttttggcgtcTGTGGTGGGGTCTTGTGTGTTGGCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGG‑CTCTTCGCGGGCACGCCTGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCAc                                                                                       <  2:368959/272‑1 (MQ=255)                                                                                                                                                                                                        cTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGcat                                                                                  >  1:716407/1‑299 (MQ=255)                                                                                                                                                                                                                                                tgtgttgtCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCAc                                          >  1:45358/1‑301 (MQ=255)                                                                                                                                                                                                                                                                     tttgtttCCCTGGATTCTGGGTTGCCTGGGCT‑GGCCTCTTCTCGGGCACGCCCGCTCCTACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGTAGCGGGCGTTCTCGCCAGGAGGCCGGACCTGCCAATGCCATGTGTGAGGCTCTGTGGGTAACTGCGACCGCATCTGTAACTGCATGAACGCCGACTCGACCAGCCCCTGTACCTCTTCCATTCCATCGATCACGGCCTGCATCACGACCGAGTTGTGGGCGTAGGGTTTGAGtttccctgtttgccggtgcccggacag                     >  1:201953/1‑276 (MQ=255)                                                                                                                                                                                                                                                                                        ggTTGCCTGGGCT‑GGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGGAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCACGGATTGGCGGGCCACGGACAGGGCGCAGAGTACGTActct  <  2:45334/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                              CCCAACTCGGCATCATCCAGCAGCGCCTTGGCCATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCT‑GG‑CTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCACGGATTGGCGGGCCACGGACAGGGCGCAGAGTACGTACTCT  >  AE015451/1419255‑1419833

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑