Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,288,032 C→A 8.5% intergenic (+22/‑58) yraH → / → yraI putative fimbrial protein YraH/putative fimbrial chaperone YraI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,288,0320CA 8.5% 116.4 / 2.9 47intergenic (+22/‑58)yraH/yraIputative fimbrial protein YraH/putative fimbrial chaperone YraI
Reads supporting (aligned to +/- strand):  ref base C (16/27);  new base A (2/2);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTA  >  NC_000913/3287891‑3288173
                                                                                                                                              |                                                                                                                                             
gAATCCAATGCTATCGTCAATGTTGACCTCATCGCTTACTACGTATCCACCGTAGTTCCATCAGAAACCACCCCCGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCt                                                                                                                                              <  2:499499/143‑1 (MQ=255)
  aTCCAATGGTATTGTCAAGGTTGACCTCACCGCTTACTACGGATCAACCGTAGTCCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGt                                                                                                                                            <  1:66550/143‑1 (MQ=255)
     cAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAg                                                                                                                                         >  2:368345/1‑143 (MQ=255)
      aaTGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGa                                                                                                                                        >  1:28746/1‑143 (MQ=255)
        ttctagtGTCTATGTTGATTTCATCGCTTATTACGTATTAACCGTAGTTCCATCAGAGATCACCCCTGTGTAAGCTGATGCTTTCGTTATCGTATCTCTTGTTTTCCTTTAATACGTTAGGGCGTTATCTGACCTGTCAGATa                                                                                                                                      <  2:46958/137‑1 (MQ=255)
        tGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAACTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGCTAGGGCGTTATCTGACCCGCCAgtga                                                                                                                                      >  2:108565/1‑140 (MQ=255)
         gttATTGTGCAGGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCATAAATCAACCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATaa                                                                                                                                     <  2:455535/141‑1 (MQ=255)
         ggTATTGTCAATGTTGACCTCTTCGCTTACTAGTTATCAACCGTAGTTCTATCATAAATCACCCCAGGGGCAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATaa                                                                                                                                     <  2:455152/143‑1 (MQ=255)
         ggTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                     <  2:49191/143‑1 (MQ=255)
              tGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGccc                                                                                                                                >  1:584768/1‑143 (MQ=255)
               gTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCt                                                                                                                               >  1:256303/1‑143 (MQ=255)
                 cactGTTGACCTCATCGCTTACTACGTATTAACCGTCGTTCCATCAGACGAATCGTTGATGAGATCTTTCGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCttt                                                                                                                             <  1:307464/140‑1 (MQ=255)
                  taTGTTGACCGTATCGCTTACTACGTATCTACCGTAGTTCCATCAGAAATCAGCCGTTGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCtttt                                                                                                                            <  2:240159/142‑1 (MQ=255)
                    tGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcc                                                                                                                          <  2:202712/143‑1 (MQ=255)
                        gACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCGCCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttccttcc                                                                                                                      <  2:596030/143‑1 (MQ=255)
                          ccTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctc                                                                                                                    <  2:220619/143‑1 (MQ=255)
                                  cTTACTACGTATCAACCGTAGTTCAATCAGAGATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGt                                                                                                            <  1:546108/143‑1 (MQ=255)
                                      ctacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTa                                                                                                        >  1:457122/1‑143 (MQ=255)
                                       tgcGTAGCAACCGTAGTGCCATCAGAAATCACCCCAGGCGAAGCAGACGTTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTAt                                                                                                       <  2:41929/141‑1 (MQ=255)
                                             tCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGtt                                                                                                 >  2:345459/1‑143 (MQ=255)
                                               aaCCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGCTTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTGCTTCATCTTTCTCGTTGTATCAGGTTGa                                                                                               >  2:153678/1‑143 (MQ=255)
                                                  cGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGATGTATCAGGATGaaaa                                                                                            >  2:13487/1‑143 (MQ=255)
                                                    tAGCTCCATCAGAAATCACCCCAGTGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTTTTTCTCGTTGTATCAGGTTGAAAAAt                                                                                          <  1:445697/143‑1 (MQ=255)
                                                      gTTCCATCAGCTACCGCCCCCGCGAAAGCAGACGCTATCGTTAACGTAACGCTGGATTATCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATaa                                                                                        <  2:332514/143‑1 (MQ=255)
                                                          caccAGACAGCACCCCAGGGGAAGCAGACGCCAACGTTAACGTAACGCTCGATTACCGTTAATACGATCGGGCGTTATCTGACCTGTCAGATAACGCCATTAGCTATTCTCTTTCTCGTTGTATCATGTTTAAAAATAACTAt                                                                                    <  2:113308/140‑1 (MQ=255)
                                                           aTCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCCCTTTCTCGTTGTCTCAGGTTGAGACGTAACTATg                                                                                   >  1:434542/1‑143 (MQ=255)
                                                            tCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGt                                                                                  >  2:545469/1‑143 (MQ=255)
                                                                             ggAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTcgcg                                                                 <  1:191946/143‑1 (MQ=255)
                                                                               gaaaccaaCGCTATCGTCAACGTAACACAGGCTTACCGTTAAAACGCTAGGGCCACGACGGAACTGTACATTGACGCCCTTTTCCTTCCTCTTTATAATAGTGTAAGGTGGATAAATAACTATGTCAAAACGAACATTCGCgg                                                                <  2:169935/136‑1 (MQ=255)
                                                                                 tgcataggCTCTCGTTCGCGTGACGCTGGATTTGCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATTACGTCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTg                                                              <  2:19923/136‑1 (MQ=255)
                                                                                      aagttgtCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATAtt                                                         <  2:256303/137‑1 (MQ=255)
                                                                                         cTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCGGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAc                                                      <  1:190066/143‑1 (MQ=255)
                                                                                                   cgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTGGCAAAAGAACTATGTCAAAACGAACATTCGCGGTGAGATTAAACTAGTTttgt                                            >  2:538227/1‑143 (MQ=255)
                                                                                                        cgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCtt                                       <  2:584768/143‑1 (MQ=255)
                                                                                                         gCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGACAATAATCGATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTc                                      >  1:150945/1‑143 (MQ=255)
                                                                                                               ttAGCGGTAGTAGGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTTGTTGTACCACATCCAGATTTTCCCGGGTCAATACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTAtt                                <  2:477585/143‑1 (MQ=255)
                                                                                                               ttACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTAtt                                <  1:368345/143‑1 (MQ=255)
                                                                                                                tgccGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTGTCCTTCCTCTGTCTCGTTGTATCAGGTTGAAAAATAACTAGGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTg                               <  2:127488/141‑1 (MQ=255)
                                                                                                                tACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTg                               >  2:507638/1‑143 (MQ=255)
                                                                                                                   cGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGcc                            <  1:163743/143‑1 (MQ=255)
                                                                                                                     ttAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTGTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAg                          <  2:385008/143‑1 (MQ=255)
                                                                                                                      tAATGCTTTAGGGCGTTATCTGACCAGTCAGATAGCGCCCTGGGCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAgg                         <  1:231504/143‑1 (MQ=255)
                                                                                                                                   cGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTGTGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCaggggg            >  1:498569/1‑143 (MQ=255)
                                                                                                                                   cGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGGGCTTGCaggagg            >  2:321775/1‑143 (MQ=255)
                                                                                                                                       aTCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGACAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCTCTTGCCGGAGGAATc        >  2:55586/1‑143 (MQ=255)
                                                                                                                                          tGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGtt     >  1:372728/1‑143 (MQ=255)
                                                                                                                                           gACCTGTCAGATAACGCCCGTTTCCTTCCTCTTGCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt    <  1:108565/143‑1 (MQ=255)
                                                                                                                                             ccTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTAAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTa  >  2:441505/1‑143 (MQ=255)
                                                                                                                                             ccTGGCAGATAACGCCCTCTTCCTTCCACTGTCTGGTGGTATCAGGTTGAAAAATAGCTATGTCAAACCGAACCTTCGCGGTGAGATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTa  <  2:245726/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTA  >  NC_000913/3287891‑3288173

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: