Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,992,054	C→T	P301L (CCA→CTA)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,992,054	0	C	T	100.0%	131.5 / NA	41	P301L (CCA→CTA)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/21);  total (20/21)

TCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTA  >  NC_000913/3991921‑3992184                                                                                                                                      |                                                                                                                                   tCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACt                                                                                                                               >  1:122515/1‑139 (MQ=255)   cccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGTGAGATTGTATCGTTTGGTCTCGATCTATACTGc                                                                                                                             >  2:913114/1‑139 (MQ=255)          aaaGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGTGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGc                                                                                                                      <  1:913114/139‑1 (MQ=255)            aGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTg                                                                                                                    >  2:980135/1‑139 (MQ=255)                 tACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCg                                                                                                               >  2:1003943/1‑139 (MQ=255)                   cTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCgtg                                                                                                             <  1:828081/139‑1 (MQ=255)                        aaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACt                                                                                                        <  1:1003943/139‑1 (MQ=255)                           acTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTa                                                                                                          >  2:1504/1‑134 (MQ=255)                           acTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTa                                                                                                          <  1:1504/134‑1 (MQ=255)                           acTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGaa                                                                                                     >  2:287236/1‑139 (MQ=255)                                       aGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCa                                                                                                                            >  2:85787/1‑104 (MQ=255)                                       aGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCa                                                                                                                            <  1:85787/104‑1 (MQ=255)                                         ccTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCatgatg                                                                                                                       <  1:48724/107‑1 (MQ=255)                                         ccTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCatgatg                                                                                                                       >  2:48724/1‑107 (MQ=255)                                           tATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCatgatg                                                                                                                       >  2:343164/1‑105 (MQ=255)                                           tATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCatgatg                                                                                                                       <  1:343164/105‑1 (MQ=255)                                                    gAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTa                                                                            <  1:287236/139‑1 (MQ=255)                                                       tACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTa                                                                                                          >  2:483607/1‑106 (MQ=255)                                                       tACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTa                                                                                                          <  1:483607/106‑1 (MQ=255)                                                       tACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTAccc                                                                         >  2:93933/1‑139 (MQ=255)                                                        aCCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCg                                                                        >  2:840102/1‑139 (MQ=255)                                                          ccGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTc                                                                      <  1:352769/139‑1 (MQ=255)                                                          ccGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTc                                                                      <  1:942385/139‑1 (MQ=255)                                                                      ctgctgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTAc                                                          <  1:980135/139‑1 (MQ=255)                                                                         ctgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGATGGCGATTGa                                                                                    >  2:593836/1‑110 (MQ=255)                                                                         ctgGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACCGAATACCTGATGGCGATTGa                                                                                    <  1:593836/110‑1 (MQ=255)                                                                                aaGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCtt                                                <  1:840102/139‑1 (MQ=255)                                                                                      tCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGAc                                                                                             >  1:534338/1‑88 (MQ=255)                                                                                      tCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGAc                                                                                             <  2:534338/88‑1 (MQ=255)                                                                                      tCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGt                                                                       <  1:818331/110‑1 (MQ=255)                                                                                      tCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGt                                                                       >  2:818331/1‑110 (MQ=255)                                                                                               gTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTa                                                                            <  1:92646/96‑1 (MQ=255)                                                                                               gTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTa                                                                            >  2:92646/1‑96 (MQ=255)                                                                                                tttGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGc                                >  2:34024/1‑139 (MQ=255)                                                                                                  tGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCgtgt                                                                                                            >  2:370392/1‑61 (MQ=255)                                                                                                  tGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCgtgt                                                                                                            <  1:370392/61‑1 (MQ=255)                                                                                                           gCGAGATTGTATCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAg                     <  1:34024/139‑1 (MQ=255)                                                                                                                      tCGTTTGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAAcgcg          <  2:452791/139‑1 (MQ=255)                                                                                                                         tttGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGa                                                                                 >  2:179829/1‑65 (MQ=255)                                                                                                                         tttGGTCTCGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGa                                                                                 <  1:179829/65‑1 (MQ=255)                                                                                                                              tctcGATCTATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTa  >  2:917682/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                   TCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCGTTTGCACGACGGCGAGATTGTATCGTTTGGTCTCGATCCATACTGCATGATGCTGGAGCGTGTTACTGAATACCTGACGGCGATTGAAGATTTTACCCGTCTGGATTTAGTACGTCGCTGCTTCTATTTAAAAGTGTGCGAAAAGCTCAGCCGTGAACGCGCCTGCGTA  >  NC_000913/3991921‑3992184

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