Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 563,871 C→T 10.5% A204T (GCC→ACC)  fimZ ← putative LuxR family transcriptional regulator FimZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913563,8710CT10.5% 93.7 / 3.9 38A204T (GCC→ACC) fimZputative LuxR family transcriptional regulator FimZ
Reads supporting (aligned to +/- strand):  ref base C (20/14);  new base T (2/2);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-02

GATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATA  >  NC_000913/563732‑564011
                                                                                                                                           |                                                                                                                                            
gATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTa                                                                                                                                           >  1:76685/1‑143 (MQ=255)
  tGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGt                                                                                                                                         >  1:54263/1‑143 (MQ=255)
                  acacTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCCCTACAAt                                                                                                                         >  2:559445/1‑143 (MQ=255)
                      tAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGaa                                                                                                                     <  1:632646/143‑1 (MQ=255)
                      tAAACTTTTTTCCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGCTCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGaa                                                                                                                     <  1:991676/143‑1 (MQ=255)
                           ttttttGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGCTAAGCTCTACAATTGAATGaaa                                                                                                                >  1:324090/1‑143 (MQ=255)
                                    aaTTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGttt                                                                                                       >  1:483992/1‑143 (MQ=255)
                                    aaTTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCtt                                                                                                       >  2:62814/1‑143 (MQ=255)
                                            aGCGTTAGCCAGAGGTGGACGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATaaa                                                                                               <  1:335084/143‑1 (MQ=255)
                                               gTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACATAGCTTGCCATAAatat                                                                                            >  2:460119/1‑143 (MQ=255)
                                                 tAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTa                                                                                          <  1:724128/143‑1 (MQ=255)
                                                  aGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTCTATTCATTCGTATAATTTGGCGTAGTCGATAATCTCTACAATTGAATGCAAACCTAGCTTACCATACATATTAg                                                                                         >  2:480260/1‑143 (MQ=255)
                                                  aGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTATACAATTGAATGCAAACCCAGATTGCCATAAATATTAg                                                                                         >  2:483885/1‑143 (MQ=255)
                                                   gCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAGGCTTTACAATTGATTGCAACCGTAGCTTGTCGTAAATATTCGa                                                                                        >  2:591808/1‑143 (MQ=255)
                                                      aGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGAttt                                                                                     <  1:366373/143‑1 (MQ=255)
                                                        cGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTAt                                                                                   >  1:622260/1‑143 (MQ=255)
                                                                 gCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATACTTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTaa                                                                          >  2:574506/1‑143 (MQ=255)
                                                                          ggCCAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTAt                                                                 <  1:163117/143‑1 (MQ=255)
                                                                           gCGAAAGCGGATGCGATATTTATGATCAATTATACATAGGAGAATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTAtt                                                                <  2:446900/143‑1 (MQ=255)
                                                                           gCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTCtt                                                                >  1:823319/1‑143 (MQ=255)
                                                                             gAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGTGTAGTCGATAAGCTCTACaa                                                                                                                          <  1:338217/83‑1 (MQ=255)
                                                                             gAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGTGTAGTCGATAAGCTCTACaa                                                                                                                          <  1:338390/83‑1 (MQ=255)
                                                                             gAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGTGTAGTCGATAAGCTCTACaa                                                                                                                          >  2:338217/1‑83 (MQ=255)
                                                                             gAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGTGTAGTCGATAAGCTCTACaa                                                                                                                          >  2:338390/1‑83 (MQ=255)
                                                                             gAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAATCTCTACAATTGAATGCAACCCTAGCTTGCCATACATATTAGATTTATGCGCACTAACTGTTTTATTGc                                                              >  1:533880/1‑143 (MQ=255)
                                                                                aGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTaa                                                           <  1:253516/143‑1 (MQ=255)
                                                                                    gATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGACATAAATATTAGATTTATGCGCAATAACTGTTTTATTGCTAAGtaa                                                       >  2:211765/1‑143 (MQ=255)
                                                                                       gCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAACCCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGtaataa                                                    <  2:240080/143‑1 (MQ=255)
                                                                                                   aTCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGTAATAACTTATCGGCAAt                                        >  1:803969/1‑143 (MQ=255)
                                                                                                   aTCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGTAATAACTTATCGGCAAt                                        >  1:6325/1‑143 (MQ=255)
                                                                                                           taattaatATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGTAATAACTTATCGGCAATTTCTttat                                <  2:315346/143‑1 (MQ=255)
                                                                                                           taagtaatATTATATTAATTCGTATAATTTGGCGTAGTCGAAAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGTAGTAACTTATCGGCAATTTCttttt                                >  2:886378/1‑141 (MQ=255)
                                                                                                                 acattatattaATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATa                          <  2:685878/141‑1 (MQ=255)
                                                                                                                  tattatattaATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTACTAAGTAATAACTTATCGGCAATTTCTTTATTAGATaa                         >  2:855984/1‑143 (MQ=255)
                                                                                                                               cGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGa            <  2:924370/143‑1 (MQ=255)
                                                                                                                                 tataATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTCTTGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATa          >  2:170071/1‑143 (MQ=255)
                                                                                                                                        tGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAat   <  2:54263/143‑1 (MQ=255)
                                                                                                                                         ggCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAata  <  1:42264/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTAACTGTTTTATTGCTAAGTAATAACTTATCGGCAATTTCTTTATTAGATAATCCGCTAACCAGATAACGTAATA  >  NC_000913/563732‑564011

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