Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,193,328	G→A	H281Y (CAT→TAT)	mrp ←	P‑loop NTPase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,193,328	0	G	A	100.0%	137.5 / NA	40	H281Y (CAT→TAT)	mrp	P‑loop NTPase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/14);  total (26/14)

AGATCCTCACGTAAGGAGATGTGTAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATGCACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCGGGGTGGTGACCACGACTGCACCGGTTACTGGAATG  >  NC_000913/2193201‑2193458                                                                                                                                |                                                                                                                                   agaTCCTCACGTAAGGAGATGTGTAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGtt                                                                                                                        <  1:205779/140‑1 (MQ=255) agaTCCTCACGTAAGGAGATGTGTAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGtt                                                                                                                        <  1:212736/140‑1 (MQ=255) agaTCCTCACGTAAGGAGATGTGTAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGtt                                                                                                                        >  1:872293/1‑140 (MQ=255)   aTCCTCACGTAAGGAGATGTGTAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTtct                                                                                                                      >  1:474531/1‑140 (MQ=255)   aTCCTCACGTAAGGAGATGTGTAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTtct                                                                                                                      >  1:387083/1‑140 (MQ=255)           gTAAGGAGATGTGTAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATa                                                                                                              >  2:256627/1‑140 (MQ=255)                    tgtgTAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCAcc                                                                                                     <  2:297409/140‑1 (MQ=255)                    tgtgTAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCAcc                                                                                                     <  1:251558/140‑1 (MQ=255)                      tgtAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCgg                                                                                                   >  2:76216/1‑140 (MQ=255)                       gtAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGa                                                                                                  <  2:264506/140‑1 (MQ=255)                              cATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGa                                                                                           >  2:270866/1‑140 (MQ=255)                              cATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGa                                                                                           >  2:168877/1‑140 (MQ=255)                                tCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGa                                                                                           >  2:972262/1‑138 (MQ=255)                                tCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGa                                                                                           <  1:972262/138‑1 (MQ=255)                                     cccAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTtctc                                                                                    <  2:313282/140‑1 (MQ=255)                                       cagcagCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGa                                                                                  >  2:919750/1‑140 (MQ=255)                                              tgtgTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATcac                                                                           >  2:89656/1‑140 (MQ=255)                                                       tATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCttt                                                                  >  1:955668/1‑140 (MQ=255)                                                          ttttCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCtttttt                                                               >  1:928040/1‑140 (MQ=255)                                                          ttttCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCtttttt                                                               >  1:116461/1‑140 (MQ=255)                                                          ttttCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCtttttt                                                               >  1:286506/1‑140 (MQ=255)                                                            ttCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCg                                                             >  1:448901/1‑140 (MQ=255)                                                               aGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCAt                                                          >  1:5002/1‑140 (MQ=255)                                                                   aGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGAt                                                      <  2:476269/140‑1 (MQ=255)                                                                     tttCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCa                                                    <  2:448901/140‑1 (MQ=255)                                                                          ctGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGcgcg                                               <  1:919750/140‑1 (MQ=255)                                                                             caccacCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATg                                            <  2:928040/140‑1 (MQ=255)                                                                               ccaccGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTc                                          >  2:285418/1‑140 (MQ=255)                                                                                  ccGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTg                                       >  1:380979/1‑140 (MQ=255)                                                                                  ccGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTg                                       >  2:752861/1‑140 (MQ=255)                                                                                  ccGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTg                                       >  2:25985/1‑140 (MQ=255)                                                                                       aCCGAAAATCTGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCgggg                                  >  1:412453/1‑140 (MQ=255)                                                                                         cGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCGGggtg                                <  2:380979/140‑1 (MQ=255)                                                                                                     cATGGTGACCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCGGGGTGGTGACCACGACt                    >  1:410311/1‑140 (MQ=255)                                                                                                            aCCACAGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCGGGGTGGTGACCACGACTGCACCgg             <  2:712226/140‑1 (MQ=255)                                                                                                               acaGTTACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCGGGGTGGTGACCACGACTGCACCGGTTa          <  1:25985/140‑1 (MQ=255)                                                                                                                   ttACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCGGGGTGGTGACCACGACTGCACCGGTTACTgg      >  2:446354/1‑140 (MQ=255)                                                                                                                   ttACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCGGGGTGGTGACCACGACTGCACCGGTTACTgg      >  2:476420/1‑140 (MQ=255)                                                                                                                   ttACTGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCGGGGTGGTGACCACGACTGCACCGGTTACTgg      >  1:454236/1‑140 (MQ=255)                                                                                                                       tGCAAATATACACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCGGGGTGGTGACCACGACTGCACCGGTTACTGGAATg  >  2:907604/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                   AGATCCTCACGTAAGGAGATGTGTAACGGCATCTGCCCCAGCAGCTGTGTATGGTATTTTTCAGCCAGTTTCTCTGCACCACCGGTACCGAAAATCGGTTCATGGTGACCACAGTTACTGCAAATATGCACACTCATGTTCTCGACAATACCGAGCACCGGAACCTCGACTTTCTCGAACATCACAATGCCTTTTTTCGCATCGATCAGCGCGATGTCCTGCGGGGTGGTGACCACGACTGCACCGGTTACTGGAATG  >  NC_000913/2193201‑2193458

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