Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	3,951,161	A→T	intergenic (‑125/+70)	PP_3486 ← / ← PP_3487	putative cytochrome c/conserved exported protein of unknown function

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	3,951,161	0	A	T	95.8%	73.6 / NA	24	intergenic (‑125/+70)	PP_3486/PP_3487	putative cytochrome c/conserved exported protein of unknown function
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (13/10);  minor base C (1/0);  total (14/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

CAACTGGCGCAGGGCAGCACGCTGGCGGGGACTGTGATATCGCTGGCGCCGACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCAT  >  AE015451/3950873‑3951448                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                cAACTGGCGCAGGGCAGCACGCTGGCGGGGACTGTGATATCGCTGGCGCCGACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAg                                                                                                                                                                                                                                                                                     <  2:60025/301‑1 (MQ=255)    gatcggcAGGGCAGCACGCTGGCGGGGACTGTGATATCGCTGGCGCCGACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGcc                                                                                                                                                                                                                                                                                  <  1:282681/296‑1 (MQ=255)       gcgcAGGGCAGAACGCTGGCGGGGACTGTGATATCGCTGGCGCCGACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCt                                                                                                                                                                                                                                                                               <  1:187622/301‑1 (MQ=255)                 ggcgcgtgggggggACTGTGATATCGCTGGGGCCGACCCGGGCTTGCAATTGGGCGTCGCTACTGGAGACGCCTTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCAAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGGGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAg                                                                                                                                                                                                                                                                     <  2:384625/291‑1 (MQ=255)                      cTGGCGGGGACTGTCGTATCGCTGGCGGCGCCCCGGGGTTGGAATTGTGCGTCGCTGCTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCAGGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGAt                                                                                                                                                                                                                                                                <  2:824169/301‑1 (MQ=255)                                       ctgccgggCGCCGATCCGGGTTGGAAATTGGGCGTGGCTACTGGAGACGCCCTCGGGGTAAAGGCGTTTGCCGGCTTGTTCGTGGTCGGGCAGGTCCAGGGCAAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATAGCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTc                                                                                                                                                                                                                                               <  2:493244/295‑1 (MQ=255)                                                   gACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGccc                                                                                                                                                                                                                                   >  1:145391/1‑301 (MQ=255)                                                                  ttGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCCCAGAAAGAGGGGAGCGCTGCGCGCCCCTCTCGCGACGCCATGCCCCTCCCACAAAAGCGGCTGCCACAGAGCCCAATTCCCCATGCGGccccacaacccccccgcc                                                                                                                                                                                                                    >  2:483874/1‑286 (MQ=255)                                                                               cTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGCGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCACGCGCAGAAAGAGTGGAGCGCGGCGCGCGCCTTTCGCGCCGCAAGACCGCGCCCACCAACGCGGACGCCAACGAGGCCAATGTCCCAGGCGGCCCCTGCGCCTCaccgaccccccaccggccg                                                                                                                                                                                                       >  2:655931/1‑289 (MQ=255)                                                                                          cTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGAGGCCCACGATACTCACCGACCGGCCAACCGCTGGGtacacagcc                                                                                                                                                                                            >  1:879460/1‑293 (MQ=255)                                                                                               ggTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCACCGGCCGGGTCAACAGCTGCATg                                                                                                                                                                                       >  1:206164/1‑301 (MQ=255)                                                                                                   aaGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATgcgcc                                                                                                                                                                                   >  1:157395/1‑300 (MQ=255)                                                                                                                             gggtcggttaaaggggcagggcTGAGGACGAAGTGGCGATGCAAAGCAGGCGGGTAAGCCAAATCGCAAGCCCGGTGGGAGGGGCATTGCGGGGCGAAAGGGCGGCGCAGCAGCCCCGGCAGTGGCAGCTCAGGCGCAGAAAGAGGGGAGCGCTGCGCGCGCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATcc                                                                                                                                                         <  2:879152/290‑1 (MQ=255)                                                                                                                                                        tcgagtatgccatgccaCGCACGGTGGGTAGCCATATCGCATGCCCTGTGGGATCTGCCGCTCGTCGCGAAAGGGGTGCGCAGCCGCCACGGCAGTGGCAGCTCAAGCGCAGAAAGAGTTGATCGCTGCTCGCTCCTTGCCCTCCGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCTGGCGAATCTTTCTCCAGCGTTTACCt                                                                                                                              <  1:327760/294‑1 (MQ=255)                                                                                                                                                           aaGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCa                                                                                                                           >  1:843417/1‑301 (MQ=255)                                                                                                                                                                                        caaagcgaaaGCCCGGTGGGAGCGGCCTGGCGGCGAGAAAGGGCTGCGCGGCAGCCCCGGAAGGGGGAGCTCAAGCGCAGAAAGAGGGGAGCGCGGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAAtg                                                                                              <  2:145340/293‑1 (MQ=255)                                                                                                                                                                                                                      ggTCGCGAAAGGTGTGCGCAGCAGCCCCGGCGGTGGCAGCTCAAGCGCAGAAAGCGGGGAGAGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCtat                                                                <  1:183527/300‑1 (MQ=255)                                                                                                                                                                                                                            gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGTCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACgccc                                                          >  2:271530/1‑298 (MQ=255)                                                                                                                                                                                                                            gAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCCAATGCTGCTCGGGACCCGGCAGGAACTGCCCTATACGGcc                                                          >  2:737301/1‑301 (MQ=255)                                                                                                                                                                                                                                        cagcagCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTc                                              >  1:456706/1‑301 (MQ=255)                                                                                                                                                                                                                                                             cTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCTGCAATCAGCTGCGTCGAATGCTGTTCGGGACTCGTCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGcgtttgcgccggg                         >  2:835460/1‑289 (MQ=255)                                                                                                                                                                                                                                                                        aaaGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGcc              >  1:337311/1‑301 (MQ=255)                                                                                                                                                                                                                                                                          agagGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAg            >  1:55391/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                    gcTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGGCGGTCAGGAACAGCCAGTTCGGCCCac  >  2:71132/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                CAACTGGCGCAGGGCAGCACGCTGGCGGGGACTGTGATATCGCTGGCGCCGACCCGGGCTTGCAATTGAGCGTCGCTACTGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAGGAACTGCCCTATACGGCCAAGCACCAGGTCGACATTGGCCTTGTCGCCGGTCAGGAACAGCCAGTTCGGCCCAT  >  AE015451/3950873‑3951448

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