Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,407,268	A→G	intergenic (+2/‑107)	accC → / → yhdT	biotin carboxylase/DUF997 domain‑containing protein YhdT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,407,268	0	A	G	100.0%	126.5 / NA	36	intergenic (+2/‑107)	accC/yhdT	biotin carboxylase/DUF997 domain‑containing protein YhdT
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/21);  total (15/21)

TGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGT  >  NC_000913/3407133‑3407386                                                                                                                                        |                                                                                                                       tGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGc                                                                                                                    >  2:1096123/1‑140 (MQ=255)  gCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGgctgct                                                                                                                   <  1:439933/140‑1 (MQ=255)  gCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACATCCAGCATGGTGGCACTATCATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGgctgct                                                                                                                   <  1:563514/140‑1 (MQ=255)             tcatcGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCaaa                                                                                                        >  2:1005249/1‑140 (MQ=255)                 cGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGc                                                                                                    >  1:1154217/1‑140 (MQ=255)                           aaaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTcc                                                                                          <  1:741081/140‑1 (MQ=255)                               ccAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGcc                                                                                      <  2:279927/140‑1 (MQ=255)                                 aaCGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCtt                                                                                    >  2:984828/1‑140 (MQ=255)                                 aaCGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCtt                                                                                    >  2:66818/1‑140 (MQ=255)                                   cGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCtttt                                                                                  >  2:198156/1‑140 (MQ=255)                                      tGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                               <  2:318526/140‑1 (MQ=255)                                      tGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                               <  2:331691/140‑1 (MQ=255)                                         tCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGCGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTttatta                                                                            >  2:835240/1‑140 (MQ=255)                                                aTCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGa                                                                     >  2:865199/1‑140 (MQ=255)                                                  ccGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTttatta                                                                            <  2:209482/131‑1 (MQ=255)                                                  ccGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTttatta                                                                            >  1:209482/1‑131 (MQ=255)                                                  ccGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATc                                                                   >  1:978531/1‑140 (MQ=255)                                                     catcatGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGAc                                                                >  2:602729/1‑140 (MQ=255)                                                     catcatGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGAc                                                                <  1:984828/140‑1 (MQ=255)                                                      atcatGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACa                                                               <  1:1096123/140‑1 (MQ=255)                                                          tgaatgaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAccc                                                           <  2:1087550/140‑1 (MQ=255)                                                          tgaatgaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAccc                                                           <  1:735346/140‑1 (MQ=255)                                                          tgaatgaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAccc                                                           <  1:340574/140‑1 (MQ=255)                                                             atgaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCa                                                        <  2:1154217/140‑1 (MQ=255)                                                                  gagaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATa                                                      <  2:1099439/137‑1 (MQ=255)                                                                  gagaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATa                                                      >  1:1099439/1‑137 (MQ=255)                                                                             gCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCg                                        <  2:978531/140‑1 (MQ=255)                                                                             gCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCg                                        <  1:865199/140‑1 (MQ=255)                                                                              cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGc                                       <  1:602729/140‑1 (MQ=255)                                                                                  gtggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATcccc                                         <  1:628567/134‑1 (MQ=255)                                                                                  gtggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATcccc                                         >  2:628567/1‑134 (MQ=255)                                                                                        cTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCAc                              >  2:257900/1‑139 (MQ=255)                                                                                        cTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCAc                              <  1:257900/139‑1 (MQ=255)                                                                                           aCATCTACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCAt                          <  2:1029504/140‑1 (MQ=255)                                                                                                   tATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACa                  >  1:122675/1‑140 (MQ=255)                                                                                                          agaAAAAACTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg           <  2:182631/140‑1 (MQ=255)                                                                                                                 aCTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAA‑‑CCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGt  <  1:138672/140‑1 (MQ=255)                                                                                                                 aCTCGGTCTTCAGGAAAAATAAGGCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAA‑‑CCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGt  <  1:1080221/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                       TGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGT  >  NC_000913/3407133‑3407386

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: ‑