Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,168,350	Δ1 bp	coding (483/633 nt)	ycfQ ←	repressor for bhsA(ycfR)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,168,348	0	G	.	100.0%	196.7 / NA	47	coding (485/633 nt)	ycfQ	repressor for bhsA(ycfR)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (23/24);  total (23/24)

GTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGCGGGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGATC  >  NC_000913/1168253‑1168476                                                                                                |                                                                                                                                 gTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCaaa                                                                                      >  1:401841/1‑139 (MQ=255)  tttttCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAg                                                                                     >  1:597479/1‑139 (MQ=255)    tttCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTg                                                                                   >  2:1032399/1‑139 (MQ=255)           cGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCtt                                                                            >  2:697605/1‑139 (MQ=255)                          tGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTg                                                             <  2:255318/139‑1 (MQ=255)                           gcgcTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGa                                                            >  2:378503/1‑139 (MQ=255)                                gATCGACATGCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCa                                                       <  2:1051155/139‑1 (MQ=255)                                 aTCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGtt                                                                                                          <  2:308817/87‑1 (MQ=255)                                 aTCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGtt                                                                                                          >  1:308817/1‑87 (MQ=255)                                 aTCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTg                                                                                   >  1:868524/1‑110 (MQ=255)                                 aTCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTg                                                                                   <  2:868524/110‑1 (MQ=255)                                     aCATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTa                                                  <  2:903931/139‑1 (MQ=255)                                         ccccTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCa                                                                                        >  1:1505094/1‑97 (MQ=255)                                         ccccTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCa                                                                                        <  2:1505094/97‑1 (MQ=255)                                             tGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGat                                          <  1:1032399/139‑1 (MQ=255)                                              gAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                         <  2:401841/139‑1 (MQ=255)                                                 taataCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATAtct                                      <  1:1203949/139‑1 (MQ=255)                                                    taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGc                                                                <  2:776479/110‑1 (MQ=255)                                                    taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGatat                                        <  2:836415/134‑1 (MQ=255)                                                    taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGatat                                        >  1:836415/1‑134 (MQ=255)                                                    taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                         >  1:1468495/1‑133 (MQ=255)                                                    taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                         <  2:1468495/133‑1 (MQ=255)                                                    taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCt                                   >  1:128270/1‑139 (MQ=255)                                                    taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCt                                   >  2:1412568/1‑139 (MQ=255)                                                    taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCt                                   >  2:18811/1‑139 (MQ=255)                                                    taCAATTAAGGAATTCTGCCAGATGTGTCACGTAACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGc                                                                >  1:776479/1‑110 (MQ=255)                                                     aCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTg                                  <  2:128270/139‑1 (MQ=255)                                                        aTTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTc                                                        >  1:740931/1‑114 (MQ=255)                                                        aTTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTc                                                        <  2:740931/114‑1 (MQ=255)                                                         ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATaaa                                                                                                  >  1:862341/1‑71 (MQ=255)                                                         ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATaaa                                                                                                  <  2:862341/71‑1 (MQ=255)                                                         ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATaaa                                                                                                  <  2:861853/71‑1 (MQ=255)                                                         ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATaaa                                                                                                  >  1:861853/1‑71 (MQ=255)                                                         ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                         <  2:1402223/128‑1 (MQ=255)                                                         ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                         >  1:1402223/1‑128 (MQ=255)                                                            aGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTgg                           >  2:1442778/1‑139 (MQ=255)                                                            aGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTgg                           <  1:378503/139‑1 (MQ=255)                                                                  tCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGagag                     >  2:1417128/1‑139 (MQ=255)                                                                    tGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAggg                   <  1:1442778/139‑1 (MQ=255)                                                                     gCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGggt                  <  2:176101/139‑1 (MQ=255)                                                                        aGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTgcg                                                                                 >  1:890539/1‑73 (MQ=255)                                                                        aGATGTGTCACGTCACAGTGAGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAATTgcg                                                                                 <  2:890539/73‑1 (MQ=255)                                                                         gATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGgtggtg              <  2:1314661/139‑1 (MQ=255)                                                                         gATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGgtggtg              <  2:597479/139‑1 (MQ=255)                                                                                     cacaGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGATc  >  2:322651/1‑139 (MQ=255)                                                                                        aGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCa                                                                                        >  2:842001/1‑50 (MQ=255)                                                                                        aGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCa                                                                                        <  1:842001/50‑1 (MQ=255)                                                                                                |                                                                                                                                 GTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGCGGGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGATC  >  NC_000913/1168253‑1168476

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