Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 307,452 T→C G627G (GGA→GGG ecpC ← putative fimbrial usher protein EcpC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913307,4520TC100.0% 134.2 / NA 37G627G (GGA→GGGecpCputative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/17);  total (20/17)

ACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTTCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGTTTTGTCATCACCGGTATCGCCG  >  NC_000913/307374‑307524
                                                                              |                                                                        
aCCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTccc                                                                        >  2:69346/1‑81 (MQ=255)
aCCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTccc                                                                        >  2:69345/1‑81 (MQ=255)
 cccTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTccc                                                                        <  1:69294/80‑1 (MQ=255)
   cTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCt                                                                     <  2:69295/81‑1 (MQ=255)
        aGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCg                                                                 >  1:69365/1‑80 (MQ=255)
           cGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAAc                                                             <  1:69296/81‑1 (MQ=255)
           cGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAAc                                                             <  1:69297/81‑1 (MQ=255)
            gACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAAc                                                             >  2:69357/1‑80 (MQ=255)
            gACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAAc                                                             >  1:69349/1‑80 (MQ=255)
            gACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACg                                                            <  2:69298/81‑1 (MQ=255)
               gCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGc                                                         <  1:69299/81‑1 (MQ=255)
                  gCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTc                                                      >  1:69361/1‑81 (MQ=255)
                     aTTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGaa                                                   >  1:69351/1‑81 (MQ=255)
                      ttGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAAt                                                  <  1:69301/81‑1 (MQ=255)
                      ttGGCGGTCAAGTTGGAATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAAt                                                  <  1:69300/81‑1 (MQ=255)
                       tGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAAtt                                                 >  1:69356/1‑81 (MQ=255)
                         gCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATtg                                                <  2:69302/80‑1 (MQ=255)
                           ggTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATtgtg                                              >  2:69377/1‑80 (MQ=255)
                             tCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCAt                                           <  1:69303/81‑1 (MQ=255)
                               aaGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATAc                                         >  1:69372/1‑81 (MQ=255)
                                   tGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGccc                                     >  2:69370/1‑81 (MQ=255)
                                    ggTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGccc                                     <  1:69304/80‑1 (MQ=255)
                                     gTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCccg                                   <  2:69305/81‑1 (MQ=255)
                                        ttGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCccgccg                                >  2:69353/1‑81 (MQ=255)
                                        ttGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCccgcc                                 <  1:69306/80‑1 (MQ=255)
                                           aCGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCt                              <  2:69307/80‑1 (MQ=255)
                                              tAGCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAgg                          >  2:69350/1‑81 (MQ=255)
                                                gCCGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGt                        >  2:69348/1‑81 (MQ=255)
                                                 ccGTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGtt                       <  1:69308/81‑1 (MQ=255)
                                                   gTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGtttt                     >  2:69355/1‑81 (MQ=255)
                                                   gTCCGCCGCGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGtttt                     >  2:69352/1‑81 (MQ=255)
                                                       gccgcGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGTTTTGtca                 >  2:69383/1‑81 (MQ=255)
                                                       gccgcGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGTTTTGtca                 <  1:69309/81‑1 (MQ=255)
                                                       gccgcGCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCACCGCCGCTGAGGGTTTTGtca                 >  1:69360/1‑81 (MQ=255)
                                                         cgcgCTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGTTTTGtcatc               <  1:69310/81‑1 (MQ=255)
                                                          gcgcTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGTTTTGtcatca              >  1:69359/1‑81 (MQ=255)
                                                            gcTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGTTTTGTCATCAcc            >  1:69364/1‑81 (MQ=255)
                                                             cTATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGTTTTGTCATCACCg           <  1:69311/81‑1 (MQ=255)
                                                              tATTGACGTTCAGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGTTTTGTCATCACCgg          >  1:69381/1‑81 (MQ=255)
                                                                         aGCGTCCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGTTTTGTCATCACCGGTATCGCCg  >  2:69379/1‑78 (MQ=255)
                                                                              |                                                                        
ACCCTGCCAGCCGACGCTGCCATTGGCGGTCAAGTTGGTATTGACGTAGCCGTCCGCCGCGCTATTGACGTTCAGCGTTCCGCTGGCGTAACGAGCGTCGAATTGTGCATACGCCCCGCCGCTGAGGGTTTTGTCATCACCGGTATCGCCG  >  NC_000913/307374‑307524

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