Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,036,781	G→A	F442F (TTC→TTT)	recJ ←	ssDNA‑specific exonuclease RecJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,036,781	0	G	A	100.0%	82.6 / NA	29	F442F (TTC→TTT)	recJ	ssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/11);  total (18/11)

GACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACA  >  NC_000913/3036535‑3037000                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                            acccGGTTCGCCCATCCCCTGCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCCCCGGGAACATCTGCCCCCCCGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGTTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAAtt                                                                                                                                                                                                                         <  2:285663/249‑1 (MQ=255)    cGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATCGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCACACCGTTGT‑TGAAAGAGTTTAAATTTAt                                                                                                                                                                                                                      >  2:24641/1‑251 (MQ=255)           cccaccccctacAAATGACGTTCGCCCCCCGGCCCCGGTTGCGGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGG‑CCAGCATCGCGCGGTAGCTGCGCCACTTCCATGGTCTTTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCGTGCAATATCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAATGAGTTTAAATTTATCCTCTTc                                                                                                                                                                                                               <  2:343816/240‑1 (MQ=255)                    ttCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACa                                                                                                                                                                                                       >  2:247783/1‑250 (MQ=255)                                  gCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCACACAGCGGCTCCGGGCACATCTGCCCAAACGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCACTAGCGAAGGGTCCAGCCACTCAGTAACCAGTTAGCCAACCCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCaacc                                                                                                                                                                                        >  1:321064/1‑248 (MQ=255)                                  gCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTACGGGAACATCTGCCCCCACGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCg                                                                                                                                                                                        >  1:163836/1‑251 (MQ=255)                                               ttttgaggaagaaggAAATGACCGACAAACAGCGGC‑CCCGGTACAATTGCCCCCCCGGGGCCAGCAGCGCGCAGCCGCGGGCCCATTTCCATGGTCATTTCGCCGGGG‑CTAAACGAACCGTCTAATCCCACTTCCCCTTGCAATGGTAAAGGGCCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑GGAAAGAGTTTAAATTTACCCTCTCCCAGCGACAAACCCGCCGCCATCGCATGACCGCCaa                                                                                                                                                                            <  2:241102/237‑1 (MQ=255)                                                gtgccagaggcgagaaaggcccGTCAAAAAGCGGCTCCGGAAACTTCCCCCCCCAGGGC‑CCACCAGCCCGCAGCAGCTGC‑CCCCCTCCCTGGTCATTTCCGCCGGGGTTTAACGGCCCCTCTGATACCCCTTCCCCTGGCAATACAGAGGGCCCCACCCACTCAGTAACCAGTCCGCCAACCCGT‑GTGTAAAAGAGTTTAAATTTACCCCCTCCCAGCGACAAACCCGCCGCCATCGCATGACCGCCaa                                                                                                                                                                            <  2:53543/229‑1 (MQ=255)                                                                  gACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTACCCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCCTGCCCGCCAAACTTCTGCCTCATGCCaagg                                                                                                                                                        >  2:251341/1‑248 (MQ=255)                                                                         aaaCAGCGGCTCCGGGAACATCTGCCCCCACGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCCGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAgtg                                                                                                                                                 >  1:132570/1‑251 (MQ=255)                                                                           aCAGCGGCTCCGGGAACATCTGCCCCCACGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAgtggg                                                                                                                                               >  1:87415/1‑249 (MQ=255)                                                                              ccGGCTCCGGAACCTTCTCCCCCCCCGGG‑CCCGCCTCGCGCAGCAGCTGCCCCCCTTCCATGTCCTTTTCGGCCGGG‑TTTACCGGCCCTTCTTATCCCCCTCCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTACCCAGTTCGCCAACCCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAACCCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCt                                                                                                                                            <  1:352956/250‑1 (MQ=255)                                                                                ggCTCCGGGAACATCTGCCCCCACGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATCCCCCTTCGCCTTGAAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGTGTAGAGTGTGTCTaa                                                                                                                                          <  1:262042/251‑1 (MQ=255)                                                                                       ggAACATCTGCCCCCACGGG‑CCAGCATCGCGAAGCAGCTCCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTcc                                                                                                                                   >  1:45011/1‑251 (MQ=255)                                                                                                                catcgcgcctcccctgagcCACTTCCATGATCATTTCGGCCGGG‑ATTACCGGACCGTCTGATACCACTTCGCATTGCAATACCGAAGCGTCCAACCACTCTGTAAACAGTTCGCCAAAACATTGT‑TGAAAGAGGTTAACTTTATACTCTTCCAGCGACACCGCAGCCGCCATCGCATGACCGCCAACCTTCAGCATCACGCCAGGGTAGGGGGTGTCTAATCGCTCCAGCGCATCACGCATATCCGGcccc                                                                                                           >  2:122652/18‑251 (MQ=255)                                                                                                                  gcgcgCAGCAGCTGTGCCATGTCCATGGTCATTCCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCTTTGCAATAGGGAAGTGTCCAGCCTCTCAGTAACCAGTTCGCCAAACGGTTGT‑TAAAAGAGTTTAAATTTATCCTCTTCCAGCAACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTg                                                                                                         <  1:110522/250‑1 (MQ=255)                                                                                                                                cgcCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCCCTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGACTGGAGCGACCgg                                                                                           >  2:218630/1‑251 (MQ=255)                                                                                                                                    aCTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAAcc                                                                                       >  2:34605/1‑251 (MQ=255)                                                                                                                                              tCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGCGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGt                                                                              >  2:258645/1‑250 (MQ=255)                                                                                                                                                       ccGGG‑CTTAACGGATAGTCTGATAACACTTCGCCTTGCACAAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCCGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGGGTGTCTAATCGCTCCAGTGCATTACGCATATGCAGCCCATGAAGGGAGCGAACGGCACCTTTTAGCGTAATGTCAccc                                                                    >  2:317485/1‑250 (MQ=255)                                                                                                                                                                    gACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCAGCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAAc                                                        <  1:218630/251‑1 (MQ=255)                                                                                                                                                                      ccGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAAc                                                                                        >  1:24922/1‑217 (MQ=255)                                                                                                                                                                      ccGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAAcgc                                                      >  2:122384/1‑251 (MQ=255)                                                                                                                                                                         tCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGAAACCAGTTCGCCAAACCGTTGT‑TGCAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGAACCGTCACCTGCTGGCGCAAACGCGAt                                                   >  1:39378/1‑251 (MQ=255)                                                                                                                                                                                          ttcaaaTAGCAAAGGGTCCAGCTCCTCAGTAACCGGTTCGCCAACCCGGTGT‑GAAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAGCCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTGCAGTGCATCACGCATATGCAGCCCCTGAATGGATCGACCGGAACCTTGCATCGTTCCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGaaa                                    <  1:159521/245‑1 (MQ=255)                                                                                                                                                                                               agtaggaaggtctcAGCCACACTGTACCCAGTTCCCAAAACCGGTGT‑TGACAGAGTTTAAATTTATCCTCTTCATGCGACATCCCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTACTCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGctcc                             <  2:61061/238‑2 (MQ=255)                                                                                                                                                                                                                 aCTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGca              >  1:185323/1‑247 (MQ=255)                                                                                                                                                                                                                    taGTATCCACTTCACCAAACCGTTGT‑TGAAAGATTTTAAATTTATCCTCTTCCAGCGCCAAACCCGCCGCCATCGCATGACCGCCAAACTGCAGCATCATGCCAGGGGAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAAt         <  2:139229/249‑1 (MQ=255)                                                                                                                                                                                                                           cAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTAAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGCTGCGCGCAGCCAGAATACCGAcc  >  2:3942/1‑249 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                            GACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGG‑CCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGG‑CTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGT‑TGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGCCATCGCATGACCGCCAAACTTCAGCATCATGCCAGGGTAGAGTGTGTCTAATCGCTCCAGCGCATCACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACA  >  NC_000913/3036535‑3037000

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