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,4270AG97.5% 145.8 / ‑3.5 40intergenic (+217/‑3)fimA/fimCtype 1 fimbrial major subunit FimA/type 1 fimbria chaperone FimC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/19);  minor base C (1/0);  total (21/19)
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.01e-01

AAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGC  >  NZ_CP009273/554340‑554512
                                                                                       |                                                                                     
aaaaGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCt                                                                                     >  1:323001/1‑90 (MQ=255)
    gggAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTatga                                                                                 <  1:48784/90‑1 (MQ=255)
    gggAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTatga                                                                                 <  2:1685325/90‑1 (MQ=255)
        atgatgTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGAc                                                                             <  2:636604/90‑1 (MQ=255)
        atgatgTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGAc                                                                             <  1:2893037/90‑1 (MQ=255)
         tgatgTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACt                                                                            <  2:1805722/90‑1 (MQ=255)
         tgatgTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACt                                                                            <  1:150560/90‑1 (MQ=255)
              tgttgtAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAAt                                                                       >  2:2566409/1‑90 (MQ=255)
              tgttgtAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAAt                                                                       >  1:141740/1‑90 (MQ=255)
                ttgtAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATaa                                                                     >  2:16095/1‑90 (MQ=255)
                  gtAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTACCTATGATGACTAAAATAAAg                                                                   >  2:2844340/1‑90 (MQ=255)
                    aaGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGtt                                                                 >  2:2395156/1‑90 (MQ=255)
                         aCCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGa                                                            <  1:1358421/90‑1 (MQ=255)
                           cAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATg                                                          <  1:2929841/90‑1 (MQ=255)
                            aaaGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGAc                                                                             >  2:2856952/1‑70 (MQ=255)
                            aaaGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGAc                                                                             <  1:2856952/70‑1 (MQ=255)
                             aaGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCt                                                        <  2:2479080/90‑1 (MQ=255)
                             aaGCAATCATTTCTTTATATTCCTCATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCt                                                        <  1:1943562/90‑1 (MQ=255)
                                  aTCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATta                                                   <  1:2275224/90‑1 (MQ=255)
                                              tatTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTaa                                       >  1:263554/1‑90 (MQ=255)
                                               atTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAAt                                      <  2:2010641/90‑1 (MQ=255)
                                               atTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAAt                                      <  1:2397477/90‑1 (MQ=255)
                                                ttCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATc                                     >  2:1888510/1‑90 (MQ=255)
                                                  ccTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCAt                                   >  1:1213364/1‑90 (MQ=255)
                                                    ttATTTTTGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCAttt                                 >  2:1519899/1‑90 (MQ=255)
                                                         tttGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGcc                            >  2:468170/1‑90 (MQ=255)
                                                          ttGCCGTCAGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCa                           >  2:1272885/1‑90 (MQ=255)
                                                                  aGGAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCAt                   <  1:1224681/90‑1 (MQ=255)
                                                                    gAATGCACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCAtgc                 <  2:2683496/90‑1 (MQ=255)
                                                                    gAATACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCAtgc                 <  1:1662529/90‑1 (MQ=255)
                                                                       tACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCAtgctgc              >  2:1782127/1‑90 (MQ=255)
                                                                       tACACAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCAtgctgc              >  2:1117594/1‑90 (MQ=255)
                                                                        acacAAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGcc             <  1:2795998/90‑1 (MQ=255)
                                                                           caAGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGa          <  2:2744050/90‑1 (MQ=255)
                                                                             aGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAgg        >  2:1140952/1‑90 (MQ=255)
                                                                             aGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAgg        >  2:2748209/1‑90 (MQ=255)
                                                                             aGGCGTATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAgg        >  2:804315/1‑90 (MQ=255)
                                                                                  tATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCg   >  1:2829938/1‑90 (MQ=255)
                                                                                  tATTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCg   >  1:1499399/1‑90 (MQ=255)
                                                                                   aTTAGCTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGc  >  2:973771/1‑90 (MQ=255)
                                                                                       |                                                                                     
AAAAGGGAATGATGTGTTGTAAGAAACCAAAGCAATCATTTCTTTATATTCCTTATTTTTGCCGTCAGGAATACACAAGGCGTATTAACTATGATGACTAAAATAAAGTTATTGATGCTCATTATATTTTATTTAATCATTTCGGCCAGCGCCCATGCTGCCGGAGGGATCGC  >  NZ_CP009273/554340‑554512

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