Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,858,955 (C)5→4 coding (1275/1380 nt) ydjK ← putative transporter YdjK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,858,9510C.93.9% 125.0 / 2.3 33coding (1279/1380 nt)ydjKputative transporter YdjK
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base . (20/11);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GCTTTTGAATTTTTCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGA  >  NC_000913/1858817‑1859088
                                                                                                                                      |                                                                                                                                         
gCTTTTGAATTTTTCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGccccag                                                                                                                                      <  1:362199/140‑3 (MQ=255)
 cTTTTGAATTTTTCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGccccaga                                                                                                                                     >  1:120661/1‑137 (MQ=255)
 cTTTTGAATTTTTCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGccccaga                                                                                                                                     <  2:343965/140‑4 (MQ=255)
 cTTTTGAATTTTTCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGccccaga                                                                                                                                     >  1:24735/1‑137 (MQ=255)
 cTTTTGAATTTTTCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATAAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGccccaga                                                                                                                                     <  2:113160/140‑4 (MQ=255)
  ttttGAATTTTTCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGccccagaa                                                                                                                                    <  2:77487/140‑5 (MQ=255)
  ttttGAATTTTTCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGccccagaa                                                                                                                                    <  2:608975/140‑5 (MQ=255)
   tttGAATTTTTCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCagaag                                                                                                                                  >  1:98590/1‑140 (MQ=255)
            ttCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGAc                                                                                                                         <  1:531466/140‑1 (MQ=255)
                          ccccGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGTTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCAt                                                                                                           <  2:501568/139‑1 (MQ=255)
                           cccGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAGACCGCCCGAGAAGACTAAAGAATCCCGTAACGCCATa                                                                                                          >  1:212510/1‑140 (MQ=255)
                               tattatTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAacta                                                                                                      <  2:24735/140‑1 (MQ=255)
                                    tttattGGCTACTGCATCAATAGTTAAACTTTCAAAGGACACACCTTTTGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGag                                                                                                 <  2:214079/140‑1 (MQ=255)
                                      tattGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGCAGAATAAAGACTCCTGTAACGCCATAACTACTGagca                                                                                               >  1:258008/1‑140 (MQ=255)
                                         tGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAgcagca                                                                                            >  1:119836/1‑140 (MQ=255)
                                         tGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAgcagca                                                                                            >  2:606384/1‑140 (MQ=255)
                                         tGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAgcagca                                                                                            >  1:454451/1‑140 (MQ=255)
                                          ggCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCAc                                                                                           <  2:485139/140‑1 (MQ=255)
                                           gCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGag                                                                                                 >  1:401563/1‑133 (MQ=255)
                                           gCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGag                                                                                                 <  2:401563/133‑1 (MQ=255)
                                              aCTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCa                                                                                       <  1:33463/140‑1 (MQ=255)
                                                     cAATATTTAAACTTTCAACGAAGGCTCTTTTGGTTTCGATTCAAATGGTGGCGATAGCATTGGCGACAAAAATTAAAACCGCCCCAGAAGAATAAAGCCTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCAt                                                                                <  2:119836/140‑1 (MQ=255)
                                                       aTACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATaa                                                                              >  1:58960/1‑140 (MQ=255)
                                                            tAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGAATCCCGTAAGGCCATAAGTACTGAGCTGCACAGCAACGGCATAAGGTGc                                                                         >  2:191952/1‑140 (MQ=255)
                                                                   tCAATGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATa                                                                  >  1:123233/1‑140 (MQ=255)
                                                                   tCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATa                                                                  >  1:480781/1‑140 (MQ=255)
                                                                   tCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATa                                                                  >  2:649296/1‑140 (MQ=255)
                                                                    cAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATAc                                                                 >  1:489478/1‑140 (MQ=255)
                                                                          agaCACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGa                                                           <  2:444906/140‑1 (MQ=255)
                                                                          agaCACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACAGa                                                           >  2:409767/1‑140 (MQ=255)
                                                                                   tGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTa                                                  <  2:454451/140‑1 (MQ=255)
                                                                                    ggTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTAc                                                 >  1:388294/1‑140 (MQ=255)
                                                                                    ggTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTAc                                                 >  2:532046/1‑140 (MQ=255)
                                                                                        tCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCa                                             >  2:431152/1‑140 (MQ=255)
                                                                                        tCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCa                                             >  2:352493/1‑140 (MQ=255)
                                                                                          gATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCa                                             >  2:492591/1‑138 (MQ=255)
                                                                                          gATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCa                                             <  1:492591/138‑1 (MQ=255)
                                                                                                    gtggCGATAGCAATTGCGACAATAATTGAAACCGCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGAAATACCACTGAATCACCCTACTGCATTCGCCAGACCg                                 >  1:267222/1‑140 (MQ=255)
                                                                                                       gCGATAGCAATTGCGACAATAATTGAAACCGCCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGa                               >  2:369620/1‑140 (MQ=255)
                                                                                                                                    cGCCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGGGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCATACCGGAGCCACGGAGATTGGCCTCTGTCGGCCAGa  >  2:387142/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                         
GCTTTTGAATTTTTCATTTTGAATCCACCCGTATTATTTATTGGCTACTGCATCAATACTTAAACTTTCAACGGAGACACCTTTGGTTTCGATTCCAATGGTGGCGATAGCAATTGCGACAATAATTGAAACCGCCCCCAGAAGAATAAAGACTCCCGTAACGCCATAACTACTGAGCAGCACTGCAACGGCATAAGGTGCGGCAATACCACTGATTCGCCCTACCGCATTCGCCAGACCGGAGCCACGGAGTTTGGCCTCTGTCGGCCAGA  >  NC_000913/1858817‑1859088

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