Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,951,007	G→A	A485A (GCC→GCT)	recD ←	exodeoxyribonuclease V subunit RecD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,951,007	0	G	A	100.0%	175.5 / NA	50	A485A (GCC→GCT)	recD	exodeoxyribonuclease V subunit RecD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/20);  total (30/20)

TCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGCCCTTCCCGCAGGGCGCACA  >  NC_000913/2950870‑2951140                                                                                                                                          |                                                                                                                                      tCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGTCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGc                                                                                                                                     >  1:693110/1‑140 (MQ=255)    tGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAat                                                                                                                                  <  1:866514/140‑1 (MQ=255)    tGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAat                                                                                                                                  >  2:866514/1‑140 (MQ=255)     gCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACTAGCaa                                                                                                                                   >  2:85180/1‑138 (MQ=255)        ctGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAatcatc                                                                                                                              <  1:712152/140‑1 (MQ=255)                   cTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGAcc                                                                                                                   >  2:884324/1‑140 (MQ=255)                    tCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCt                                                                                                                  <  2:164376/140‑1 (MQ=255)                     cGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCtt                                                                                                                 >  2:235023/1‑140 (MQ=255)                      ggTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTc                                                                                                                <  1:884324/140‑1 (MQ=255)                              aGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAAc                                                                                                        <  1:235023/140‑1 (MQ=255)                                          gtcGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGAt                                                                                            >  1:723731/3‑140 (MQ=255)                                          gccggccgGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGGCGACCTTCGTACCAACGAGAGTGCGGAt                                                                                            >  2:618434/1‑140 (MQ=255)                                          gCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATGTCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGAt                                                                                            >  1:1001081/1‑140 (MQ=255)                                          gCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTATGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGAt                                                                                            >  2:480890/1‑140 (MQ=255)                                          gCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGac                                                                                            >  1:814045/1‑139 (MQ=255)                                          gCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGAt                                                                                            >  2:208952/1‑140 (MQ=255)                                          gCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGAt                                                                                            >  1:105754/1‑140 (MQ=255)                                                        cAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTg                                                                              <  2:566952/140‑1 (MQ=255)                                                               gACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTATGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCt                                                                       <  1:480890/140‑1 (MQ=255)                                                               gACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCt                                                                       <  1:158129/140‑1 (MQ=255)                                                               gACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCt                                                                       <  2:1001081/140‑1 (MQ=255)                                                               gACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGACGACGATGAATTTTGCGCTTCt                                                                       <  2:814045/140‑1 (MQ=255)                                                                aCGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTg                                                                      >  1:320432/1‑140 (MQ=255)                                                                    gcgTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGc                                                                  <  2:1012290/140‑1 (MQ=255)                                                                      gTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCAt                                                                >  2:1021193/1‑140 (MQ=255)                                                                         cccTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATaaa                                                             <  1:1021193/140‑1 (MQ=255)                                                                              cccGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCt                                                                       >  1:88135/1‑125 (MQ=255)                                                                              cccGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCt                                                                       <  2:88135/125‑1 (MQ=255)                                                                              cccGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCt                                                        >  1:373996/1‑140 (MQ=255)                                                                                cgcgATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACTATGAATTTTGCGCTTCTGTTGCATAAACTGCTCa                                                      >  1:633093/1‑140 (MQ=255)                                                                                  cgATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGc                                                         >  2:799404/1‑135 (MQ=255)                                                                                  cgATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGc                                                         <  1:799404/135‑1 (MQ=255)                                                                                      ccAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGc                                                <  2:373996/140‑1 (MQ=255)                                                                                           gcAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCAtt                                           <  1:101445/140‑1 (MQ=255)                                                                                           gcAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCAtt                                           <  2:723731/140‑1 (MQ=255)                                                                                               tACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAGCTGCTCAATTCGCTCATTCAGt                                       >  1:765843/1‑140 (MQ=255)                                                                                               tACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCGTTCAGt                                       >  1:654374/1‑140 (MQ=255)                                                                                               tACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGt                                       >  2:701682/1‑140 (MQ=255)                                                                                               tACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGt                                       >  2:71797/1‑140 (MQ=255)                                                                                               tACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGt                                       >  1:94527/1‑140 (MQ=255)                                                                                               tACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGt                                       >  1:271101/1‑140 (MQ=255)                                                                                                       tCGCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCAc                               <  2:693110/140‑1 (MQ=255)                                                                                                         gCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGc                             >  1:202267/1‑140 (MQ=255)                                                                                                         gCCATTAAACAACCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGc                             >  1:19255/1‑140 (MQ=255)                                                                                                                 acaacCCAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGc                     <  1:701682/140‑1 (MQ=255)                                                                                                                     cccAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGCCCtt                 >  1:618309/1‑140 (MQ=255)                                                                                                                     cccAAGCGCGCTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGCCCtt                 >  1:366089/1‑140 (MQ=255)                                                                                                                               cTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGCCCTTCCCGCAGGgc       >  2:1055360/1‑140 (MQ=255)                                                                                                                               cTGTCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGCCCTTCCCGCAGGgc       >  2:884/1‑140 (MQ=255)                                                                                                                                  tCATTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGCCCTTCCCGCAGGGCGca    >  1:811557/1‑140 (MQ=255)                                                                                                                                    aTTACGAGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGCCCTTCCCGTAGGGCGcaca  <  1:961226/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      TCGTGCTCTGGCAAGCGACTCGGTTGCACAGACTTAATATTGCCGTCCGGCATCGCAAACCAGACGCGCGTCCCCTGCCCGCGATCCAGCGCAATACCGATATCGCCATTAAACAACCCAAGCGCGCTGTCATTACGGGCAATCATCACCGGTCGACCTTCGTACCAACGAGAGTGCGGATGACGATGAATTTTGCGCTTCTGTTGCATAAACTGCTCAATTCGCTCATTCAGTCCAGCCACGCCAAACGGCCCTTCCCGCAGGGCGCACA  >  NC_000913/2950870‑2951140

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: ‑