Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,356,378 (A)8→7 coding (213/2382 nt) yhcD → fimbrial biogenesis outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,356,3710A.96.5% 221.9 / 2.1 58coding (206/2382 nt)yhcDfimbrial biogenesis outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base . (46/10);  total (48/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

GTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAA  >  NZ_CP009273/3356284‑3356455
                                                                                       |                                                                                    
gtgtAGATCCTAATATATCGCATTATTTTTCCCGTTCAACCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaa                                                                                    >  2:461924/1‑90 (MQ=255)
   tAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaa                                                                                 <  2:318370/90‑1 (MQ=255)
    aGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaa                                                                                >  2:61586/1‑90 (MQ=255)
         cTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCaa                                                                          >  2:79598/1‑90 (MQ=255)
          tAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAc                                                                         >  2:915292/1‑90 (MQ=255)
          tAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAc                                                                         >  1:459539/1‑90 (MQ=255)
          tAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAc                                                                         >  1:532268/1‑90 (MQ=255)
            atatatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACAt                                                                       <  2:829429/90‑1 (MQ=255)
            atatatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACAt                                                                       >  1:737554/1‑90 (MQ=255)
               tatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAc                                                                         <  1:386286/85‑1 (MQ=255)
               tatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAc                                                                         >  2:386286/1‑85 (MQ=255)
               tatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGc                                                                    >  2:41731/1‑90 (MQ=255)
               tatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGc                                                                    >  1:320980/1‑90 (MQ=255)
               tatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGc                                                                    >  2:15394/1‑90 (MQ=255)
                atCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCt                                                                   >  1:252375/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:232007/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:297905/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:262485/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:917398/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  1:830070/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  <  1:517775/90‑1 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:559225/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:770527/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  1:388672/1‑90 (MQ=255)
                   gCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTAcg                                                                >  1:746805/1‑90 (MQ=255)
                   gCATTATTTTTCCCGTCCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTAcg                                                                >  2:524255/1‑90 (MQ=255)
                        attTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCtt                                                           <  2:513683/90‑1 (MQ=255)
                        attTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCtt                                                           >  2:692403/1‑90 (MQ=255)
                        attTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCtt                                                           <  1:232865/90‑1 (MQ=255)
                             tCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTtgatg                                                      >  2:437212/1‑90 (MQ=255)
                                cGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGaaa                                                   >  1:682051/1‑90 (MQ=255)
                                  ttCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAAt                                                 >  2:112546/1‑90 (MQ=255)
                                      gCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCg                                             >  2:81728/1‑90 (MQ=255)
                                           ttttttGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                        >  1:431691/1‑90 (MQ=255)
                                           ttttttGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                        <  2:362845/90‑1 (MQ=255)
                                           ttttttGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                        <  1:614410/90‑1 (MQ=255)
                                           ttttttGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                        <  2:508286/90‑1 (MQ=255)
                                                    aGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGAt                               >  1:773432/1‑90 (MQ=255)
                                                    aGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGAt                               >  1:45114/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTc                     >  2:620000/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  1:594855/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  1:239391/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  2:8807/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  1:388260/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  1:436270/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  1:856747/1‑90 (MQ=255)
                                                                  tGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGca                 >  2:738893/1‑90 (MQ=255)
                                                                    aTAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGcaac               >  1:680691/1‑90 (MQ=255)
                                                                       gTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAg            >  2:399675/1‑90 (MQ=255)
                                                                       gTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAg            >  1:366880/1‑90 (MQ=255)
                                                                        tATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTc                                    <  2:424188/65‑1 (MQ=255)
                                                                        tATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTc                                    >  1:424188/1‑65 (MQ=255)
                                                                         aTCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGcc          <  2:12229/90‑1 (MQ=255)
                                                                           cAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCg         >  1:724258/1‑90 (MQ=255)
                                                                           cAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCgg        >  1:565378/1‑90 (MQ=255)
                                                                             gTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGtt      >  2:189686/1‑90 (MQ=255)
                                                                                 aTGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa  >  1:504158/1‑90 (MQ=255)
                                                                                 aTGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa  >  1:674431/1‑90 (MQ=255)
                                                                                 aTGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa  >  1:717026/1‑90 (MQ=255)
                                                                                 aTGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa  >  2:180715/1‑90 (MQ=255)
                                                                                 aTGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa  >  1:157084/1‑90 (MQ=255)
                                                                                       |                                                                                    
GTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAA  >  NZ_CP009273/3356284‑3356455

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