Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 308,082 Δ1 bp coding (1251/2526 nt) ecpC ← putative fimbrial usher protein EcpC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913308,0800C.97.6% 164.8 / ‑3.1 41coding (1253/2526 nt)ecpCputative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (18/22);  minor base A (0/1);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGG  >  NC_000913/308003‑308152
                                                                             |                                                                        
aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccac                                                                       >  1:71218/1‑79 (MQ=255)
aTGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccac                                                                       >  1:71198/1‑79 (MQ=255)
  gAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTcaccaccg                                                                     >  1:71207/1‑77 (MQ=255)
    gCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGcc                                                                  <  2:71161/81‑1 (MQ=255)
     cTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCt                                                                 <  2:71162/81‑1 (MQ=255)
      tGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCt                                                                 >  1:71220/1‑80 (MQ=255)
        tCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGAt                                                              >  1:71225/1‑81 (MQ=255)
          aCTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGAtt                                                             <  2:71163/80‑1 (MQ=255)
          aCTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTa                                                            >  1:71219/1‑81 (MQ=255)
           cTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTAt                                                           <  1:71164/81‑1 (MQ=255)
              gCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCg                                                         <  2:71165/80‑1 (MQ=255)
                cAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTAt                                                      >  1:71229/1‑81 (MQ=255)
                cAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTAt                                                      >  1:71223/1‑81 (MQ=255)
                  gCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATcc                                                    <  2:71166/81‑1 (MQ=255)
                     tatTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCata                                                 <  1:71168/81‑1 (MQ=255)
                     tatTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCata                                                 <  1:71167/81‑1 (MQ=255)
                            cAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGt                                           <  2:71169/80‑1 (MQ=255)
                                     gTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAg                                 >  1:71221/1‑81 (MQ=255)
                                      ttGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGc                                <  2:71170/81‑1 (MQ=255)
                                      ttGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGc                                >  2:71241/1‑81 (MQ=255)
                                        gATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCt                               <  2:71171/80‑1 (MQ=255)
                                         aTCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTaa                             <  2:71172/81‑1 (MQ=255)
                                            gcacccccAAGCGGCAGCGTCAGACGGGTTTAACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGc                           <  1:71173/77‑1 (MQ=255)
                                            gCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCg                          <  2:71174/81‑1 (MQ=255)
                                               cccccAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTAc                       >  2:71239/1‑81 (MQ=255)
                                                 cccAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACt                      >  1:71235/1‑80 (MQ=255)
                                                  ccAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTc                     >  1:71232/1‑80 (MQ=255)
                                                  ccAAGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCa                    <  2:71175/81‑1 (MQ=255)
                                                    aaGCGGCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAAt                  >  2:71244/1‑81 (MQ=255)
                                                        ggCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc              >  2:71208/1‑81 (MQ=255)
                                                        ggCAGCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc              <  1:71176/81‑1 (MQ=255)
                                                            gCGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGa          <  2:71177/81‑1 (MQ=255)
                                                             cGTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAg         <  2:71178/81‑1 (MQ=255)
                                                              gtgaGACGGGTTTCACCACCGACTGAGTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAgg        >  2:71233/4‑81 (MQ=255)
                                                              gTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAgg        >  1:71237/1‑81 (MQ=255)
                                                              gTCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAgg        <  1:71179/81‑1 (MQ=255)
                                                               tCAGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGt       <  1:71180/81‑1 (MQ=255)
                                                                 aGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTc      <  1:71183/80‑1 (MQ=255)
                                                                 aGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCg     <  1:71182/81‑1 (MQ=255)
                                                                 aGACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCg     <  2:71181/81‑1 (MQ=255)
                                                                  gACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTc      >  2:71228/1‑79 (MQ=255)
                                                                  gACGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTc      >  2:71231/1‑79 (MQ=255)
                                                                    cGGGTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAgg  >  1:71271/1‑81 (MQ=255)
                                                                     gggTTTCACCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAgg  >  1:71236/1‑80 (MQ=255)
                                                                             |                                                                        
ATGAGCTGTCACTGGCCAGCATATTTTGCAGGTTAACGTTGATCGCCCCCCCAAGCGGCAGCGTCAGACGGGTTTCACCCACCGCCTGATTATCGTATCCATATCCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGG  >  NC_000913/308003‑308152

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