Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 308,755 C→T S193N (AGC→AAC)  ecpC ← ECP production outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913308,7550CT100.0% 119.8 / NA 33S193N (AGC→AAC) ecpCECP production outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/6);  total (27/6)

CCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGAC  >  NC_000913/308608‑308900
                                                                                                                                                   |                                                                                                                                                 
ccGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTgg                                                                                                                                                >  2:693838/1‑151 (MQ=255)
 cGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTgg                                                                                                                                                >  2:319263/1‑150 (MQ=255)
   gCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGAtg                                                                                                                                             <  2:478614/151‑1 (MQ=255)
               ccGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCccgc                                                                                                                                 >  1:286849/1‑151 (MQ=255)
                 ggCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCccgccg                                                                                                                               >  2:765016/1‑151 (MQ=255)
                    aaaaTCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGtt                                                                                                                             >  2:114009/1‑150 (MQ=255)
                      aaTCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGATCAGTCTAGCACCACATGATGTTCGCGCTGTGCAGTAACGCTGTTCAGCGACGGATAGTTGGATGGGTTGCTCCCGCCGTTACg                                                                                                                          >  1:256542/1‑151 (MQ=255)
                          gcgTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAAc                                                                                                                      >  2:146145/1‑151 (MQ=255)
                           cgTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACt                                                                                                                     >  1:139460/1‑151 (MQ=255)
                            gTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTg                                                                                                                    >  2:200271/1‑151 (MQ=255)
                               cATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGgtt                                                                                                                 >  1:71940/1‑151 (MQ=255)
                                aTACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGgt                                                                                                                  >  2:338269/1‑149 (MQ=255)
                                    aTCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTta                                                                                                             >  1:341933/1‑150 (MQ=255)
                                         tttATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACg                                                                                                       >  2:278247/1‑151 (MQ=255)
                                             tataATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGccc                                                                                                    >  2:233300/1‑150 (MQ=255)
                                                  ttCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGtt                                                                                               >  1:532183/1‑150 (MQ=255)
                                                   tCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTtat                                                                                             >  1:289218/1‑151 (MQ=255)
                                                     ccTGTCCTGTTGCCCGCTACCGCTCCCGTACAGCGGGCCGGCGACCCCCACATGTTGTGCGCGCAGTGCAGTATCTTTATTCAGCGAAAAATAGTGGGTTGTGTTGCCCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAg                                                                                           <  1:1245/150‑1 (MQ=255)
                                                           cTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGAAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTc                                                                                        >  2:632175/1‑148 (MQ=255)
                                                              ttGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGAtt                                                                                   >  2:67864/1‑150 (MQ=255)
                                                                  ccGCTACCTCACCCGTACCGGGAGCTGTCGACCACCACGTGATGTTCGCGCAGTGCATTACCTTTATTGCGCTACATATAGTGGGATTTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTAct                                                                                <  2:238990/149‑1 (MQ=255)
                                                                          gATCCCGAACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCACGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCACAGATTACTGCTGAgg                                                                         >  2:421988/1‑148 (MQ=255)
                                                                                       gAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGGTAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGTTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGGGTTAACACTGGAc                                                          >  1:49475/1‑150 (MQ=255)
                                                                                          ccgccgacccccaattgattTTCGCGCAGTGAAGAAACTTTATTCGGCGAAAAATAGTGGGTTGTTTTCCCCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGcc                                                      <  1:713128/136‑1 (MQ=255)
                                                                                                     cacaTGATGTGCGCGGAGTGCAGTATTGTTATTCAGCGACAGATGTTTGGTTGTTTTCCCCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGAt                                                  <  1:90364/144‑1 (MQ=255)
                                                                                                        atgatgTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCt                                         >  2:280739/1‑150 (MQ=255)
                                                                                                        atgatgTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCt                                         >  1:174321/1‑150 (MQ=255)
                                                                                                                   ggggtgcagtATTGTTATGCGGCGAAAGATAGTGGGATGTTTTGCCCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAg                                  <  1:306537/143‑1 (MQ=255)
                                                                                                                            tAACGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGcc                    >  2:98575/1‑151 (MQ=255)
                                                                                                                             aaCGTTATTCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGccc                   >  2:506418/1‑151 (MQ=255)
                                                                                                                                    ttCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCtt            >  2:213003/1‑151 (MQ=255)
                                                                                                                                     tCAGCGACAGATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTc           >  2:730758/1‑151 (MQ=255)
                                                                                                                                              gATAGTTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGAc  >  2:606722/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGGCAAATCGGTGACCGGCAAAATCGCGTTCATACATCGCTTTATATAATTCACTGTCCTGTTGACCGCTACCGATCCCGTACAGCGAGCCGTCGAGCACCACATGATGTTCGCGCAGTGCAGTAACGTTATTCAGCGACAGATAGCTGGATGTGTTGCTCCCGCCGTTACGCAACTGGTTGTTATAGACGCCCAAGTTATAGCTCAGATTACTGCTGAGGGTGTTAACACTGGACTGCCCGATGTCTTCGCTACGTGAGCGTAGTACGGTGCCCAGCGCTTCGCGCTTGAC  >  NC_000913/308608‑308900

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