Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,357,671 T→G intergenic (‑3/+131) ymjA ← / ← puuP DUF2543 domain‑containing protein YmjA/putrescine:H(+) symporter PuuP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,357,6710TG96.9% 104.6 / NA 32intergenic (‑3/+131)ymjA/puuPDUF2543 domain‑containing protein YmjA/putrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (14/17);  minor base A (1/0);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATCAAACGGGTGAGCAGTAACTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCATCCGACTA  >  NC_000913/1357535‑1357801
                                                                                                                                        |                                                                                                                                  
ttCATCAAACGGGTGAGCAGTAACTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTg                                                                                                                                 >  2:98598/1‑140 (MQ=255)
ttCATCAAACGGGTGAGCAGTAACTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTg                                                                                                                                 >  1:593979/1‑140 (MQ=255)
          gggTGAGCAGTAACTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGGGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGt                                                                                                                       <  2:107106/140‑1 (MQ=255)
          gggTGAGCAGTAACTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGt                                                                                                                       <  2:614800/140‑1 (MQ=255)
              gAGCAGTAACTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTGATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAAt                                                                                                                   <  2:720989/140‑1 (MQ=255)
                   gTAACTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCt                                                                                                              <  2:509059/140‑1 (MQ=255)
                    tAACTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTa                                                                                                             <  2:191164/140‑1 (MQ=255)
                      aCTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATTGTGGTTCATTAATTGCCCCTGTGGTTAATGCGCTAAg                                                                                                           >  1:392510/1‑140 (MQ=255)
                      aCTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGcccc                                                                                                                             <  1:566156/122‑1 (MQ=255)
                      aCTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGcccc                                                                                                                             >  2:566156/1‑122 (MQ=255)
                           aaGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATaa                                                                                                      <  1:98598/140‑1 (MQ=255)
                              tagtGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCGTTAGTTGACCCGGGGGTTAATGCGCTAAGCATAAggg                                                                                                   >  1:226946/1‑140 (MQ=255)
                              tagtGGGCCAGCGGTGTGCGTTAAGCCTCTGTGACGGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAggg                                                                                                   <  2:392510/140‑1 (MQ=255)
                                    gCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGgatg                                                                                             >  2:7696/1‑140 (MQ=255)
                                    gCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTAGTTCATTAATTGACACTGGGGTTAAAGCGCTAAGCATAACGGCGgatg                                                                                             >  2:404315/1‑140 (MQ=255)
                                     ccAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTCATGCGCTAAACATAAGGGCGgatga                                                                                            >  1:309489/1‑140 (MQ=255)
                                        gCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAg                                                                                         >  1:781082/1‑140 (MQ=255)
                                        gCGGTGTGCGTTCAGCCCCTGCCACAGGCTCGGCGCATTTGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAg                                                                                         <  2:409196/140‑1 (MQ=255)
                                        gCGGTGTGCGTTCAGCCACTGCGACAGGCTCTGCGCATTGGGTCGCGTGCTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAg                                                                                         <  2:309489/140‑1 (MQ=255)
                                            tgtgCGTTCAGCCTCCGCGACACGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACt                                                                                     <  1:251501/140‑1 (MQ=255)
                                              tgCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCg                                                                                                 <  1:739578/126‑1 (MQ=255)
                                              tgCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCg                                                                                                 >  2:739578/1‑126 (MQ=255)
                                              tgCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCg                                                                                   >  1:782690/1‑140 (MQ=255)
                                              tgCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCg                                                                                   >  1:511516/1‑140 (MQ=255)
                                                  ttCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGAt                                                                               <  1:595632/140‑1 (MQ=255)
                                                   tCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAACCATTAGCGCGGATTATAGCACTCGCGATc                                                                              >  1:352059/1‑140 (MQ=255)
                                                    ctgCCTCTGAGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTATAGCGGGCTATCGTGGTACATCAGTTGCCCCTAGGTTTAATGCGCTAAGCATAAGGGCGGATGATGGCTCTCGTGATcc                                                                             <  2:860158/138‑1 (MQ=255)
                                                     aGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCg                                                                            <  2:593979/140‑1 (MQ=255)
                                                          ctGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCGAAGAATACGGGCGGATGATAGGACTCGTGATCCGGGTCa                                                                       >  2:235450/1‑140 (MQ=255)
                                                                   gCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCttt                                                              >  2:808789/1‑140 (MQ=255)
                                                                                  tCGCGTACTCATCGGAAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGACGATAGCACTCGTGGTCCGGGTCATACCTCTTTTTATAGTTAAAggcg                                               <  2:147124/140‑1 (MQ=255)
                                                                                                    tatCAAAATATTTTAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTcagg                             <  1:808789/140‑1 (MQ=255)
                                                                                                      tCAAAATATTTTAGCGGGATATAGTGGTTCAATAGGTGCGCCTGGGGTTAATTCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTcaggcg                           <  2:352059/140‑2 (MQ=255)
                                                                                                               tttAGCGGGATATCGTGGTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATgc                  >  2:858024/1‑140 (MQ=255)
                                                                                                                            cGTGGTTCATTAGTTACCCCTGAGGATAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCATCCGa     <  1:404315/140‑1 (MQ=255)
                                                                                                                               ggTTCATTAGTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCATCCGACTa  <  2:511516/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
TTCATCAAACGGGTGAGCAGTAACTGGAAGTAGTGGGCCAGCGGTGTGCGTTCAGCCTCTGCGACAGGCTCTGCGCATTCGGTCGCGTACTCATCGGCAATATCAAAATATTTTAGCGGGATATCGTGGTTCATTATTTGCCCCTGGGGTTAATGCGCTAAGCATAAGGGCGGATGATAGCACTCGTGATCCGGGTCATACATCTTTTTATAGTTAATGGCGCGGCGCGTTACCCTCAGGCAGGATAATGCGCCGCGCATCCGACTA  >  NC_000913/1357535‑1357801

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