Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 149,183 A→G C39C (TGT→TGC)  yadM ← fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273149,1830AG97.9% 177.0 / ‑3.7 48C39C (TGT→TGC) yadMfimbrial protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (22/25);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

AAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGTAG  >  NZ_CP009273/149064‑149315
                                                                                                                       |                                                                                                                                    
aaaTCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAg                                                                                                                   <  1:757343/139‑1 (MQ=255)
      gCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTg                                                                                                             >  1:1158736/1‑139 (MQ=255)
        cTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCAACATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTa                                                                                                           <  1:478484/139‑1 (MQ=255)
            gaTTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAAct                                                                                                       <  2:155064/139‑1 (MQ=255)
                cTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTa                                                                                                    <  1:566710/138‑1 (MQ=255)
                             cgcgATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGt                                                                                       <  2:610718/138‑1 (MQ=255)
                             cgcgATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGt                                                                                       <  1:420388/138‑1 (MQ=255)
                             cgcgATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGt                                                                                       <  2:515335/138‑1 (MQ=255)
                                      aTCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCa                                                                                   <  2:981247/133‑1 (MQ=255)
                                       tCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACt                                                                             <  2:472103/138‑1 (MQ=255)
                                       tCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTg                                                                            <  1:777900/139‑1 (MQ=255)
                                          aGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGATGTAAACTCTACATTAATATCCGTTGCAAATACTGAgg                                                                         <  2:43873/139‑1 (MQ=255)
                                            aCCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGt                                                                                                                              >  1:1118749/1‑84 (MQ=255)
                                            aCCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGt                                                                                                                              <  2:1118749/84‑1 (MQ=255)
                                            aCCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGt                                                                       >  2:1139799/1‑139 (MQ=255)
                                            aCCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGt                                                                       >  2:919988/1‑139 (MQ=255)
                                                aTCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGat                                                                    >  2:49205/1‑138 (MQ=255)
                                                aTCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGat                                                                    >  2:534877/1‑138 (MQ=255)
                                                aTCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGat                                                                    >  2:296994/1‑138 (MQ=255)
                                                      gggATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataata                                                             <  2:87177/139‑1 (MQ=255)
                                                           tctcAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAAct                                                                                                       >  2:534149/1‑92 (MQ=255)
                                                           tctcAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAAct                                                                                                       <  1:534149/92‑1 (MQ=255)
                                                            ctcAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCa                                                       >  2:1273942/1‑139 (MQ=255)
                                                               aaTGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCAt                                                      <  1:792396/137‑1 (MQ=255)
                                                                aTGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAg                                                   <  1:130535/139‑1 (MQ=255)
                                                                     tAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATc                                               <  1:1036305/138‑1 (MQ=255)
                                                                     tAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATc                                               >  2:1036305/1‑138 (MQ=255)
                                                                        ttattaTTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTa                                           >  2:1216268/1‑139 (MQ=255)
                                                                             attGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAgg                                                                         <  2:1127145/104‑1 (MQ=255)
                                                                             attGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAgg                                                                         >  1:1127145/1‑104 (MQ=255)
                                                                             attGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTa                                      >  2:1255705/1‑139 (MQ=255)
                                                                             attGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTa                                      >  2:101985/1‑139 (MQ=255)
                                                                                 ccATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataat                                                              <  2:64465/111‑1 (MQ=255)
                                                                                 ccATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGataataat                                                              >  1:64465/1‑111 (MQ=255)
                                                                                 ccATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTAtgtg                                  >  2:793463/1‑139 (MQ=255)
                                                                                    tcatTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAg                                                                                                                   <  1:605719/55‑1 (MQ=255)
                                                                                    tcatTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAg                                                                                                                   >  2:605719/1‑55 (MQ=255)
                                                                                    tcatTCGTGACGTTATTACCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACACTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCg                               <  2:1158736/139‑1 (MQ=255)
                                                                                                  attaCCAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCAtt                 >  1:992338/1‑139 (MQ=255)
                                                                                                      ccAGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGAc             <  1:843394/139‑1 (MQ=255)
                                                                                                        aGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGt           <  1:101985/139‑1 (MQ=255)
                                                                                                        aGTAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGt           <  1:1255705/139‑1 (MQ=255)
                                                                                                          tAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGAc             >  1:131449/1‑135 (MQ=255)
                                                                                                          tAAGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGAc             <  2:131449/135‑1 (MQ=255)
                                                                                                            aGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTa                                                 <  1:1141733/97‑1 (MQ=255)
                                                                                                            aGTGTGATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTa                                                 >  2:1141733/1‑97 (MQ=255)
                                                                                                                tgATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTg     >  2:1061013/1‑137 (MQ=255)
                                                                                                                 gATGTTGCAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGTAg  >  2:111916/1‑139 (MQ=255)
                                                                                                                       |                                                                                                                                    
AAATCAGCCTGAGATTCTGTCGTTTTATTCGCGATCTTATCCAGACCCATCTTAGGGATTCTCAATGTGTAGTTATTATTGCCATCATTCGTGACGTTATTACCAGTAAGTGTGATGTTACAGGTTGTCGCTTTGACAGTGGCTGTAAACTCTACATTAATATCCGTTGCAAATACTGAGGGTGATAATAATATTCCCATCAGTATCACTATTTTATGTGGCGTTGTTTTTATCATTTGACGTCCCTTGTAG  >  NZ_CP009273/149064‑149315

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