Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 249,530 (A)5→4 pseudogene (541/1713 nt) lfhA ← flagellar system protein, promoterless fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913249,5260A.97.7% 173.1 / ‑2.1 44pseudogene (545/1713 nt)lfhAflagellar system protein, promoterless fragment
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (19/24);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACG  >  NC_000913/249391‑249657
                                                                                                                                       |                                                                                                                                   
tCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCGAGCCGGTGAAAACAGGCAGGGAGCTGAACAGCaaaac                                                                                                                                 >  1:226773/1‑139 (MQ=255)
tCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCACCAACCGGTAAAACCACGCAGAGCGCTGAACAGCaaaac                                                                                                                                 >  2:285185/1‑139 (MQ=255)
      cTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAgg                                                                                                                          >  1:443655/1‑140 (MQ=255)
       tGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAggg                                                                                                                         >  2:827052/1‑139 (MQ=255)
       tGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGt                                                                                                                         >  1:766765/1‑140 (MQ=255)
         aCCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGc                                                                                                                       >  1:279762/1‑140 (MQ=255)
         aCCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCCCCAGCCGGGAAAACCAAGCAGGGCGCTGAATAGCAAAACGGCAGGGGc                                                                                                                       >  2:287264/1‑140 (MQ=255)
            tgctgcTCGCTGGTTTCGGTGAAGGGGCGTGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGc                                                                                                                    <  1:287264/140‑1 (MQ=255)
               tgctCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCAtt                                                                                                                 <  1:662258/140‑1 (MQ=255)
                 ctCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGTAGGGCGCAGAAAAGAAAAACGGCAGGTGCGGCATTcc                                                                                                               >  1:331951/1‑140 (MQ=255)
                   cGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGGGCGGCATTCCCg                                                                                                             >  2:674535/1‑140 (MQ=255)
                        gTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTacc                                                                                                        >  1:766904/1‑140 (MQ=255)
                              gTGATGTTGGGGGTCACAGTTTCGAGGCTTTTCCCATCCTCCGCCGCCTCCAGCGGCTGTTTCCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGcc                                                                                                  <  2:226773/140‑1 (MQ=255)
                              gTGATGGTGCGGGTCAGGGTTTGGAGGCTTTGCTCGTCCTCCTCCGCCCCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGcc                                                                                                  <  2:302953/140‑1 (MQ=255)
                              gTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGcc                                                                                                  <  2:766765/140‑1 (MQ=255)
                                  tGGTGCGGGTCAGCGTTTCGAGGTTTTTCTCTCCCGCCTCGGCCGCCGGGGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCa                                                                                              <  2:660394/140‑1 (MQ=255)
                                  tGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCGTTCCCGGGACCACCGCCAGcc                                                                                              >  2:116536/1‑139 (MQ=255)
                                  tGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCa                                                                                              >  1:54418/1‑140 (MQ=255)
                                  tGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCa                                                                                              <  1:674535/140‑1 (MQ=255)
                                  tGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCa                                                                                              <  1:742061/140‑1 (MQ=255)
                                   ggTGGGGGTCAGCGTTTCGGGGTTTTTCTCTCCCGCCTCCGCCGCCAGGGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCAc                                                                                             <  2:443655/140‑1 (MQ=255)
                                    gTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAACCGGCAGGCGCGCCAGTCCCGGTACCACCGCCAGCACg                                                                                            >  1:296819/1‑140 (MQ=255)
                                      gCGGGTCAGCGTTTCGGGGCTTGTCTCTCCCCCCCCCGCCGCCGGCGGTTGTTTGCTCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGaa                                                                                          <  2:306583/140‑1 (MQ=255)
                                         ggTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCCGGTGCGGGATTCCCGGTACCACCGCCAGCACGAACAt                                                                                       >  1:260879/1‑140 (MQ=255)
                                            cAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAAt                                                                                    <  2:752272/140‑1 (MQ=255)
                                             aGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGTGGCGGCATTCCCGGGAACACCACCAGCACGAACATAATc                                                                                   >  2:239580/1‑140 (MQ=255)
                                             aGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTACCGGTACCACCGCTAGCACGAACATAATc                                                                                   >  1:267317/1‑140 (MQ=255)
                                               cGTTTCGAGGCTTTTGTCTTCCACCTCCGACGCCAGCTGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATccc                                                                                 <  1:285185/140‑1 (MQ=255)
                                                    ggAGGCTTTTCTCTCCCCCCCCCGCGGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTgg                                                                            <  1:777518/139‑1 (MQ=255)
                                                    cggggCTTTTCTCTTCCGCCCCTGCCGCCAGGGGCGGTTTGCTCATCCCCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTgg                                                                            <  2:279762/137‑1 (MQ=255)
                                                    cGAGGTTTTTCTCTCCCGCCCCCGCCGCCGGCGGCTTTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTgg                                                                            <  2:766904/140‑1 (MQ=255)
                                                    cGAGGCTTTTCTTTTCCCCCTCCCCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTgg                                                                            <  2:568562/140‑1 (MQ=255)
                                                    cGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCCGGATTCCCGGTACCCCCTCCAGCACGAACATAATccccgggg                                                                            >  2:783762/1‑135 (MQ=255)
                                                      ggctttttCTCTTCCGCCCCCGGCGCCAGCGCCTGATTGATCACCCGCGAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCg                                                                          <  2:331951/136‑1 (MQ=255)
                                                         cTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGgtg                                                                       >  1:783627/1‑140 (MQ=255)
                                                              ctctTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTagag                                                                  <  1:827052/140‑1 (MQ=255)
                                                                      cctccgccGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACg                                                          >  1:484212/1‑140 (MQ=255)
                                                                      cctccgccGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACg                                                          >  1:494487/1‑140 (MQ=255)
                                                                                    ggCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCagcagc                                            <  2:783627/140‑1 (MQ=255)
                                                                                            gCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACg                                    >  1:118521/1‑140 (MQ=255)
                                                                                              tCATCCGCCAGCCGGTAAAACCTAGCAGGGCCCTGGACATCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCa                                  <  2:286609/140‑1 (MQ=255)
                                                                                                   cGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCa                             <  2:494487/140‑1 (MQ=255)
                                                                                                   cGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCa                             <  2:484212/140‑1 (MQ=255)
                                                                                                       aGCCGGTAAAACCAACCAGGGCGCTCAACTGCAAAACAGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATggg                         <  1:239580/140‑1 (MQ=255)
                                                                                                                          ggCGCTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCt      >  1:457283/1‑140 (MQ=255)
                                                                                                                              cTGAACAGCAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGAcg  <  1:616128/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
TCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACG  >  NC_000913/249391‑249657

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