Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,208,517 C→T P129L (CCT→CTT) 
*201* (TAG→TAA
tfaP →
tfaE ←
e14 prophage, uncharacterized protein
e14 prophage, putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,208,5170CT100.0% 187.8 / NA 62P129L (CCT→CTT) 
*201* (TAG→TAA
tfaP
tfaE
e14 prophage, uncharacterized protein
e14 prophage, putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/33);  total (29/33)

GCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCC  >  NC_000913/1208385‑1208646
                                                                                                                                    |                                                                                                                                 
gCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGcc                                                                                                                             <  2:120969/139‑1 (MQ=255)
    aCAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCg                                                                                                                         <  2:876595/139‑1 (MQ=255)
    aCAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCg                                                                                                                         <  1:762303/139‑1 (MQ=255)
       gTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCa                                                                                                                      >  1:470081/1‑139 (MQ=255)
        tGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAg                                                                                                                     >  2:696617/1‑139 (MQ=255)
         ggCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGt                                                                                                                    >  1:743099/1‑139 (MQ=255)
           cAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGtta                                                                                                                  <  2:470081/139‑2 (MQ=255)
                 aaGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAg                                                                                                           <  2:141538/139‑1 (MQ=255)
                        cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACAt                                                                                                    >  1:631327/1‑139 (MQ=255)
                              gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                              <  2:631327/139‑1 (MQ=255)
                              gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                              <  1:410848/139‑1 (MQ=255)
                                aGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGcgcg                                                                                            <  1:252104/139‑1 (MQ=255)
                                 gACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCgg                                                                                           <  2:1046156/139‑1 (MQ=255)
                                  aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGttaggc                                                                                                               <  1:1008778/119‑5 (MQ=255)
                                  aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGttaggc                                                                                                               >  2:1008778/1‑115 (MQ=255)
                                  aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGt                                                                                          <  2:743099/139‑1 (MQ=255)
                                  aCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGt                                                                                          >  2:990121/1‑139 (MQ=255)
                                      aaaGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTg                                                                                      >  1:1086620/1‑139 (MQ=255)
                                       aagaagATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTgg                                                                                     >  2:1063633/1‑139 (MQ=255)
                                             aTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCg                                                                                                                         >  1:1120258/1‑98 (MQ=255)
                                             aTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCg                                                                                                                         <  2:1120258/98‑1 (MQ=255)
                                                 gCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGc                                                                                                                              <  2:1478987/89‑1 (MQ=255)
                                                 gCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGc                                                                                                                              >  1:1478987/1‑89 (MQ=255)
                                                            ataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                <  2:748869/139‑1 (MQ=255)
                                                            ataATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                >  1:748869/1‑139 (MQ=255)
                                                                   gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTcc                                                                                                                          >  2:524100/1‑75 (MQ=255)
                                                                   gtgGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTcc                                                                                                                          <  1:524100/75‑1 (MQ=255)
                                                                     ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                              >  2:1049859/1‑100 (MQ=255)
                                                                     ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                              <  1:1049859/100‑1 (MQ=255)
                                                                     ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                >  1:1134930/1‑130 (MQ=255)
                                                                     ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt                                                                <  2:1134930/130‑1 (MQ=255)
                                                                     ggCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAAt                                                       >  2:669597/1‑139 (MQ=255)
                                                                      gCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTg                                                                                         >  1:997927/1‑104 (MQ=255)
                                                                      gCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTg                                                                                         <  2:997927/104‑1 (MQ=255)
                                                                      gCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATc                                                      >  2:1307302/1‑139 (MQ=255)
                                                                           aTTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGca                                                 <  1:669597/139‑1 (MQ=255)
                                                                                  ggAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGtt                                          <  1:990121/139‑1 (MQ=255)
                                                                                         cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTAt                                                                                  >  2:764021/1‑92 (MQ=255)
                                                                                         cTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTAt                                                                                  <  1:764021/92‑1 (MQ=255)
                                                                                          tGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTgg                                  <  1:1083626/139‑1 (MQ=255)
                                                                                                     gATACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                       <  1:1063633/139‑1 (MQ=255)
                                                                                                      aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                              <  1:947973/67‑1 (MQ=14)
                                                                                                      aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                              >  2:947973/1‑67 (MQ=14)
                                                                                                      aTACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcct                      >  1:836560/1‑139 (MQ=255)
                                                                                                       tACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                              <  1:366346/66‑1 (MQ=14)
                                                                                                       tACATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                              >  2:366346/1‑66 (MQ=14)
                                                                                                        aCATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCgg                                                                                                <  2:1303251/63‑1 (MQ=14)
                                                                                                        aCATCCAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCgg                                                                                                >  1:1303251/1‑63 (MQ=14)
                                                                                                            ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                >  2:1261984/1‑139 (MQ=255)
                                                                                                            ccAGTGCGCCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                >  2:1103740/1‑139 (MQ=255)
                                                                                                                 gcgcCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTa           >  2:202976/1‑139 (MQ=255)
                                                                                                                  cgcCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTAc          <  2:836560/139‑1 (MQ=255)
                                                                                                                   gcCAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACg         <  1:202976/139‑1 (MQ=255)
                                                                                                                    ccAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                              >  2:548653/1‑53 (MQ=14)
                                                                                                                    ccAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGgcg                                                                                              <  1:548653/53‑1 (MQ=14)
                                                                                                                    ccAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTACCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCg                   <  1:153810/128‑1 (MQ=17)
                                                                                                                    ccAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTACCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCg                   >  2:153810/1‑128 (MQ=17)
                                                                                                                     cAGATATTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc       <  1:1307302/139‑1 (MQ=255)
                                                                                                                        atatTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCgg    <  1:1103740/139‑1 (MQ=255)
                                                                                                                          atTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGt                                                                                          >  2:43627/1‑51 (MQ=14)
                                                                                                                          atTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGt                                                                                          <  1:43627/51‑1 (MQ=14)
                                                                                                                          atTGAATGGCTTACGCCTCCGGCAGTTAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGTCTGCGTCAGTTTACGTCCGGcc  >  2:1524415/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                 
GCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCC  >  NC_000913/1208385‑1208646

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