Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,584,502 T→C D662G (GAT→GGT)  ydeP ← putative oxidoreductase YdeP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,584,5020TC96.4% 88.0 / ‑5.1 28D662G (GAT→GGT) ydePputative oxidoreductase YdeP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (9/18);  total (9/19)
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

GCGTGAGCTGCGCTTACCGTCTGGCGTAAGCGCAATAAGATTAACTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGATCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCG  >  NC_000913/1584369‑1584637
                                                                                                                                     |                                                                                                                                       
gcgTGAGTTGCGCTTACCGTCTGGCGTAAGCGCAATAAGATTAACTCTTGCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATcca                                                                                                                                   <  2:450482/140‑1 (MQ=255)
            cTTACCGTCTGGCGTAAGCGCAATAAGATTAACTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAAtcg                                                                                                                       <  2:83039/140‑1 (MQ=255)
                 cGTCTGGCGTAAGCGCAATAAGATTAACTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCgtg                                                                                                                  <  2:70417/140‑1 (MQ=255)
                      ggCGTAAGCGCAATAAGATTAACTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCgtgt                                                                                                                 >  1:132387/1‑136 (MQ=255)
                      ggCGTAAGCGCAATAAGATTAACTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCgtgt                                                                                                                 <  2:132387/136‑1 (MQ=255)
                             gcgcAATAAGATTAACTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCg                                                                                                      <  2:223355/140‑1 (MQ=255)
                                cAATAAGATTAACTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTgg                                                                                                   <  2:297595/140‑1 (MQ=255)
                                         tAACTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACt                                                                                          >  2:246798/1‑140 (MQ=255)
                                          aaCTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTg                                                                                         <  2:324902/140‑1 (MQ=255)
                                           aCTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGt                                                                                        >  1:342342/1‑140 (MQ=255)
                                                    ccGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGAc                                                                               <  2:90761/140‑1 (MQ=255)
                                                         ttttACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCt                                                                          <  2:75199/140‑1 (MQ=255)
                                                                   cAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTaaa                                                                >  2:497933/1‑140 (MQ=255)
                                                                     aaTTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACg                                                              <  1:18994/140‑1 (MQ=255)
                                                                          tAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGaa                                                         <  1:497933/140‑1 (MQ=255)
                                                                              ttGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGATCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAgg                                                     <  2:200160/140‑1 (MQ=255)
                                                                              ttGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAgg                                                     <  1:125642/140‑1 (MQ=255)
                                                                                        cTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTaa                                           <  2:307636/140‑1 (MQ=255)
                                                                                         tCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAAc                                          >  2:512650/1‑140 (MQ=255)
                                                                                                gACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCtt                                   <  1:115772/140‑1 (MQ=255)
                                                                                                gACCACATCTCGTTGACCGAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCtt                                   >  1:174849/1‑140 (MQ=255)
                                                                                                                  gAATACCCCTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAg                 <  2:18629/140‑1 (MQ=255)
                                                                                                                         ccTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTAcc          >  1:355795/1‑140 (MQ=255)
                                                                                                                          cTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCt         >  2:80816/1‑140 (MQ=255)
                                                                                                                          cTCGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCt         >  1:146470/1‑140 (MQ=255)
                                                                                                                            cGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGa       <  1:406118/140‑1 (MQ=255)
                                                                                                                            cGATAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGa       <  1:246798/140‑1 (MQ=255)
                                                                                                                              aTAGCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACg     >  2:364905/1‑140 (MQ=255)
                                                                                                                                 gCGACCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCg  <  1:63645/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
GCGTGAGCTGCGCTTACCGTCTGGCGTAAGCGCAATAAGATTAACTCTTTCGCCGTTTTTTACACGGCAAATTTTAGCTTGTTTAGCACTCATAAAGACCACATCTCGTTGACCGAATACCCCTCGATAGCGATCATCCATACCATAAATCGTCGTGTTGTACTGATCGTGGCTGCGTACTGTCGCCATGACCAGCTTACTGTTAAACGCTGAAGAGGGATCTTCTAACAGCCCTTTGCTGGTAATGAAATTAGCCTTACCTGACGGCG  >  NC_000913/1584369‑1584637

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑