Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,448,105 G→A pseudogene (376/2646 nt) yfcU ← fimbrial usher protein YfcU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,448,1050GA96.8% 231.4 / ‑2.1 63pseudogene (376/2646 nt)yfcUfimbrial usher protein YfcU
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (31/30);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

CGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATACCAGTAA  >  NZ_CP009273/2447970‑2448238
                                                                                                                                       |                                                                                                                                     
cGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACtt                                                                                                                                    <  2:683799/139‑1 (MQ=255)
  gagaTGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACtt                                                                                                                                    >  1:286931/1‑137 (MQ=255)
        ccATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCAtt                                                                                                                            >  1:1186324/1‑139 (MQ=255)
            catcCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTc                                                                                                                                   >  1:18267/1‑128 (MQ=255)
            catcCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTc                                                                                                                                   <  2:18267/128‑1 (MQ=255)
                 cAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGt                                                                                                                   <  2:1155812/139‑1 (MQ=255)
                                      tCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCg                                                                                              >  1:66835/1‑139 (MQ=255)
                                              ggTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGtct                                                                                      <  2:138298/139‑1 (MQ=255)
                                                  tatTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTcttctt                                                                                  <  2:447590/139‑1 (MQ=255)
                                                    ttCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGcc                                                                                            >  1:421429/1‑127 (MQ=255)
                                                    ttCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGcc                                                                                            <  2:421429/127‑1 (MQ=255)
                                                    ttCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCtttt                                                                                <  2:201675/139‑1 (MQ=255)
                                                        aGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAAc                                                                            >  1:1055751/1‑139 (MQ=255)
                                                         ggTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAAc                                                                            >  2:1204595/1‑138 (MQ=255)
                                                         ggTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAAc                                                                            >  1:1048066/1‑138 (MQ=255)
                                                         ggTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAAc                                                                            >  2:228585/1‑138 (MQ=255)
                                                              gCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGCTTTTTCGCCACGTCTTCTTTTAAACCAAACt                                                                      >  2:183667/1‑139 (MQ=255)
                                                               ccTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTg                                                                     >  1:153840/1‑139 (MQ=255)
                                                                   cGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTgcgc                                                                  >  1:753758/1‑138 (MQ=255)
                                                                   cGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTacgc                                                                  >  2:209818/1‑138 (MQ=255)
                                                                         tGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGt                                                           >  1:159677/1‑139 (MQ=255)
                                                                           aaaTGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTc                                                         <  2:614863/139‑1 (MQ=255)
                                                                                aCTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGtc                                                                                       <  2:52585/104‑1 (MQ=255)
                                                                                aCTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGtc                                                                                       >  1:52585/1‑104 (MQ=255)
                                                                                aCTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGgtg                                                    >  1:1277378/1‑139 (MQ=255)
                                                                                aCTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGgtg                                                    <  2:840705/139‑1 (MQ=255)
                                                                                aCTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGgt                                                     <  2:947946/138‑1 (MQ=255)
                                                                                   aaTGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGgtgt                                                   >  2:5061/1‑137 (MQ=255)
                                                                                    aTGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAg                                                <  2:961163/139‑1 (MQ=255)
                                                                                         gACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAg                                           >  2:98940/1‑139 (MQ=255)
                                                                                         gACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAg                                           >  2:814271/1‑139 (MQ=255)
                                                                                         gACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAg                                           <  1:115163/139‑1 (MQ=255)
                                                                                          aCTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTcttctt                                                                                  >  1:129023/1‑99 (MQ=255)
                                                                                          aCTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTcttctt                                                                                  <  2:129023/99‑1 (MQ=255)
                                                                                            tGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTgcgc                                                                  <  1:4476/113‑1 (MQ=255)
                                                                                            tGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTgcgc                                                                  >  2:4476/1‑113 (MQ=255)
                                                                                                   aaTCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCtttt                                 >  2:743904/1‑139 (MQ=255)
                                                                                                        gCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTAcc                                                                                                                        <  2:887972/47‑1 (MQ=255)
                                                                                                        gCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTAcc                                                                                                                        >  1:887972/1‑47 (MQ=255)
                                                                                                        gCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTaa                            <  1:98940/139‑1 (MQ=255)
                                                                                                        gCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTaa                            >  2:1259468/1‑139 (MQ=255)
                                                                                                        gCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTaa                            >  1:136218/1‑139 (MQ=255)
                                                                                                        gCCTTAATTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTaa                            >  2:488529/1‑139 (MQ=255)
                                                                                                             aaTTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTc                        <  2:1055751/138‑1 (MQ=255)
                                                                                                             aaTTTCCACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTc                        <  1:1259468/138‑1 (MQ=255)
                                                                                                                   cACGCCTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCa                 >  1:781112/1‑139 (MQ=255)
                                                                                                                        cTTCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAg             >  1:1029420/1‑138 (MQ=255)
                                                                                                                          tCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGccac                                                               >  1:541175/1‑86 (MQ=255)
                                                                                                                          tCCAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGccac                                                               <  2:541175/86‑1 (MQ=255)
                                                                                                                           ccAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCa           <  1:358646/137‑1 (MQ=255)
                                                                                                                            ccgTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAGCCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATa         <  1:548759/136‑1 (MQ=255)
                                                                                                                            cAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATa         <  2:136218/138‑1 (MQ=255)
                                                                                                                            cAGTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATa         <  1:1243099/138‑1 (MQ=255)
                                                                                                                              gTTGACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATACCa      <  2:1228266/139‑1 (MQ=255)
                                                                                                                                 gACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATACCAGt    >  2:1246580/1‑138 (MQ=255)
                                                                                                                                 gACCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATACCAGt    <  2:1186324/138‑1 (MQ=255)
                                                                                                                                  aCCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATACCAGt    >  2:1013247/1‑137 (MQ=255)
                                                                                                                                  aCCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATACCAGt    <  1:1013247/137‑1 (MQ=255)
                                                                                                                                  aCCGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATACCAGTaa  <  1:684956/139‑1 (MQ=255)
                                                                                                                                    cGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGcc                                                                                            >  2:845970/1‑47 (MQ=255)
                                                                                                                                    cGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGcc                                                                                            <  1:845970/47‑1 (MQ=255)
                                                                                                                                    cGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATACCAGt    <  2:598927/135‑1 (MQ=255)
                                                                                                                                    cGGACTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATACCAGt    >  1:598927/1‑135 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CGGAGATGCCATCATCCCAACGAGAAGGTGGATCCCAGTCGGGCCAGGTATATTCGAGGTAAGCCTGCGGTAGTGAAATGACTAATGCGGACTGGCTTAAATCAGCCTTAATTTCCACGCCTTCCAGTTGACCGGGCTTCAGGCATTTACCATCGTGGCTCCATTGCAGATTTTTCGCCACGTCTTCTTTTAAACCAAACTGCGCCACCAGTTCCGGTGTCAGACAAGCATAGCTTTTGCTAACGTCATCTTCACCAGCATACCAGTAA  >  NZ_CP009273/2447970‑2448238

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