Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	3,951,159	T→C	100%	intergenic (‑123/+72)	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,159	0	T	C	100.0%	126.9 / NA	36	intergenic (‑123/+72)	PP_3486/PP_3487	putative cytochrome c/conserved exported protein of unknown function
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/15);  total (21/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAG  >  AE015451/3950952‑3951374                                                                                                                                                                                                                |                                                                                                                                                                                                                        tGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTGCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCCTATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGc                                                                                                                                                                                                      <  1:181741/227‑1 (MQ=255)           cTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGc                                                                                                                                                                                    >  1:776298/1‑235 (MQ=255)                                           ggtcggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAggc                                                                                                                                                              >  1:159931/1‑225 (MQ=255)                                               ggtcAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGc                                                                                                                                                                                          >  1:393247/1‑193 (MQ=255)                                                   aGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACt                                                                                                                                                 >  1:485549/1‑230 (MQ=255)                                                      tCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGCCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCcaca                                                                                                                                                                                <  2:617632/196‑1 (MQ=255)                                                                  gcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCCTGATACCCACCCACCGGCCAg                                                                                                                                   >  2:355967/1‑229 (MQ=255)                                                                    tGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCCCCCACAAAAGCGGATGCCACAGAGGCCAATTcccc                                                                                                                                                                  >  1:654803/1‑196 (MQ=255)                                                                           gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCGcc                                                                                                                                                                                                   >  2:179837/1‑156 (MQ=255)                                                                                              cGCAGGTAAGCCATAGCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGc                                                                                                                                                                                                         <  2:789253/131‑1 (MQ=255)                                                                                                       gCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGGGAAAGGGCTGCGCCGCAGCCCCGGCAGTGGCAGCTCAAGCGCCGAAAGAGGGGAGCGCTGCGCGCCCCTTTCGCGACGCAAGGc                                                                                                                                                                                                         >  2:604467/1‑122 (MQ=37)                                                                                                          atatCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTggg                                                                                                      >  1:789396/1‑218 (MQ=255)                                                                                                              cGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACaga                                                                                                                                                                              >  1:216358/1‑142 (MQ=255)                                                                                                                cAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCGGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAg                                                                                                                                                                <  2:219960/154‑1 (MQ=255)                                                                                                                     ccggtgGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCCGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGa                                                                         <  1:203554/233‑1 (MQ=255)                                                                                                                     ccTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                                                                                                       >  1:482601/1‑190 (MQ=255)                                                                                                                     ccTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCACCAggc                                                                                                                                                              >  1:48083/1‑151 (MQ=255)                                                                                                                           ggAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCCAGGCCGCGCCCACAAAAGCGGATGCcac                                                                                                                                                                                 >  2:164005/1‑126 (MQ=255)                                                                                                                                  ccTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGCCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGc                                                            >  1:269412/1‑236 (MQ=255)                                                                                                                                    ttGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGACGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCCGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTg                                                                   >  1:460418/1‑227 (MQ=255)                                                                                                                                      gCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTgg                                                                                         >  1:698492/1‑203 (MQ=255)                                                                                                                                                 gggCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCggtc                                                                           >  2:443085/1‑204 (MQ=255)                                                                                                                                                               ccccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCACTTCCCCAGGCGGCCCATGATACTCACCAACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGc                                                                           >  1:295833/1‑192 (MQ=255)                                                                                                                                                               ccccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                                                                                                          >  2:181639/1‑145 (MQ=255)                                                                                                                                                               ccccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCa                                 <  1:213345/234‑1 (MQ=255)                                                                                                                                                                cccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCt                              <  1:181641/236‑1 (MQ=255)                                                                                                                                                                                      cAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGcc                                                                                                                                                                                   <  1:761823/65‑1 (MQ=255)                                                                                                                                                                                      cAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACAt                                             <  1:355969/199‑1 (MQ=255)                                                                                                                                                                                      cAGAAAGAGGGGAGCGCTGCGCGCCCCGTTCGCGCCGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTg                                                                                    <  2:116001/160‑1 (MQ=255)                                                                                                                                                                                      cAGAAAGAGCGGCGCGCTGCGCGCTCCTTTCGCGCCGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCGATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCgg            <  1:226974/232‑1 (MQ=255)                                                                                                                                                                                           agagGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCg                                                                              >  1:285294/1‑161 (MQ=255)                                                                                                                                                                                             agGGGAGCGCGGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAg  <  1:787145/235‑1 (MQ=255)                                                                                                                                                                                                   gcgcTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGc                                                                                                     <  2:776296/130‑1 (MQ=255)                                                                                                                                                                                                    cgcTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTc       <  1:756679/223‑1 (MQ=255)                                                                                                                                                                                                          gcgcTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGa          >  1:736224/1‑214 (MQ=255)                                                                                                                                                                                                          gcgcTCCTGTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAg  <  2:740840/222‑1 (MQ=255)                                                                                                                                                                                                             ccctTTTCGGGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTcccc                                                                                                                                                                  <  1:604469/57‑1 (MQ=25)                                                                                                                                                                                                                |                                                                                                                                                                                                                        TGGAGACGCCCTCGCGGTAAAGGCGTTTGCCTGCTTGTTCGTGGTCGGTCAGGTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCGGGCGAATCTTGCTCCAGCGTTTACCTGCAACATCACCGGCAATCAGCTGCGTCGAATGCTGCTCGGGACTCGGCAG  >  AE015451/3950952‑3951374

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