Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 635,915 A→G Y219H (TAT→CAT)  ybdN ← putative PAPS reductase/DUF3440 domain‑containing protein YbdN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913635,9150AG98.0% 174.2 / ‑5.1 51Y219H (TAT→CAT) ybdNputative PAPS reductase/DUF3440 domain‑containing protein YbdN
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (21/29);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGG  >  NC_000913/635778‑636051
                                                                                                                                         |                                                                                                                                        
ttCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGcc                                                                                                                                        <  2:365740/140‑1 (MQ=255)
ttCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGcc                                                                                                                                        <  2:138414/140‑1 (MQ=255)
  cGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAg                                                                                                                                      <  1:103396/140‑1 (MQ=255)
       aCGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCa                                                                                                                                 >  1:818588/1‑140 (MQ=255)
          tGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGat                                                                                                                              >  2:510137/1‑140 (MQ=255)
            aTAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGatat                                                                                                                            >  1:771390/1‑140 (MQ=255)
             tAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATc                                                                                                                           >  1:351911/1‑140 (MQ=255)
                agCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGc                                                                                                                        <  2:870614/140‑1 (MQ=255)
                agCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGc                                                                                                                        <  1:558357/140‑1 (MQ=255)
                  ccATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCa                                                                                                                      >  2:86217/1‑140 (MQ=255)
                       aTCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTc                                                                                                                 <  2:329198/140‑1 (MQ=255)
                       aTCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTc                                                                                                                 <  1:912150/140‑1 (MQ=255)
                       aTCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTc                                                                                                                 <  2:771390/140‑1 (MQ=255)
                       aTCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTc                                                                                                                 <  1:466446/140‑1 (MQ=255)
                       aTCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTc                                                                                                                 <  2:920282/140‑1 (MQ=255)
                                   gCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATg                                                                                                     >  2:895390/1‑140 (MQ=255)
                                    cTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATgg                                                                                                    >  1:337535/1‑140 (MQ=255)
                                        gggCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTa                                                                                                <  2:110659/140‑1 (MQ=255)
                                              aaaGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa                                                                                          <  1:510137/140‑1 (MQ=255)
                                              aaaGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa                                                                                          <  1:895390/140‑1 (MQ=255)
                                                     cGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTa                                                                                   >  2:463986/1‑140 (MQ=255)
                                                       cAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATg                                                                                 <  2:828064/140‑1 (MQ=255)
                                                       cAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATg                                                                                 <  2:583706/140‑1 (MQ=255)
                                                         aaTTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGAc                                                                               <  2:818588/140‑1 (MQ=255)
                                                             cGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGcc                                                                           <  1:249335/140‑1 (MQ=255)
                                                                      cGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTg                                                                      >  2:953379/1‑136 (MQ=255)
                                                                      cGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTg                                                                      <  1:953379/136‑1 (MQ=255)
                                                                        cAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCg                                                                <  1:63414/140‑1 (MQ=255)
                                                                        cAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCg                                                                <  2:911362/140‑1 (MQ=255)
                                                                          gagGAACGCCTCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAg                                                             >  1:332138/1‑140 (MQ=255)
                                                                             gAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTc                                                           <  1:220743/140‑1 (MQ=255)
                                                                             gAACGCCTCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTcc                                                          <  2:332138/140‑1 (MQ=255)
                                                                              aaCGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTcc                                                          <  2:165236/140‑1 (MQ=255)
                                                                               aCGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCa                                                         <  1:86217/140‑1 (MQ=255)
                                                                                                   gTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAAc                                     <  2:351911/140‑1 (MQ=255)
                                                                                                    ttATACAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACg                                    <  2:126098/140‑1 (MQ=255)
                                                                                                      ataCAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACgtt                                  >  2:125253/1‑140 (MQ=255)
                                                                                                       taCAGTGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACgttg                                 >  2:506200/1‑140 (MQ=255)
                                                                                                            tGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGtttat                            >  1:102756/1‑140 (MQ=255)
                                                                                                            tGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGtttat                            >  1:14223/1‑140 (MQ=255)
                                                                                                            tGGGTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGtttat                            <  1:506200/140‑1 (MQ=255)
                                                                                                              ggTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACgttg                                 >  1:299853/1‑133 (MQ=255)
                                                                                                              ggTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACgttg                                 <  2:299853/133‑1 (MQ=255)
                                                                                                              ggTTGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGtttattt                          >  2:7213/1‑140 (MQ=255)
                                                                                                                 tGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTaag                       >  1:78050/1‑139 (MQ=255)
                                                                                                                 tGCAGAGACTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTaaa                       >  2:405870/1‑140 (MQ=255)
                                                                                                                        aCTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACtggcga                >  2:857574/1‑140 (MQ=255)
                                                                                                                        aCTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACtggcga                >  1:564782/1‑140 (MQ=255)
                                                                                                                        aCTTTGATGATTTGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACtggcga                >  2:682544/1‑140 (MQ=255)
                                                                                                                                  tttGCATGCCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACaaa      <  2:564782/140‑1 (MQ=255)
                                                                                                                                      cATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACgg  <  2:614564/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
TTCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAAACTGGCGATGGCGACAAAACGG  >  NC_000913/635778‑636051

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