Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007779	4,154,211	T→G	K29T (AAG→ACG)	crp ←	DNA‑binding transcriptional dual regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	4,154,211	0	T	G	100.0%	255.9 / NA	62	K29T (AAG→ACG)	crp	DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (36/26);  total (36/26)

GGATCATTTCTTTACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCTTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATGCGCGGTTATCC  >  NC_007779/4154113‑4154308                                                                                                   |                                                                                                  ggATCATTTCTTTACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTg                                                                                                 >  3:180630/1‑101 (MQ=255)   aTCATTTCTTTACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGct                                                                                               <  3:725331/101‑1 (MQ=255)      aTTTCTTTACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCtt                                                                                            >  2:67224/1‑101 (MQ=255)        ttCTTTACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTgg                                                                                          <  2:107853/101‑1 (MQ=255)         tCTTTACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGa                                                                                         >  1:1384883/1‑101 (MQ=255)           tttACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATg                                                                                       >  1:148750/1‑101 (MQ=255)            ttACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATgg                                                                                      <  3:68940/101‑1 (MQ=255)             tACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATggg                                                                                     <  3:1166915/101‑1 (MQ=255)             tACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATggg                                                                                     >  3:60649/1‑101 (MQ=255)              aCCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGt                                                                                    >  1:589446/1‑101 (MQ=255)                  tctTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACtt                                                                                >  3:607519/1‑101 (MQ=255)                      cGTCATTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGa                                                                            >  1:1524899/1‑101 (MQ=255)                         cTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAAtg                                                                         >  1:269820/1‑101 (MQ=255)                            tGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTgg                                                                      >  3:1482118/1‑101 (MQ=255)                               tCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCaa                                                                   >  2:84954/1‑101 (MQ=255)                                 aGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATg                                                                 <  1:1276082/101‑1 (MQ=255)                                  gCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATga                                                                <  1:1419838/101‑1 (MQ=255)                                    aCTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATgaga                                                              >  3:1211761/1‑101 (MQ=255)                                         cacaGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGa                                                         >  3:721689/1‑101 (MQ=255)                                          acaGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGaa                                                        >  3:528083/1‑101 (MQ=255)                                            agagCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAAcc                                                      <  3:10339/101‑1 (MQ=255)                                              agCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCAt                                                    <  3:152960/101‑1 (MQ=255)                                                ccTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTc                                                  <  1:283035/101‑1 (MQ=255)                                                ccTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTc                                                  <  1:17807/101‑1 (MQ=255)                                                    tAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCgaga                                              <  4:57704/101‑1 (MQ=255)                                                     aaCGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGagag                                             >  1:1494944/1‑101 (MQ=255)                                                     aaCGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGagag                                             <  1:990326/101‑1 (MQ=255)                                                       cGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTc                                           >  1:503965/1‑101 (MQ=255)                                                        gATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCg                                          >  1:430395/1‑101 (MQ=255)                                                              gtaCAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTg                                    >  1:531395/1‑101 (MQ=255)                                                               taCAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGt                                   >  3:1020563/1‑101 (MQ=255)                                                               taCAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGt                                   >  3:1458013/1‑101 (MQ=255)                                                                aCAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGtt                                  >  3:621126/1‑101 (MQ=255)                                                                aCAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGtt                                  >  3:868562/1‑101 (MQ=255)                                                                 cAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGttt                                 <  1:112638/101‑1 (MQ=255)                                                                      tttCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGt                            <  3:101104/101‑1 (MQ=255)                                                                      tttCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGt                            <  3:150068/101‑1 (MQ=255)                                                                      tttCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGt                            <  3:224047/101‑1 (MQ=255)                                                                       ttCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGtt                           <  1:188022/101‑1 (MQ=255)                                                                        tCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGttt                          <  1:1433883/101‑1 (MQ=255)                                                                        tCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGttt                          >  3:693469/1‑101 (MQ=255)                                                                        tCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGttt                          >  3:198836/1‑101 (MQ=255)                                                                        tCCGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGGCTGTTTGCGGttt                          <  4:62684/101‑1 (MQ=255)                                                                          cGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGc                        >  1:1078925/1‑101 (MQ=255)                                                                          cGCTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGc                        >  3:91014/1‑101 (MQ=255)                                                                            cTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCa                      >  2:25703/1‑101 (MQ=255)                                                                            cTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCa                      <  1:688544/101‑1 (MQ=255)                                                                            cTTTTTCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCa                      >  3:1560515/1‑101 (MQ=255)                                                                             tttttCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCaa                     >  3:789575/1‑101 (MQ=255)                                                                              ttttCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAg                    <  1:1154267/101‑1 (MQ=255)                                                                              ttttCACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAg                    <  1:1467110/101‑1 (MQ=255)                                                                                  cACCCTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCAcc                <  1:1200155/101‑1 (MQ=255)                                                                                      cTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATgc            >  3:1238344/1‑101 (MQ=255)                                                                                      cTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATgc            >  1:379130/1‑101 (MQ=255)                                                                                      cTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATgc            >  1:698315/1‑101 (MQ=255)                                                                                      cTGGTGAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATgc            >  1:562252/1‑101 (MQ=255)                                                                                           gAATAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCTGTTGGCCAAACACCAGGCGCGGt       <  4:160988/101‑1 (MQ=255)                                                                                             aTAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATGCGCGGTTa     >  1:1562122/1‑101 (MQ=255)                                                                                              tAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATGCGCGGTTAt    >  3:718437/1‑101 (MQ=255)                                                                                              tAAGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATGCGCGGTTAt    >  1:471718/1‑101 (MQ=255)                                                                                                aGCGTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATGCGCGGGTATcc  <  4:48172/101‑1 (MQ=255)                                                                                                aGCGTGCTCTTGGATGGGTACCTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATgctctttgatcc  <  4:129580/101‑11 (MQ=255)                                                                                                   |                                                                                                  GGATCATTTCTTTACCCTCTTCGTCTTTGATCAGCACTGCCACAGAGCCTTTAACGATGTAGTACAGCGTTTCCGCTTTTTCACCCTGGTGAATAAGCTTGCTCTTGGATGGGTACTTATGAATGTGGCAATGAGACAAGAACCATTCGAGAGTCGGGTCTGTTTGCGGTTTGCCAAGCACCATGCGCGGTTATCC  >  NC_007779/4154113‑4154308

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