Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,885,132	C→A	I19I (ATC→ATA)	yidC →	membrane protein insertase YidC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,885,132	0	C	A	100.0%	95.7 / NA	28	I19I (ATC→ATA)	yidC	membrane protein insertase YidC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/12);  total (16/12)

ATTAAAATGCCACCCTTTACACCCTGGTGGTGACGATCCCGTCCCGCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATCTGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTT  >  NC_000913/3884998‑3885264                                                                                                                                       |                                                                                                                                     aTTAAAATGCCACCCTTTACACCCTGGTGGTGACGATCCCGTCCCGCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCa                                                                                                                                 <  2:160707/140‑1 (MQ=255)      aaTGCCACCCTTTACACCCTGGTGGTGACGATCCCGTCCCGCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCt                                                                                                                            <  2:204336/140‑1 (MQ=255)       aTGCCACCCTTTACACCCTGGTGGTGACGATCCCGTCCCGCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTg                                                                                                                           >  2:140057/1‑140 (MQ=255)             cccTTTACACCCTGGTGGTGACGATCCCGTCCCGCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCa                                                                                                                     >  2:155582/1‑140 (MQ=255)                        cTGGTGGTGACGATCCCGTCCCGCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAAcc                                                                                                          >  1:153200/1‑140 (MQ=255)                            tggtgACGATCCCGTCCCGCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCa                                                                                                      >  1:433032/1‑140 (MQ=255)                                gACGATCCCGTCCCGCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCt                                                                                                  <  1:204736/140‑1 (MQ=255)                                 aCGATCCCGTCCCGCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTc                                                                                                 >  1:144686/1‑140 (MQ=255)                                     tcccgtcccgCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGc                                                                                             >  2:439977/1‑140 (MQ=255)                                             cgcccgGACCATTTGATACCAGATAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGa                                                                                     <  1:411274/140‑1 (MQ=255)                                             cgcccgGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGa                                                                                     <  2:306683/140‑1 (MQ=255)                                                              aCCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAGAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGacc                                                                    >  1:467363/1‑140 (MQ=255)                                                                cAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGaccac                                                                  >  2:99619/1‑140 (MQ=255)                                                                cAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGaccac                                                                  >  2:169088/1‑140 (MQ=255)                                                                  gagaACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCg                                                                >  1:43908/1‑140 (MQ=255)                                                                  gagaACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCg                                                                >  1:345901/1‑140 (MQ=255)                                                                   agaACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGc                                                               >  1:388857/1‑140 (MQ=255)                                                                   agaACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCACGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGc                                                               >  1:460317/1‑140 (MQ=255)                                                                                 ggATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGcgccg                                                 <  1:99619/140‑1 (MQ=255)                                                                                  gATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGcgccgc                                                >  1:367859/1‑140 (MQ=255)                                                                                          acgcaaTCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAgg                                        <  2:345901/140‑1 (MQ=255)                                                                                                      aGTCATCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCa                            <  1:159796/140‑1 (MQ=255)                                                                                                           tCGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGc                       <  2:393819/140‑1 (MQ=255)                                                                                                            cGCTTTGCTGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGcc                      <  2:350720/140‑1 (MQ=255)                                                                                                                    tGTTCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGaaa              <  1:169088/140‑1 (MQ=255)                                                                                                                      ttCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACt            <  2:433032/140‑1 (MQ=255)                                                                                                                       tCGTGTCTTTCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTg           >  1:74240/1‑140 (MQ=255)                                                                                                                                tCATGATATGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGtt  >  2:444562/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                     ATTAAAATGCCACCCTTTACACCCTGGTGGTGACGATCCCGTCCCGCCCGGACCATTTGATACCAGAGAACACTAACGATGGATTCGCAACGCAATCTTTTAGTCATCGCTTTGCTGTTCGTGTCTTTCATGATCTGGCAAGCCTGGGAGCAGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTT  >  NC_000913/3884998‑3885264

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