Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,797,320 C→A P37H (CCC→CAC)  ygaV → putative DNA‑binding transcriptional regulator YgaV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,797,3200CA95.7% 149.8 / 0.7 46P37H (CCC→CAC) ygaVputative DNA‑binding transcriptional regulator YgaV
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (25/19);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

ACTTAAATAAATATGACTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCCCGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAACGCATTCTATATTCCATTAA  >  NC_000913/2797199‑2797455
                                                                                                                         |                                                                                                                                       
aCTTAAATAAATATGACTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGgaga                                                                                                                       >  1:447265/1‑140 (MQ=255)
         aataTGACTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGAcgcg                                                                                                              <  1:108051/140‑1 (MQ=255)
         aataTGACTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGAcgcg                                                                                                              <  2:162067/140‑1 (MQ=255)
         aataTGACTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGAcgcg                                                                                                              <  1:118041/140‑1 (MQ=255)
           tatGACTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCa                                                                                                            >  1:380280/1‑140 (MQ=255)
               aCTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTAc                                                                                                        >  2:151488/1‑140 (MQ=255)
                cTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTAcc                                                                                                       >  2:446956/1‑140 (MQ=255)
                cTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTAcc                                                                                                       >  2:425993/1‑140 (MQ=255)
                cTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTAcc                                                                                                       >  2:43384/1‑140 (MQ=255)
                   aaCTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCgg                                                                                                     <  2:158850/139‑1 (MQ=255)
                   aaCTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCgg                                                                                                     >  1:158850/1‑139 (MQ=255)
                   aaCTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGa                                                                                                    >  2:56851/1‑140 (MQ=255)
                   aaCTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGa                                                                                                    >  1:191472/1‑140 (MQ=255)
                     cTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACt                                                                                                  >  2:158383/1‑140 (MQ=255)
                        gcgcAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAg                                                                                               >  2:403449/1‑140 (MQ=255)
                                   ggCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGa                                                                                    <  1:293367/140‑1 (MQ=255)
                                         tGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATcac                                                                              >  1:31055/1‑140 (MQ=255)
                                         tGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATcac                                                                              >  1:355264/1‑140 (MQ=255)
                                         tGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATcac                                                                              >  1:388831/1‑140 (MQ=255)
                                         tGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATcac                                                                              >  1:113222/1‑140 (MQ=255)
                                              aaCAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCAt                                                                         >  2:207486/1‑140 (MQ=255)
                                                  ggCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATctcg                                                                     <  2:61492/140‑1 (MQ=255)
                                                     aGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATctcgctc                                                                  >  2:406152/1‑140 (MQ=255)
                                                             tATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCgg                                                          >  2:228600/1‑140 (MQ=255)
                                                             tATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTACGACATCACAGCATCTCGCTCGTATGCggg                                                         >  2:437848/1‑140 (MQ=255)
                                                              aTTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCggg                                                         <  1:406152/140‑1 (MQ=255)
                                                                  aaaGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGa                                                     <  1:151488/140‑1 (MQ=255)
                                                                  aaaGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGa                                                     >  1:106034/1‑140 (MQ=255)
                                                                           aGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTAt                                            >  1:14885/1‑140 (MQ=255)
                                                                              cACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGa                                         <  1:446956/140‑1 (MQ=255)
                                                                                         aTTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACggg                              <  1:403449/140‑1 (MQ=255)
                                                                                         aTTGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACggg                              <  2:31055/140‑1 (MQ=255)
                                                                                           tGCTGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGAt                            >  2:82705/1‑140 (MQ=255)
                                                                                           tGCTGATTCTGTGCATGCTTAGCGGCTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGAt                            >  1:307694/1‑140 (MQ=255)
                                                                                              tGATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGcc                         <  1:81764/140‑1 (MQ=255)
                                                                                               gATTCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGccc                        >  1:396572/1‑140 (MQ=255)
                                                                                                  tCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAAc                     <  1:425993/140‑1 (MQ=255)
                                                                                                  tCTGTGCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAAc                     <  2:447265/140‑1 (MQ=255)
                                                                                                     gtgCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAACGCa                  <  2:191472/140‑1 (MQ=255)
                                                                                                     gtgCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAACGCa                  <  1:82705/140‑1 (MQ=255)
                                                                                                     gtgCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAACGCa                  <  1:56851/140‑1 (MQ=255)
                                                                                                     gtgCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAACGCa                  <  1:325203/140‑1 (MQ=255)
                                                                                                       gCATGCTTAGCGGTTCCCACGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAACGCAtt                <  2:372768/140‑1 (MQ=255)
                                                                                                                     tCCCCCGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAACGCa                  >  2:365939/1‑124 (MQ=255)
                                                                                                                     tCCCCCGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAACGCa                  <  1:365939/124‑1 (MQ=255)
                                                                                                                     tCCCACGGCACCAGCGCGGGAGAGCTGACGCGAATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAACGCATTCTATATTCCATTaa  >  1:423496/1‑140 (MQ=255)
                                                                                                                         |                                                                                                                                       
ACTTAAATAAATATGACTGAACTCGCGCAATTACAGGCCAGTGCCGAACAGGCAGCGGCCTTATTGAAAGCAATGAGCCACCCTAAACGATTGCTGATTCTGTGCATGCTTAGCGGTTCCCCCGGCACCAGCGCGGGAGAGCTGACGCGCATTACCGGACTGAGTGCCTCTGCGACATCACAGCATCTCGCTCGTATGCGGGACGAAGGGCTTATCGACAGCCAACGGGATGCCCAACGCATTCTATATTCCATTAA  >  NC_000913/2797199‑2797455

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