Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 554,427 A→G intergenic (+217/‑3) fimA → / → fimC type 1 fimbrial major subunit FimA/type 1 fimbria chaperone FimC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273554,4270AG100.0% 67.6 / NA 21intergenic (+217/‑3)fimA/fimCtype 1 fimbrial major subunit FimA/type 1 fimbria chaperone FimC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (6/15);  total (6/15)

TGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGAT  >  NZ_CP009273/554349‑554509
                                                                              |                                                                                  
tgatgTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACt                                                                         <  2:234711/90‑1 (MQ=255)
     tgttgtAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAAt                                                                    >  2:488989/1‑90 (MQ=255)
           aaGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGtt                                                              >  1:251166/1‑90 (MQ=255)
              aaaCCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTAtt                                                           <  1:328517/90‑1 (MQ=255)
                aCCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGa                                                         >  1:211713/1‑90 (MQ=255)
                aCCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGa                                                         <  1:392371/90‑1 (MQ=255)
                aCCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGa                                                         <  2:22591/90‑1 (MQ=255)
                    aaGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCt                                                     <  2:457738/90‑1 (MQ=255)
                                      atTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAAt                                   <  1:426016/90‑1 (MQ=255)
                                      atTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAAt                                   <  2:277605/90‑1 (MQ=255)
                                          cTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCAtt                               >  2:35595/1‑90 (MQ=255)
                                                tttGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGcc                         <  1:67279/90‑1 (MQ=255)
                                                 ttGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCa                        <  2:248089/90‑1 (MQ=255)
                                                         aGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCAt                <  2:7038/90‑1 (MQ=255)
                                                         aGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCAt                <  2:417025/90‑1 (MQ=255)
                                                                cacaAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCg         >  1:486405/1‑90 (MQ=255)
                                                                  caAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGa       <  2:481646/90‑1 (MQ=255)
                                                                    aGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGc                                                      <  2:23455/41‑1 (MQ=255)
                                                                    aGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGc                                                      >  1:23455/1‑41 (MQ=255)
                                                                     ggCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAggg    <  2:167580/90‑1 (MQ=255)
                                                                       cGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGAt  <  2:211713/90‑1 (MQ=255)
                                                                              |                                                                                  
TGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGAT  >  NZ_CP009273/554349‑554509

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