Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,086,636 C→T 7.1% A143T (GCC→ACC)  plaP ← putrescine importer, low affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,086,6360CT 7.1% 178.8 / 3.8 56A143T (GCC→ACC) plaPputrescine importer, low affinity
Reads supporting (aligned to +/- strand):  ref base C (29/23);  new base T (2/2);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

ACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGA  >  NC_000913/2086490‑2086782
                                                                                                                                                  |                                                                                                                                                  
aCGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCgg                                                                                                                                                  >  3:411889/1‑149 (MQ=255)
        tCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAg                                                                                                                                          >  1:121087/1‑149 (MQ=255)
        tCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAg                                                                                                                                          >  3:355371/1‑149 (MQ=255)
             aGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGccac                                                                                                                                     <  4:411926/149‑1 (MQ=255)
                     aGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGcca                                                                                                                             <  4:279575/149‑1 (MQ=255)
                      gcgcCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGccac                                                                                                                            <  2:428410/149‑1 (MQ=255)
                                cttcaAATACGCCATAAACAACTATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATc                                                                                                                  <  3:374850/149‑1 (MQ=255)
                                cttcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATc                                                                                                                  >  3:224875/1‑149 (MQ=255)
                                cttcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATc                                                                                                                  >  4:406614/1‑149 (MQ=255)
                                 ttcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATcc                                                                                                                 >  3:273183/1‑149 (MQ=255)
                                  tcaAATACGCCATCAACAGCCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCa                                                                                                                <  1:365676/149‑1 (MQ=255)
                                  tcaAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCa                                                                                                                <  3:17748/149‑1 (MQ=255)
                                         cGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACgatggg                                                                                                         <  4:355371/149‑1 (MQ=255)
                                          gCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACgatggga                                                                                                        >  3:381432/1‑149 (MQ=255)
                                            cATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACgatgggatg                                                                                                      >  1:399894/1‑149 (MQ=255)
                                            cATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACgatgggatg                                                                                                      <  4:273183/149‑1 (MQ=255)
                                                   aacCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCa                                                                                               >  1:1542/1‑149 (MQ=255)
                                                    acCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCAc                                                                                              >  4:362241/1‑149 (MQ=255)
                                                           ccAGAATCACCGCGATCAGCACTACCTGCGACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGTGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCt                                                                                       >  2:353355/1‑149 (MQ=255)
                                                           ccAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGTGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCt                                                                                       >  4:407924/1‑149 (MQ=255)
                                                           ccAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCt                                                                                       <  3:362241/149‑1 (MQ=255)
                                                             aGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAGGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTc                                                                                     >  4:370835/1‑149 (MQ=255)
                                                             aGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTc                                                                                     <  3:262358/149‑1 (MQ=255)
                                                             aGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTc                                                                                     >  2:257409/1‑149 (MQ=255)
                                                              gAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCa                                                                                    >  4:73853/1‑149 (MQ=255)
                                                                   aCCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATa                                                                               <  2:432424/149‑1 (MQ=255)
                                                                   aCCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATa                                                                               <  4:381432/149‑1 (MQ=255)
                                                                         aTCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGAtttt                                                                         >  4:109766/1‑149 (MQ=255)
                                                                          tCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTc                                                                        <  4:61285/149‑1 (MQ=255)
                                                                            aGCACTACCTGCAACACGACGATTACAGTGTTGAGGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGc                                                                      <  3:370835/149‑1 (MQ=255)
                                                                             gCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGcc                                                                     >  1:39224/1‑149 (MQ=255)
                                                                                     tGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAAt                                                             >  2:454546/1‑149 (MQ=255)
                                                                                        aaCACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTc                                                                                     >  2:267779/1‑122 (MQ=255)
                                                                                        aaCACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTc                                                                                     <  1:267779/122‑1 (MQ=255)
                                                                                           acgacgATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGat                                                       >  1:65843/1‑149 (MQ=255)
                                                                                           acgacgATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGat                                                       >  4:112972/1‑149 (MQ=255)
                                                                                                          ttGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAg                                        >  4:370729/1‑149 (MQ=255)
                                                                                                           tGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGa                                       <  1:454546/149‑1 (MQ=255)
                                                                                                               gTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGt                                   >  4:463844/1‑149 (MQ=255)
                                                                                                                  cGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGa                                <  1:257409/149‑1 (MQ=255)
                                                                                                                     tACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCa                             >  2:394648/1‑149 (MQ=255)
                                                                                                                            ttCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGagaaga                      <  3:406614/149‑1 (MQ=255)
                                                                                                                                gACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCaa                  >  4:251184/1‑149 (MQ=255)
                                                                                                                                gACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCaa                  >  2:354021/1‑149 (MQ=255)
                                                                                                                                gACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCaa                  >  2:381577/1‑149 (MQ=255)
                                                                                                                                      gCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGtt                                                          <  2:313180/103‑1 (MQ=255)
                                                                                                                                      gCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGtt                                                          >  1:313180/1‑103 (MQ=255)
                                                                                                                                          gTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATaa        >  4:8682/1‑149 (MQ=255)
                                                                                                                                           ttAAAGGTGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCATCCCACCATaaa       <  1:353355/149‑1 (MQ=255)
                                                                                                                                           ttAAAGGTGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATaaa       <  3:407924/149‑1 (MQ=255)
                                                                                                                                             aaaGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGtt                                                          <  2:143662/96‑1 (MQ=255)
                                                                                                                                             aaaGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGtt                                                          >  1:143662/1‑96 (MQ=255)
                                                                                                                                             aaaGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGc     <  3:251184/149‑1 (MQ=255)
                                                                                                                                              aaGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGcc    <  2:65843/149‑1 (MQ=255)
                                                                                                                                                ggCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGa  <  4:224875/149‑1 (MQ=255)
                                                                                                                                                ggCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGa  >  1:192007/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATTTCAGACTACGCAGGTTAAAGGCGGTCATAAAGGCCACCAGCGCCACCACAAACATCCACGATGGGATGGAAGGCACCAGAGCTTCAAAATAGATTTTCGCCAGTAGAATGTTGATCATCGGCGCGAACAGATAGTCGAGCAGAGAAGACCAACCCACCATAAAGCCGA  >  NC_000913/2086490‑2086782

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