Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,271,455	G→A	T132I (ACC→ATC)	chaA ←	calcium/sodium:proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,271,455	0	G	A	100.0%	159.8 / NA	50	T132I (ACC→ATC)	chaA	calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/23);  total (27/23)

TTAATGCTACCAGTAACGCCTGACCGGTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGGTGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATCATAATGATTGAATAGAGCGTATCACGCATTAGCGTTGGCGCGGCGTCGCCGGTTGCCATTAAAGCTG  >  NC_000913/1271321‑1271581                                                                                                                                       |                                                                                                                               ttAATGCTACCAGTAACGCCTGACCGGTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGc                                                                                                                            <  2:1065899/139‑1 (MQ=255)          ccAGTAACGCCTGACCGGTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACg                                                                                                                   >  1:368999/1‑139 (MQ=255)             gTAACGCCTGACCGGTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGAcc                                                                                                                >  1:1550974/1‑139 (MQ=255)                 cGCCTGACCGGTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGccc                                                                                                            <  2:368999/139‑1 (MQ=255)                    cTGACCGGTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCcaac                                                                                                         >  1:697212/1‑139 (MQ=255)                        ccGGTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGccc                                                                                                            <  2:1579393/132‑1 (MQ=255)                        ccGGTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGccc                                                                                                            >  1:1579393/1‑132 (MQ=255)                        ccGGTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACaata                                                                                                     >  2:1521832/1‑139 (MQ=255)                          ggTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACaataat                                                                                                   <  2:1030878/139‑1 (MQ=255)                            ttGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGa                                                                                                 >  1:908757/1‑139 (MQ=255)                                      cgccgcAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCg                                                                                                               <  1:225983/115‑1 (MQ=255)                                      cgccgcAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCg                                                                                                               >  2:225983/1‑115 (MQ=255)                                                     cATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATAGATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACa                                                                        <  2:908757/139‑1 (MQ=255)                                                        tGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGt                                                                            <  2:861627/132‑1 (MQ=255)                                                        tGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGt                                                                            >  1:861627/1‑132 (MQ=255)                                                        tGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATc                                                                     >  2:1424448/1‑139 (MQ=255)                                                              tACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATg                                                                                                  <  1:269428/104‑1 (MQ=255)                                                              tACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATg                                                                                                  >  2:269428/1‑104 (MQ=255)                                                                    tACGATTATTGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGt                                                                            <  2:1000616/120‑1 (MQ=255)                                                                    tACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATCAAAAGCCAACCAGCCCACCGGt                                                                            >  1:1000616/1‑120 (MQ=255)                                                                         ttATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAAc                                                                                        <  2:937407/103‑1 (MQ=255)                                                                         ttATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAAc                                                                                        >  1:937407/1‑103 (MQ=255)                                                                         ttATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGt                                                                            >  2:407097/1‑115 (MQ=255)                                                                         ttATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGt                                                                            <  1:407097/115‑1 (MQ=255)                                                                                   ggggAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGa                                                                                                                  <  1:793381/67‑1 (MQ=255)                                                                                   ggggAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGa                                                                                                                  >  2:793381/1‑67 (MQ=255)                                                                                   ggggAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGcc                                                                                   <  2:1114440/98‑1 (MQ=255)                                                                                   ggggAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGcc                                                                                   >  1:1114440/1‑98 (MQ=255)                                                                                   ggggAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGt                                                                            >  1:452634/1‑105 (MQ=255)                                                                                   ggggAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGt                                                                            <  2:452634/105‑1 (MQ=255)                                                                                   ggggAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATCATAATGATTGAATa                                                       <  2:138081/126‑1 (MQ=255)                                                                                   ggggAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATCATAATGATTGAATa                                                       >  1:138081/1‑126 (MQ=255)                                                                                            cgcAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGc                                                                                                                            <  2:1228921/48‑1 (MQ=255)                                                                                            cgcAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGc                                                                                                                            >  1:1228921/1‑48 (MQ=255)                                                                                             gcAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAAc                                                                                        >  2:478705/1‑83 (MQ=255)                                                                                             gcAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAAc                                                                                        <  1:478705/83‑1 (MQ=255)                                                                                                    aaTACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGc                                                                                    >  1:1457495/1‑80 (MQ=255)                                                                                                    aaTACTGCTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGc                                                                                    <  2:1457495/80‑1 (MQ=255)                                                                                                           cTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGt                                                                            <  1:517998/81‑1 (MQ=255)                                                                                                           cTTGATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGt                                                                            >  2:517998/1‑81 (MQ=255)                                                                                                              gATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCAc                                                                                >  1:1298505/1‑74 (MQ=255)                                                                                                              gATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCAc                                                                                <  2:1298505/74‑1 (MQ=255)                                                                                                              gATACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATCATAATGATTGAATAGAGCGTATCACGCATTAGCGTTGGCGCGGCGTCGCCGGtt               <  1:1521832/139‑1 (MQ=255)                                                                                                                tACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGaa                                                                                                <  1:91632/56‑1 (MQ=255)                                                                                                                tACCAAACAGATTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGaa                                                                                                >  2:91632/1‑56 (MQ=255)                                                                                                                          aTTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATCATAATGATTGAATagag                                                    >  2:1096178/1‑90 (MQ=255)                                                                                                                          aTTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATCATAATGATTGAATagag                                                    <  1:1096178/90‑1 (MQ=255)                                                                                                                          aTTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATCATAATGATTGAATAGAGCGTATCACGCATTAGCGTTGGCGCGGCGTCGCCGGTTGCCATTAAAGCt   >  2:1437920/1‑139 (MQ=255)                                                                                                                          aTTCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATCATAATGATTGAATAGAGCGTATCACGCATTAGCGTTGGCGCGGCGTCGCCGGTTGCCATTAAAGCt   >  2:1375251/1‑139 (MQ=255)                                                                                                                           ttCATATATTGGATGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATCATAATGATTGAATAGAGCGTATCACGCATTAGCGTTGGCGCGGCGTCGCCGGTTGCCATTAAAGCTg  >  2:1473214/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                               TTAATGCTACCAGTAACGCCTGACCGGTTGAAAAATTCGCCGCAGGCAGAGCCATTGGAAATACCAGTACGATTATCGCCAGGGGGAACAGCGCAATTAAATACTGCTTGATACCAAACAGATTCATATATTGGGTGGCAAACTTACGACCGCCCAACAATAATGAAAAGCCAACCAGCCCACCGGTAACAATCATAATGATTGAATAGAGCGTATCACGCATTAGCGTTGGCGCGGCGTCGCCGGTTGCCATTAAAGCTG  >  NC_000913/1271321‑1271581

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