Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,172,139	A→G	I609T (ATT→ACT)	mfd ←	transcription‑repair coupling factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,172,139	0	A	G	100.0%	149.5 / NA	41	I609T (ATT→ACT)	mfd	transcription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/24);  total (17/24)

AGGCACCAGCACCGCCACCTGCTTGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAATGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGTTGCGCGTAGATATCCAGCAATTCCGCCGCCACATCACGCACTT  >  NC_000913/1171994‑1172286                                                                                                                                                  |                                                                                                                                                    aGGCACCAGCACCGCCACCTGCTTGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGcc                                                                                                                                                >  2:261568/1‑151 (MQ=255) aGGCACCAGCACCGCCACCTGCTTGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGcc                                                                                                                                                >  1:71826/1‑151 (MQ=255)    cACCAGCACCGCCACCTGCTTGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTgc                                                                                                                                             <  1:124761/151‑1 (MQ=255)     aCCAGCACCGCCACCTGCTTGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTgcg                                                                                                                                            <  2:412291/151‑1 (MQ=255)        aGCACCGCCACCTGCTTGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCt                                                                                                                                         >  1:309886/1‑151 (MQ=255)         gCACCGCCACCTGCTTGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTg                                                                                                                                        <  2:243087/151‑1 (MQ=255)           aCCGCCACCTGCTTGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTg                                                                                                                                        <  2:79617/149‑1 (MQ=255)            ccGCCACCTGCTTGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGAt                                                                                                                                      <  2:372444/150‑1 (MQ=255)                       ttGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGttt                                                                                                                          >  1:677023/1‑151 (MQ=255)                                   aCTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGaa                                                                                                                <  1:466506/149‑1 (MQ=255)                                   aCTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGaa                                                                                                                <  1:533056/149‑1 (MQ=255)                                          ggAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCa                                                                                                       <  1:634519/151‑1 (MQ=255)                                          ggAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCa                                                                                                       <  2:270453/151‑1 (MQ=255)                                                      gcATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATa                                                                                           <  2:287785/151‑1 (MQ=255)                                                       cATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATAc                                                                                          >  2:255367/1‑151 (MQ=255)                                                            ccACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTc                                                                                     >  1:574764/1‑151 (MQ=255)                                                                ttCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACg                                                                                  >  1:658789/1‑150 (MQ=255)                                                                    gTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGAt                                                                                <  2:556000/148‑1 (MQ=255)                                                                    gTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGAt                                                                                <  1:294330/148‑1 (MQ=255)                                                                    gTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCg                                                                              <  2:295481/150‑1 (MQ=255)                                                                        tACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGttt                                                                         >  1:203165/1‑151 (MQ=255)                                                                          ccAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTa                                                                        >  2:167063/1‑150 (MQ=255)                                                                              aGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAAcgc                                                                   <  1:478226/151‑1 (MQ=255)                                                                              aGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAAcgc                                                                   <  1:422393/151‑1 (MQ=255)                                                                                  aaCATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGa                                                                 >  1:666220/1‑149 (MQ=255)                                                                                    cATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGc                                                              <  2:495212/150‑1 (MQ=255)                                                                                       cgccgcACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCctc                                                          >  2:22517/1‑151 (MQ=255)                                                                                       cgccgcACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCct                                                           <  2:266299/150‑1 (MQ=255)                                                                                       cgccgaACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCct                                                           >  2:340089/1‑150 (MQ=255)                                                                                          cgcACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCttt                                                       <  2:117538/151‑1 (MQ=255)                                                                                           gcACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTg                                                      >  2:604517/1‑151 (MQ=255)                                                                                                    aCGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGt                                             <  1:595586/151‑1 (MQ=255)                                                                                                        tCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGTTgc                                          <  2:142608/150‑1 (MQ=255)                                                                                                              gCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGTTGCGCGt                                      >  1:115865/1‑148 (MQ=255)                                                                                                               ccAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGTTGCGCGTAGat                                  <  2:52011/151‑1 (MQ=255)                                                                                                                cAGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGTTGCGCGTAGat                                  >  1:658489/1‑150 (MQ=255)                                                                                                                 aGCGGCTGACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATACGGCGTGGTTTCACACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCGTCTTTGGCGGCGCGTTGCGCGTAGa                                   <  1:62114/148‑1 (MQ=255)                                                                                                                        gACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGTTGCGCGTAGATATCCAGCaa                         <  2:75625/151‑1 (MQ=255)                                                                                                                        gACACATGTCGCTAAGTACCGCATTAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGTTGCGCGTAGATATCCAGCa                          <  2:88147/150‑1 (MQ=255)                                                                                                                                               ttAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGTTGCGCGTAGATATCCAGCAATTCCGCCGCCACATCACGCACtt  >  2:585850/1‑151 (MQ=255)                                                                                                                                               ttAGTGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGTTGCGCGTAGATATCCAGCAATTCCGCCGCCACATCACGCACtt  >  2:391867/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AGGCACCAGCACCGCCACCTGCTTGTGGTTATCTACTGCCAGGAAAGCTGCGCGCATCGCCACTTCTGTTTTACCAAAGCCAACATCGCCGCACACCAGACGATCCATTGCCAGCGGCTGACACATGTCGCTAAGTACCGCATTAATGGCCTGCGCCTGATCCGGCGTGGTTTCAAACGGGAAGCTGTCGCAGAACAACTGATACTGCTCACGATCGTGTTTAAACGCGAAGCCCTCTTTGGCGGCGCGTTGCGCGTAGATATCCAGCAATTCCGCCGCCACATCACGCACTT  >  NC_000913/1171994‑1172286

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