Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,985,816	T→C	intergenic (‑224/‑31)	yqeF ← / → yqeG	putative acyltransferase/putative transporter YqeG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,985,816	0	T	C	100.0%	125.3 / NA	34	intergenic (‑224/‑31)	yqeF/yqeG	putative acyltransferase/putative transporter YqeG
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/13);  total (21/13)

GCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAG  >  NC_000913/2985602‑2986060                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     gCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTTATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccc                                                                                                                                                                                                                  >  1:253584‑M1/1‑240 (MQ=255)     cAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGCAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgac                                                                                                                                                                                                              >  1:57527‑M1/1‑236 (MQ=255)                          aTGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcat                                                                                                                                                                                         >  1:380469‑M1/1‑215 (MQ=255)                            gATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataa                                                                                                                                                                                       <  1:446004‑M1/251‑39 (MQ=255)                            gATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcata                                                                                                                                                                                        <  1:317165‑M1/250‑38 (MQ=255)                                          gTTATCTTAATTGATTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGTGCAGCAGTCAGAGACATAACTGGCACGTACGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTCATCTCCCATACTTAACCTGGTTTATGGGAAATTGCCCTCCATTCTGTTTAATTTGATGATGATACGTTCaaaaaaatccccgacatttaccggggttgtgagcataacgtaatgcttattt                                                                                                                                                                         >  2:60459‑M1/1‑199 (MQ=255)                                                   aTTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacg                                                                                                                                                                >  1:304898‑M1/1‑190 (MQ=255)                                                   aTTGTTTAAAAAAAGTGATTTTTATCAAAAAGGAAATATGACTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAAAGGGGTCGATAGGGCGTTAATCTACCATACATAACCTGGTTTATGGTAAAATGCCCTACATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgggcataacgtcatgcttattttaccggacc                                                                                                                                                                >  2:401323‑M1/1‑190 (MQ=255)                                                    ttGTTTCAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgc                                                                                                                                                               >  2:141416‑M1/1‑189 (MQ=255)                                                      gTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGGAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcat                                                                                                                                                             <  1:13932‑M1/251‑65 (MQ=255)                                                             aaaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGCTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccggg                                                                                                                                                      >  1:418828‑M1/1‑180 (MQ=255)                                                              aaaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccggga                                                                                                                                                     <  1:326742‑M1/251‑73 (MQ=255)                                                               aaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgc                                                                                                                                                                                                        <  1:215836‑M1/199‑22 (MQ=255)                                                               aaGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgc                                                                                                                                                                                                        >  2:215836‑M1/1‑178 (MQ=255)                                                                    aTTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGCTAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAACTTGCCCTCCATTCTGTTTAATTTGTAGATGCTACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacaga                                                                                                                                                >  2:137101‑M1/1‑173 (MQ=255)                                                                     tttttATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGATTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATATCCCATACTTCACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaataaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccggggacggaat                                                                                                                                              >  2:38041‑M1/1‑172 (MQ=255)                                                                            cacaAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcg                                                                                                                                       >  1:371145‑M1/1‑165 (MQ=255)                                                                                aaGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaa                                                                                                                                                                                  <  1:136235‑M1/204‑44 (MQ=255)                                                                                aaGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaa                                                                                                                                                                                  >  2:136235‑M1/1‑161 (MQ=255)                                                                                   gAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatgg                                                                                                                                >  1:247690‑M1/1‑158 (MQ=255)                                                                                                         agagACAGAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcat                                                                                                          <  2:188257‑M1/251‑116 (MQ=255)                                                                                                             aCATAACTGGCACGTAGGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgtttgtgaacagcattacc                                                                                                      >  2:49851‑M1/1‑132 (MQ=255)                                                                                                                                     aaCCACTAACCCTCCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccaga                                                                              <  2:293954‑M1/251‑144 (MQ=255)                                                                                                                                         aCTAACCCACCAATCGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccataaggtcgatacccattcccagacctg                                                                          >  2:307136‑M1/1‑104 (MQ=255)                                                                                                                                                          ggggTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacgggcgataccaattcccagacctgctgtcggaggtaagcc                                                          >  1:289493‑M1/1‑87 (MQ=255)                                                                                                                                                               aGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgtcggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagcataacgtcgcggatggtatggctgttggttaacagccttaacataacgtccatacgcattcccagaccgtctgtctgcggtcgagcatattc                                                    >  2:109151‑M1/1‑82 (MQ=255)                                                                                                                                                               aGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgtt                                                     <  2:81377‑M1/250‑169 (MQ=255)                                                                                                                                                                      gCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcg                                             <  1:215070‑M1/251‑177 (MQ=255)                                                                                                                                                                           aaTCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCgaaaaaacccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccagacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgac                                        <  2:148832‑M1/251‑182 (MQ=255)                                                                                                                                                                                cccATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgta                                     <  1:322323‑M1/249‑185 (MQ=255)                                                                                                                                                                                       ttAACCTGGTTTCTGGTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcgggaagccatgttccagtgcggtgacgtcatctttatc                            >  2:435588‑M1/1‑58 (MQ=255)                                                                                                                                                                                                      gTAAATTGCCCTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgcttc             <  1:335624‑M1/251‑209 (MQ=255)                                                                                                                                                                                                                cTCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggggacgtaatcttcatcgtagaacatcgcttcgtcgtcaccc   >  1:317627‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                 tCCATTCTGTTTAATTTGTAGATGATACGTTCtaaaaaaaccccgacatttgccggggttgtgagcataacgtaatgcttattttaccggacgcatcgccgggaacagaataacgtcgcggatggtatggctgttggtgaacagcattaccatacggtcgataccaattcccagacctgctgtcggcggtaagccatgttccagtgcggtgacgtaatcttcatcgtagaacatcgcttcgtcgtcacctg  >  1:11046‑M1/1‑32 (MQ=255)                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     GCACCAATTAACAGTGCAATAGTAAATGATTGTTAATTACAGTTATCTTAATTGTTTAAAAAAAGTGATTTTTATCACAAAGGAAATATGCCTGAGCAGCAGTCAGAGACATAACTGGCACGTAAGGTTTGCAACCACTAACCCACCAATAGAGGGGTAGATAGGGCGTTAATCTCCCATACTTAACCTGGTTTATGGTAAATTGCCCTCCATTTTGTTTAATTTGTAGATGATACGTTCAGATAATGTCTAATATTTGGTCTAAAGAAGAAACTCTGTGGAGTTTCGCGCTCTACGGCACAGCCGTTGGTGCAGGCACGCTCTTCCTTCCTATTCAGTTAGGTTCGGCAGGGGCTGTGGTCCTGTTTATTACTGCTCTGGTCGCCTGGCCTTTAACATATTGGCCACATAAAGCCTTGTGCCAGTTCATCCTCTCATCAAAAACATCAGCAGGTGAAG  >  NC_000913/2985602‑2986060

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