Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,796,916	T→C	Y38C (TAT→TGT)	yqaE ←	cyaR sRNA‑regulated protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,796,916	0	T	C	100.0%	185.8 / NA	51	Y38C (TAT→TGT)	yqaE	cyaR sRNA‑regulated protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/30);  total (21/30)

TAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATATACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGAAACCCATATGTAC  >  NC_000913/2796786‑2797034                                                                                                                                   |                                                                                                                       tAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGc                                                                                                                <  2:632718/139‑1 (MQ=255)    cTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAgt                                                                                                             <  2:1028989/139‑1 (MQ=255)      aCCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAgtgt                                                                                                           <  2:223610/139‑1 (MQ=255)      aCCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAgtgt                                                                                                           >  1:1156410/1‑139 (MQ=255)       ccTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTc                                                                                                          >  1:1281198/1‑139 (MQ=255)        cTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCa                                                                                                         >  1:287571/1‑139 (MQ=255)          cATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAAc                                                                                                       >  1:1341194/1‑139 (MQ=255)           aTTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAg                                                                                                                 >  2:973539/1‑128 (MQ=255)           aTTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAg                                                                                                                 <  1:973539/128‑1 (MQ=255)            ttCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCaa                                                                                                        >  1:862666/1‑136 (MQ=255)            ttCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCaa                                                                                                        <  2:862666/136‑1 (MQ=255)            ttCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAg                                                                                                     >  1:1096780/1‑139 (MQ=255)              cGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGaa                                                                                                   <  1:1333958/139‑1 (MQ=255)                             tCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACg                                                                                    <  2:1341194/139‑1 (MQ=255)                                     tgtAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCg                                                                            <  2:287571/139‑1 (MQ=255)                                      gtAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGa                                                                           <  1:792745/139‑1 (MQ=255)                                         gCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAAcc                                                                        <  1:806184/139‑1 (MQ=255)                                               gAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTAc                                                                  >  2:1259901/1‑139 (MQ=255)                                                  cgcgGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGa                                                               >  1:617744/1‑139 (MQ=255)                                                       cgcgTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCa                                                                                                               <  1:1249553/86‑1 (MQ=255)                                                       cgcgTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCa                                                                                                               >  2:1249553/1‑86 (MQ=255)                                                         cgTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCAc                                                        >  1:33420/1‑139 (MQ=255)                                                         cgTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGCTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCAc                                                        >  1:267954/1‑139 (MQ=255)                                                          gTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGcc                                                                                  >  1:938346/1‑112 (MQ=255)                                                          gTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGcc                                                                                  <  2:938346/112‑1 (MQ=255)                                                           ttAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGc                                                      <  2:1281198/139‑1 (MQ=255)                                                                    cgGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAgcggcg                                             <  2:617744/139‑1 (MQ=255)                                                                            aGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAatg                                     <  2:1156410/139‑1 (MQ=255)                                                                               acaGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAatgatg                                  <  2:723508/139‑1 (MQ=255)                                                                                 aGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCg                                                                            >  2:391991/1‑95 (MQ=255)                                                                                 aGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCg                                                                            <  1:391991/95‑1 (MQ=255)                                                                                   gTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCaa                                                                                                        >  1:575870/1‑65 (MQ=255)                                                                                   gTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCaa                                                                                                        <  2:575870/65‑1 (MQ=255)                                                                                                  tGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCTGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGa               <  2:267954/139‑1 (MQ=255)                                                                                                  tGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGa               <  1:1259901/139‑1 (MQ=255)                                                                                                  tGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGa               <  1:712780/139‑1 (MQ=255)                                                                                                  tGCACCCAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGa               <  2:33420/139‑1 (MQ=255)                                                                                                       ccAGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGAAAccc          <  2:220166/139‑1 (MQ=255)                                                                                                         aGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGtgat                             <  2:1305811/118‑1 (MQ=255)                                                                                                         aGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGtgat                             >  1:1305811/1‑118 (MQ=255)                                                                                                         aGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGa               <  1:1035897/132‑1 (MQ=255)                                                                                                         aGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGa               >  2:1035897/1‑132 (MQ=255)                                                                                                         aGAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGAAACCCat        <  2:1096780/139‑1 (MQ=255)                                                                                                          gAACGCGTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGTCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGAAACCCata       >  2:858789/1‑139 (MQ=255)                                                                                                               cgTGAATCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGTCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGAAACCCATATGTAc  <  1:858789/139‑1 (MQ=255)                                                                                                                     tCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGa               <  2:316102/120‑1 (MQ=255)                                                                                                                     tCAAACCGGGGATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGa               >  1:316102/1‑120 (MQ=255)                                                                                                                              ggATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAgcgg                                               >  2:1068378/1‑79 (MQ=255)                                                                                                                              ggATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAgcgg                                               <  1:1068378/79‑1 (MQ=255)                                                                                                                              ggATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAatg                                     <  2:1306663/89‑1 (MQ=255)                                                                                                                              ggATACACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAatg                                     >  1:1306663/1‑89 (MQ=255)                                                                                                                                   |                                                                                                                       TAGCTACCTCATTCGTTAACGGAAACACTCAGAAAGTGTAGCAGAGGAACGCGGCGCGTTAACGCCGCGGGGATGAGCACAGGTTTAATCGCGGGTTTGCACCCAGAACGCGTGAATCAAACCGGGGATATACCCAAGCAGTGTCAACAGAATATTAATAATGAACGCCCAACCGAACCCTTTACCGAGCAGCACGCCGAGCGGCGGCAGAATGATGGTGATGACGATTCTCCAGAAACCCATATGTAC  >  NC_000913/2796786‑2797034

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