Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,310,694 C→T D159N (GAC→AAC)  pnp ← polynucleotide phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,310,6940CT100.0% 116.1 / NA 32D159N (GAC→AAC) pnppolynucleotide phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/16);  total (16/16)

CTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGC  >  NC_000913/3310619‑3310769
                                                                           |                                                                           
cTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGa                                                                        <  2:558950/81‑1 (MQ=255)
    aGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGt                                                                     >  2:559001/1‑80 (MQ=255)
        ccGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAAcc                                                                 <  1:558952/80‑1 (MQ=255)
             aacaacCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTAcgcg                                                           <  2:558952/81‑1 (MQ=255)
                aacCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGgc                                                        <  1:558954/81‑1 (MQ=255)
                   cAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGgcagc                                                     <  2:558954/81‑1 (MQ=255)
                      aTCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCAc                                                   >  1:558999/1‑80 (MQ=255)
                       tCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCAcc                                                  >  1:559016/1‑80 (MQ=255)
                          aGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATc                                              <  1:558956/81‑1 (MQ=255)
                              tGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGc                                          >  1:559018/1‑81 (MQ=255)
                                     tttCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCAtt                                      <  2:558956/78‑1 (MQ=255)
                                       tCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGa                                    <  1:558958/78‑1 (MQ=255)
                                        cAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACg                                 >  1:559003/1‑80 (MQ=255)
                                             cGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATAc                           >  2:559021/1‑81 (MQ=255)
                                             cGTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATAc                           >  2:559016/1‑81 (MQ=255)
                                              gTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATAcc                          <  2:558958/81‑1 (MQ=255)
                                              gTCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATAc                           >  1:559010/1‑80 (MQ=255)
                                               tCCTGAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATAc                           >  1:559026/1‑79 (MQ=255)
                                                   gAGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACcaga                       >  2:559008/1‑79 (MQ=255)
                                                    aGTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACcagacag                    <  2:558959/81‑1 (MQ=255)
                                                     gTCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACcagaca                     >  1:559014/1‑79 (MQ=255)
                                                      tCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCagacagag                  <  2:558960/81‑1 (MQ=255)
                                                      tCGGGTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCagacagag                  <  2:558961/81‑1 (MQ=255)
                                                         ggTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCagacagagaca               <  1:558963/81‑1 (MQ=255)
                                                          gTTCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCagacagagacag              >  1:559013/1‑81 (MQ=255)
                                                           ttCAGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAgc             >  1:559011/1‑81 (MQ=255)
                                                              aGTACGTACTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCt          >  2:559024/1‑81 (MQ=255)
                                                               gtacgtacTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTg         >  2:559020/1‑81 (MQ=255)
                                                                tacgtacTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGc        <  1:558964/81‑1 (MQ=255)
                                                                tacgtacTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGc        >  1:559012/1‑81 (MQ=255)
                                                                      cTGGTTATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGc  <  2:558965/81‑1 (MQ=255)
                                                                      cTGGTTATTGATGTAACCTACGCGGGCAGAACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAg   <  2:558964/80‑1 (MQ=255)
                                                                           |                                                                           
CTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGC  >  NC_000913/3310619‑3310769

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