Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,993,145	C→A	100%	R665S (CGT→AGT)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,145	0	C	A	100.0%	146.3 / NA	43	R665S (CGT→AGT)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/23);  total (20/23)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGC  >  NC_000913/3993015‑3993286                                                                                                                                   |                                                                                                                                              aaCTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGACTGGAGTCTTTCCAGCACCCGcc                                                                                                                             >  2:238601/1‑149 (MQ=255)         tcGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAAcc                                                                                                                     >  1:150587/1‑149 (MQ=255)           ggCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCgg                                                                                                                   >  1:125098/1‑149 (MQ=255)            gCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCgg                                                                                                                   >  3:45019/1‑148 (MQ=255)                aaTGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGt                                                                                                               <  4:160925/148‑1 (MQ=255)                 aTGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGttt                                                                                                             <  2:23131/149‑1 (MQ=255)                     aTCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAg                                                                                                         >  1:4923/1‑149 (MQ=255)                      tCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAgg                                                                                                        <  1:119168/149‑1 (MQ=255)                             gccgcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTgcg                                                                                                 <  2:125098/149‑1 (MQ=255)                              ccgcACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTgcg                                                                                                 >  4:114204/1‑148 (MQ=255)                                       cAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTgg                                                                                        <  3:114204/148‑1 (MQ=255)                                            aGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCaaa                                                                                   <  3:278689/148‑1 (MQ=255)                                            aGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCaaa                                                                                   <  4:309600/148‑1 (MQ=255)                                               gTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCt                                                                                >  3:124550/1‑148 (MQ=255)                                                     gTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgt                                                                          >  3:326923/1‑148 (MQ=255)                                                             ttATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTAGTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGGTTCAAGGCGCTGCGCGTGTCTGGTCAAACCTGGGGGTTGttcttc                                                                  <  4:248629/148‑1 (MQ=255)                                                                 aGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAAc                                                              >  4:231781/1‑148 (MQ=255)                                                                   ccAGCATCTGCGCGGCTTAATTCGTACCCGCGTGCAGCAACTGGTTTCTGAGTATATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGGTGCGCGTTTCTGGTCAAACCTGGGGGTTGTCCTTCGAACGc                                                            <  3:229245/148‑1 (MQ=255)                                                                    cAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCt                                                          >  2:170216/1‑149 (MQ=255)                                                                       cATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGa                                                        <  3:123713/148‑1 (MQ=255)                                                                        aTCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGaa                                                       <  4:45019/148‑1 (MQ=255)                                                                            gcgcGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTa                                                   <  3:128507/148‑1 (MQ=255)                                                                            gcgcGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTAt                                                  >  2:305457/1‑149 (MQ=255)                                                                            gcgcGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTAt                                                  <  1:305457/149‑1 (MQ=255)                                                                                   ttAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgttgtt                                                                      >  2:164943/1‑122 (MQ=255)                                                                                   ttAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGgttgtt                                                                      <  1:164943/122‑1 (MQ=255)                                                                                            aCTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTgg                                   >  4:289672/1‑148 (MQ=255)                                                                                                gcgTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCACTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGaaaa                               <  4:330958/148‑1 (MQ=255)                                                                                                  gTGCAGCAACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGc                            <  1:342922/149‑1 (MQ=255)                                                                                                        caACTGGTTTCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGa                      <  2:201910/149‑1 (MQ=255)                                                                                                               tttCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTa                <  3:79852/148‑1 (MQ=255)                                                                                                               tttCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTAt               <  1:248827/149‑1 (MQ=255)                                                                                                                ttCTGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAACCCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCTTGAATGTATCGGTAAAGAACCTGGAAAACGCCATCGAGTTTTAt               <  3:242725/148‑1 (MQ=255)                                                                                                                   tGAGTGTATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGgcg           >  1:339748/1‑149 (MQ=255)                                                                                                                       tgtATTGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGa        >  4:285783/1‑148 (MQ=255)                                                                                                                           ttGAATTTAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTg                                    >  2:217970/1‑116 (MQ=255)                                                                                                                           ttGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTg                                    <  1:217970/116‑1 (MQ=255)                                                                                                                           ttGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCg   >  2:373199/1‑149 (MQ=255)                                                                                                                            tGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTa                >  3:390604/1‑135 (MQ=255)                                                                                                                            tGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTa                <  4:390604/135‑1 (MQ=255)                                                                                                                            tGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCg   >  3:214835/1‑148 (MQ=255)                                                                                                                            tGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCg   >  3:396617/1‑148 (MQ=255)                                                                                                                            tGAATTGAGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGc  <  2:150587/149‑1 (MQ=255)                                                                                                                                   |                                                                                                                                              AACTATTCTCGGCAAAATGCATCAGGACGCCGCACCGCCAGATAGCGTGGAAGTCTTCTGTTATAGCCAGCATCTGCGCGGCTTAATTCGTACTCGCGTGCAGCAACTGGTTTCTGAGTGTATTGAATTGCGTCTTTCCAGCACCCGCCAGGAAACCGGGCGTTTCAAGGCGCTGCGCGTTTCTGGTCAAACCTGGGGGTTGTTCTTCGAACGCCTGAATGTATCGGTACAGAAACTGGAAAACGCCATCGAGTTTTATGGCGCGATTTCGC  >  NC_000913/3993015‑3993286

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