Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,452,315 C→T 16.3% A8T (GCC→ACC)  yfcS ← putative fimbrial chaperone YfcS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,452,3150CT16.3% 52.6 / 3.4 28A8T (GCC→ACC) yfcSputative fimbrial chaperone YfcS
Reads supporting (aligned to +/- strand):  ref base C (13/9);  new base T (2/2);  total (16/12)
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.82e-01

GATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCATTAAATATCAAACGGGTACGATCTGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGTCATCAGGAACAACCGG  >  NC_000913/2452174‑2452451
                                                                                                                                             |                                                                                                                                        
gATCGTTTTTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCATTAAATATCAAACGGGTACGATCTGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCa                                                                                                                                         <  1:298712/143‑1 (MQ=255)
 aTCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCATTAAATATCAAACGGGTACGATCTGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCag                                                                                                                                        >  1:56919/1‑143 (MQ=255)
                       cTGATTGATTTGTCGCTTTCATTAAATATCAAACGGGTACGATCTGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGa                                                                                                                  <  2:136506/143‑1 (MQ=255)
                             gattTGTCGCTTTCATTAAATATCAAACGGGTACGATCTGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTc                                                                                                            <  1:202086/143‑1 (MQ=255)
                               tttGTCGCTTTCATTAAATATCAAACGGGTACGATCTGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCt                                                                                                          <  2:244927/143‑1 (MQ=255)
                                             tAAATATCAAACGGGTACGATCTGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAAACAGACATGGATCATTCCTTAATCACGTTATTg                                                                                            >  1:223574/1‑143 (MQ=255)
                                                  atCAAACGTGTACGATCTGGCGTTACTGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAg                                                                                       <  1:19984/143‑1 (MQ=255)
                                                  atCAAACGGGGACGATCTGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAg                                                                                       >  2:231680/1‑143 (MQ=255)
                                                     aaaCGGGTACGATCTGGCGTTACAGACGCCAGGGCAGACAAACTAGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGAAGAGCAAAGTAAATCAGACATGGAACATTCCTAAATAACGTGATTGGTAAGcag                                                                                    >  2:18821/1‑141 (MQ=255)
                                                                         tACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTGAATCACGTTATTGGTTAGCTGAAACAGAAAAAGGGGGATCg                                                                >  2:163979/1‑143 (MQ=255)
                                                                           cAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGTCAGCGTCATAGCACCCAATTTTGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCggt                                                              <  2:117159/143‑1 (MQ=255)
                                                                                     ggCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGCGGACGGGTGGGTGTCAcc                                                    >  2:301941/1‑143 (MQ=255)
                                                                                         gACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGTAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGtt                                                >  1:257388/1‑143 (MQ=255)
                                                                                         gACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGTAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGtt                                                >  1:100161/1‑143 (MQ=255)
                                                                                                tGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGCTCATTCCTCAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGGCCGCTGTTACGGCTg                                         >  2:129869/1‑143 (MQ=255)
                                                                                                 cttAGCGCTAAGCATTTACGCCAGCGTCTTCTGACTCATTTTGGCAGATCAACGTAAATCAGACATGTATCATTCCTTACTCACGTTATTGGTTAGCTGAAACAGAAACTGGGGGTTCTTTGGGTGTCACCTGTTGCGTCTGc                                        <  2:43072/141‑1 (MQ=255)
                                                                                                  gTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGCAACAGGGGGATCGGTGGGTGGCACCTGCTGCCTCTGca                                       >  1:228507/1‑142 (MQ=255)
                                                                                                     gcgcTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCt                                    >  1:210307/1‑143 (MQ=255)
                                                                                                      cgcTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACCGGGGGATCGGTGGGTGTCCCCTGTTGCGTCTGCTCCTg                                   >  1:18040/1‑143 (MQ=255)
                                                                                                         tAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATCGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGCATCGGTGGGTGTCACCTGCTGCGTCTCCTCCTGCAt                                >  2:62984/1‑143 (MQ=255)
                                                                                                           aGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGa                              >  1:244200/1‑143 (MQ=255)
                                                                                                                   cGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCgtt                      >  2:191320/1‑143 (MQ=255)
                                                                                                                     atagtgTCATAGCACGCAATGTGGGAGAGCAAGGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTAt                    <  1:282043/138‑1 (MQ=255)
                                                                                                                                gCACCCAATTTGGTAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGTCATCAGGaa         <  2:100161/143‑1 (MQ=255)
                                                                                                                                gCACCCAATTTGGTAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGTCATCAGGaa         <  2:257388/143‑1 (MQ=255)
                                                                                                                                    ccAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACACGGGGATCTGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGTCATCAGGaacaac     <  2:150543/143‑1 (MQ=255)
                                                                                                                                    ccAAGTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTGATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGTCATCAGGaacaac     <  2:56919/143‑1 (MQ=255)
                                                                                                                                       aTTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTCCAACAGACATCGGCGGTTCTGTGGGTGTCACCTGTTGCATCTACAAATGAATAAAAGAAATTATGTCATAAGGAACAAACgg  >  1:111679/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                        
GATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCATTAAATATCAAACGGGTACGATCTGGCGTTACAGACGCCAGGGCAGACAAACTGGTAGCGCTAAGCAGTAACGCCAGCGTCATAGCACCCAATTTGGCAGAGCAAAGTAAATCAGACATGGATCATTCCTTAATCACGTTATTGGTTAGCTGAAACAGAAACAGGGGGATCGGTGGGTGTCACCTGTTGCGTCTGCTCCTGAATGACAGGCTTTATGTCATCAGGAACAACCGG  >  NC_000913/2452174‑2452451

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