Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,992,784	C→A	Y544* (TAC→TAA)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,992,784	0	C	A	100.0%	120.3 / NA	35	Y544* (TAC→TAA)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/12);  total (23/12)

CTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAG  >  NC_000913/3992658‑3992926                                                                                                                               |                                                                                                                                               cTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATc                                                                                                                          >  4:23513/1‑149 (MQ=255)  tCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCt                                                                                                                         >  6:130530/1‑149 (MQ=255)         cccGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGAt                                                                                                                   >  1:166296/1‑148 (MQ=255)          ccGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTa                                                                                                                 >  5:103613/1‑149 (MQ=255)             tttGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCg                                                                                                              <  5:130530/149‑1 (MQ=255)                 tataTTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCgtgg                                                                                                                                       >  8:70595/1‑119 (MQ=255)                 tataTTAAAGCTAACCTAATTGACGAATTCCCTAACTTGCAGGAGATGGTCGCAGAACTGTCGCACCATTTCCCGCTGCTCTTACCTGCACCGCCACAGAAGGCTCTCTAAAGCCCGTCTGAGATCCGACATCTGGCGATTATCGTTaa                                                                                                          >  4:11852/1‑149 (MQ=255)                  atatTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGccc                                                                                                                                           >  4:34305/1‑115 (MQ=255)                          ggTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAata                                                                                                 >  2:19652/1‑149 (MQ=255)                               cGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGAcc                                                                                            <  2:166296/149‑1 (MQ=255)                                         gATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCCTTTCCCGCTGCGCTTCCCTGCACCGACACCGAAGGCGCTCTAACGCCCGTGTGAGATCAGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAgcggc                                                                                  >  7:47329/1‑149 (MQ=255)                                                        cAGGAGATGGTCGCCGACGTGTCGCAACATTTCCCTCTGCGCTTACCTGCAACGACACCGAAGTAGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACACGAAAGCGGAGTTCCGCAATCAggt                                                                   >  6:77835/1‑149 (MQ=255)                                                                 gTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTCGaa                                                                                                   <  6:93160/108‑1 (MQ=255)                                                                      cGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGCCCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGAtt                                                     >  4:84817/1‑149 (MQ=255)                                                                       gACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGAttt                                                    <  6:56586/149‑1 (MQ=255)                                                                               gCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCt                                                                                                       <  3:68529/90‑1 (MQ=255)                                                                               gCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCt                                                                                                       >  4:68529/1‑90 (MQ=255)                                                                               gCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGc                                            >  4:102035/1‑149 (MQ=255)                                                                               gCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATAATCGTTAACCTGGAATATGACCCGACAGCGGCGGTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGc                                            >  3:78983/1‑149 (MQ=255)                                                                                cACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCt                                           >  7:96590/1‑149 (MQ=255)                                                                                 aCCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTg                                          >  8:23301/1‑149 (MQ=255)                                                                                   cATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGa                                        >  5:2173/1‑149 (MQ=255)                                                                                    aTTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGAt                                       <  6:103613/149‑1 (MQ=255)                                                                                     tttCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATg                                      >  2:8856/1‑149 (MQ=255)                                                                                          cGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCtt                                  >  1:119375/1‑148 (MQ=255)                                                                                           gCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCa                                <  1:19652/149‑1 (MQ=255)                                                                                                gcTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCttt                           <  2:119375/149‑1 (MQ=255)                                                                                                    aCCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGc                        >  2:163720/1‑148 (MQ=255)                                                                                                     ccTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCga                      <  3:102035/149‑1 (MQ=255)                                                                                                     ccTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGAGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCga                      <  5:77835/149‑1 (MQ=255)                                                                                                            cgacaccgaAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCa                >  1:4515/1‑148 (MQ=255)                                                                                                                  cgaAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGc         >  1:159091/1‑149 (MQ=255)                                                                                                                       gcgcTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTggt    <  1:163720/149‑1 (MQ=255)                                                                                                                           tctAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTggtag  >  6:77821/1‑147 (MQ=255)                                                                                                                           tctAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTggtag  <  5:77821/147‑1 (MQ=255)                                                                                                                               |                                                                                                                                               CTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAG  >  NC_000913/3992658‑3992926

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