Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,820,029 A→G V381A (GTG→GCG)  ligB ← DNA ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,820,0290AG100.0% 93.0 / NA 27V381A (GTG→GCG) ligBDNA ligase B
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/12);  total (15/12)

AGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACT  >  NC_000913/3819953‑3820104
                                                                            |                                                                           
aGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGcgccacgcca                                                                          <  1:684796/80‑1 (MQ=255)
   aGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCAcatc                                                                      <  1:684797/81‑1 (MQ=255)
                 tCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGaa                                                        >  2:684833/1‑81 (MQ=255)
                        gTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCt                                                  >  2:684842/1‑80 (MQ=255)
                        gTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCt                                                  >  1:684825/1‑80 (MQ=255)
                              gCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCgg                                            >  2:684852/1‑80 (MQ=255)
                                ggTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGc                                           >  2:684844/1‑79 (MQ=255)
                                ggTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCaa                                         <  1:684798/81‑1 (MQ=255)
                                     tCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCt                                     <  2:684799/80‑1 (MQ=255)
                                      cTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGa                                   >  2:684831/1‑81 (MQ=255)
                                      cTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGa                                   >  1:684839/1‑81 (MQ=255)
                                      cTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGa                                   >  1:684847/1‑81 (MQ=255)
                                         gCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGAc                                  <  2:684801/79‑1 (MQ=255)
                                         gCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGAc                                  <  2:684802/79‑1 (MQ=255)
                                         gCGGTGTCGGTTTTGTACGTTCTGCACCACGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGAc                                  <  2:684800/79‑1 (MQ=255)
                                          cGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACgag                               <  2:684803/81‑1 (MQ=255)
                                           ggTGTCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACgag                               >  2:684830/1‑80 (MQ=255)
                                            gtgtCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGaa                             >  2:684840/1‑81 (MQ=255)
                                              gtCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATc                           >  2:684836/1‑81 (MQ=255)
                                              gtCGGTTTTGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATc                           >  2:684851/1‑81 (MQ=255)
                                                   ttttGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATc                      <  1:684804/81‑1 (MQ=255)
                                                    tttGTACGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACTCTGACCGGCAAGGCTGACGAGAATCTGATCa                     >  2:684846/1‑81 (MQ=255)
                                                          cGTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAgg                <  2:684805/80‑1 (MQ=255)
                                                           gTTCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGgcg              <  1:684806/81‑1 (MQ=255)
                                                             tCTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGgcgc             >  1:684837/1‑80 (MQ=255)
                                                              cTGCACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCa            <  1:684807/80‑1 (MQ=255)
                                                                 cACCGCGCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAatat        <  1:684808/81‑1 (MQ=255)
                                                                    cgcgCCACGCCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATcc      >  1:684845/1‑80 (MQ=255)
                                                                       gccacgccacATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACt  >  2:684866/1‑81 (MQ=255)
                                                                            |                                                                           
AGAAGCAAAGTAGCAGGTCAACGAGTTAAAGCGGTTTTCTGGCGGTGTCGGTTTTGTACGTTCTGCACCGCGCCACACCACATCATCAATGCGAGGAATACCCTGACCGGCAAGGCTGACGAGAATCTGATCACCAGGCGCAATATCCCACT  >  NC_000913/3819953‑3820104

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: