Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,383,956 A→T 8.9% intergenic (+32/+39) elaD → / ← yfbK protease ElaD/IPR002035/DUF3520 domain‑containing protein YfbK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,383,9560AT 8.9% 124.9 / 5.1 46intergenic (+32/+39)elaD/yfbKprotease ElaD/IPR002035/DUF3520 domain‑containing protein YfbK
Reads supporting (aligned to +/- strand):  ref base A (15/25);  new base T (2/2);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

ATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAC  >  NC_000913/2383815‑2384094
                                                                                                                                             |                                                                                                                                          
aTTGTTACTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAACACTATATCAGCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCaaa                                                                                                                                           <  1:664297/143‑1 (MQ=255)
    tcaTTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAACTAAAGCACTCACAACCCAAAAGAAGTTCTGCCCGTAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAatagg                                                                                                                                       <  1:197663/141‑3 (MQ=255)
      aTTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCt                                                                                                                                     <  2:95375/143‑1 (MQ=255)
        ttAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGa                                                                                                                                   <  2:29308/143‑1 (MQ=255)
        ttAGATAAAGATGGCGATGCAAACATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGa                                                                                                                                   <  1:353436/143‑1 (MQ=255)
         tAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGAc                                                                                                                                  <  2:34355/143‑1 (MQ=255)
          aGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCATATTTGCTGACt                                                                                                                                 <  2:950638/143‑1 (MQ=255)
              aaaTATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTc                                                                                                                             >  1:316302/1‑143 (MQ=255)
                  atGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAAc                                                                                                                         <  2:427768/143‑1 (MQ=255)
                     gCGATGAAAATATCAATCATTATTATAGAAACTTAGAACTAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATa                                                                                                                      <  2:446635/143‑1 (MQ=255)
                      cGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAg                                                                                                                     >  1:206697/1‑143 (MQ=255)
                        aTGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCAGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGtt                                                                                                                   <  1:889916/143‑1 (MQ=255)
                             aaTATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTa                                                                                                              >  2:370236/1‑143 (MQ=255)
                                 tcaatcaTTATTATAGATACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACa                                                                                                          <  1:585361/143‑1 (MQ=255)
                                   aatcaTTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTATTTCGACCAGTCAGCAAATATGCTGACTGGTCGAACAGAGTTTGCTAAACAGt                                                                                                        <  1:1049564/143‑1 (MQ=255)
                                      caTTATTATGGAACCGTAGAGATACAACACTCACAACCCAAAAGTAGAGCATCCCTCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGtcat                                                                                                     <  2:294482/143‑1 (MQ=255)
                                         tattatAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCac                                                                                                  >  2:302506/1‑143 (MQ=255)
                                           ttatAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactg                                                                                                >  2:79416/1‑143 (MQ=255)
                                              taGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgact                                                                                             >  2:360526/1‑143 (MQ=255)
                                               aGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCactgactg                                                                                            >  2:945990/1‑143 (MQ=255)
                                                       aGAAATAAAACACTCGCAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATgtcagt                                                                                    <  1:121610/143‑1 (MQ=255)
                                                          aaTAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCACATATCCTGACTGGTCGAACATACTTTGCTAAACAGTCATCACTGACTGCTCTCAGTcac                                                                                 >  1:976295/1‑143 (MQ=255)
                                                          aaTAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTcac                                                                                 >  1:540702/1‑143 (MQ=255)
                                                          aaTAAAACACCCACAACCCAAAAATAGAGCATCCGGCAAAAGATTGAGTTAACATCATCGTTAATTAAATCGACCAGTCATCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTcac                                                                                 <  1:241423/143‑1 (MQ=255)
                                                               aaCACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCAt                                                                            <  2:932696/143‑1 (MQ=255)
                                                                      acaACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGCCCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTCAACAGTCATCACTGACTGATGTCAGTCACACCAGCCACCAg                                                                     >  2:159534/1‑143 (MQ=255)
                                                                        aaCCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTACCATCATTGTTACTTAATTAGACCAGTCAGCAACTATGCTGACTGGTCGAACATAGTTTGATGAGATCCGGTCGCTGACTGATGTCAGTCACACCATCCGCCAGtt                                                                   <  2:917985/143‑1 (MQ=255)
                                                                        aaCCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCCAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGCTGTCAGTCCCACCACTCGCCAGtt                                                                   >  2:396673/1‑143 (MQ=255)
                                                                                         tCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAAt                                                  >  1:259837/1‑143 (MQ=255)
                                                                                           cGGCAAAAGAGTGAGTTAACATCATGGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTc                                                <  2:714202/143‑1 (MQ=255)
                                                                                             gCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTAGAACATAGTTTGCTAAACAGTCATCAATGACTGATGTCAGTCACACCACCCGCAAGCTCAATCACGCGAACAAATTcgg                                              >  2:787335/1‑141 (MQ=255)
                                                                                                 atgactgagtTAACATCATTGTAAACTAATGCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCt                                          <  1:81990/138‑1 (MQ=255)
                                                                                                      cgagtAAACATCATTTTTAATTAATCCGACCAGTCAGCAACTATGCTGAATGGTCTAACATAGTTTGCTAAACAGTTATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAAc                                     <  2:138787/142‑1 (MQ=255)
                                                                                                       gagtTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAAcc                                    <  1:747116/143‑1 (MQ=255)
                                                                                                         gctAACAGCATGGTTAATCAATTCGACCAGTCAGAAAATCTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCt                                  <  1:396673/141‑1 (MQ=255)
                                                                                                          ttAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAAcccgg                                 >  1:190364/1‑141 (MQ=255)
                                                                                                          ttAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtg                                 <  1:113043/143‑1 (MQ=255)
                                                                                                            aaCATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCtgtg                               <  1:737905/143‑1 (MQ=255)
                                                                                                                   tgcttaaataaTTCGAACAGTCAGCAAATATGCTGACTGGTCGAACATAGTGTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTc                        <  2:600375/140‑1 (MQ=255)
                                                                                                                   ctGTTAATTAGGTCGACCAGTCAGCACATATGCTGACTCGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTc                        <  2:497928/142‑1 (MQ=255)
                                                                                                                          ttaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCtttt                 <  1:370236/143‑1 (MQ=255)
                                                                                                                          ttaattCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCtttt                 <  1:175239/143‑1 (MQ=255)
                                                                                                                                    cAGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAg       >  1:521155/1‑143 (MQ=255)
                                                                                                                                    cAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAg       >  1:964598/1‑143 (MQ=255)
                                                                                                                                    cAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGAACTTCCCCTTTTGCCTGCTGAg       >  2:149961/1‑143 (MQ=255)
                                                                                                                                     aGTCAGCATATTTGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGCTCTTCCCCTTTTGCCTGCTGAGc      >  2:335173/1‑143 (MQ=255)
                                                                                                                                         aGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAc  <  2:540702/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                          
ATTGTTATTTAGATAAATATGGCGATGCAAATATCAATCATTATTATAGAAACTTAGAAATAAAACACTCACAACCCAAAAATAGAGCATCCGGCAAAAGAGTGAGTTAACATCATTGTTAATTAATTCGACCAGTCAGCAAATATGCTGACTGGTCGAACATAGTTTGCTAAACAGTCATCACTGACTGATGTCAGTCACACCATCCGCCAGTTCAATCAGGCGAATAAATTCCGCCCTGTAACCCTGTGGATCTTCCCCTTTTGCCTGCTGAGCCCAC  >  NC_000913/2383815‑2384094

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: