Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,182,771 G→C 9.2% A20G (GCC→GGC)  potD ← spermidine preferential ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,182,7710GC 9.2% 119.0 / 3.8 48A20G (GCC→GGC) potDspermidine preferential ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base C (2/3);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

GGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAAT  >  NC_000913/1182629‑1182913
                                                                                                                                              |                                                                                                                                              
ggCGTACCTGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAgcgg                                                                                                                                                <  2:106656/143‑1 (MQ=255)
ggCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAgcgg                                                                                                                                                <  1:126088/143‑1 (MQ=255)
  cGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAATGTGAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAgcggcg                                                                                                                                              <  1:867968/143‑1 (MQ=255)
    tACATGGTTTCGTTCGACTCGTAAGTCGAAGAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCt                                                                                                                                            <  2:482955/143‑1 (MQ=255)
           tttgtttctACTCGTAGGTCGGTTAGATAACCTTATTACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCACGccg                                                                                                                                     <  2:18100/138‑2 (MQ=255)
              cgttcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCATGAGCGCCGCTCACGccgatg                                                                                                                                  <  2:468514/143‑5 (MQ=255)
                 tcgACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGcc                                                                                                                               >  2:37682/1‑143 (MQ=255)
                          aaGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTAATCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGACCAGTGCCAGAGCAccc                                                                                                                      >  2:642277/1‑143 (MQ=255)
                            gTCGAATAGATAACCTTAATCCCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACccgc                                                                                                                    <  1:918554/143‑1 (MQ=255)
                               gAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCcgccgc                                                                                                                 >  1:608948/1‑143 (MQ=255)
                               gAATAGATAACCTTAATACCGAGTGCTAAAGAGAAATGTACAAGCAGTGGTGGCGGCACCCACAAAATCAAATTCTAAAAAAACAGCGAGTCGTTGTCATCGACCTGAGTGCAGCTCACGCGCAGTTCAACAGCACccgccag                                                                                                                 >  2:1246618/1‑141 (MQ=255)
                                 aTAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGa                                                                                                               <  1:48097/143‑1 (MQ=255)
                                  tAGATAACCTTACGACCGGTTTGTTTGGTGCACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAg                                                                                                              <  2:1281370/143‑1 (MQ=255)
                                             aaTACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTg                                                                                                   <  1:1123085/143‑1 (MQ=255)
                                              aTACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGa                                                                                                  <  2:1251698/143‑1 (MQ=255)
                                                  cGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCCCGTCCTCGGTCCAGTTTTAGACATACCGCGTGTTGTTGTCATAGGCGTGCGCGGCGCCCATCCCCAGTGCCATAGCAACCGCCGCGAGCAGGTcgcgctaccct                                                                                              >  2:970476/1‑136 (MQ=255)
                                                    gTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCAttt                                                                                            >  2:1221878/1‑143 (MQ=255)
                                                          ttGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGACGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGCGGCGTGACCATTTTTTCAt                                                                                      >  1:136229/1‑143 (MQ=255)
                                                            ggTGAACTGTTCAAGCCGTCCTGGCGGCACGTACTCGGTCCAGCTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGCCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCAttt                                                                                    <  1:1153801/143‑1 (MQ=255)
                                                            ggTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCAttt                                                                                    >  2:156077/1‑143 (MQ=255)
                                                                  cTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCAGCCGAAGTGCCGGAGCACCCGCCGCGAGCAGGTGGCGTGAGCAGTTTTTCATTTAACGTc                                                                              >  1:727579/1‑143 (MQ=255)
                                                                                   gcggcACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCCCCCGCCGCGAGCAGGGGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCCTTCg                                                             >  1:650981/1‑143 (MQ=255)
                                                                                     ggcACGGACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGtt                                                           <  1:1000647/143‑1 (MQ=255)
                                                                                      gcgtgttcTCGGTCCGGTTGGAGAAATACAGCGTGTTGTTGGCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGTGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGttt                                                          <  2:236440/136‑1 (MQ=255)
                                                                                      gcACGTGGTCGGTCCGGTTTTGGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGTGCACCGGCCGCGAGCAGGTGGCGTGACCATTTTTGCATTTAACGTCCCCTGTGTTACCTGTCGttt                                                          <  2:196239/143‑1 (MQ=255)
                                                                                         cggaCTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTAt                                                       <  1:1164201/140‑1 (MQ=255)
                                                                                             cTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACg                                                   <  1:791813/143‑1 (MQ=255)
                                                                                              tCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGa                                                  <  2:596714/143‑1 (MQ=255)
                                                                                              tCGGTCCAGTTGTAGAAATACAGCGAGTTGTTATCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCCGCGGGCGGCACCAGTTTATCATCTAGAGACCGCAGTGCTACCTTTCACTTGATCAcaa                                                  >  2:750867/1‑141 (MQ=255)
                                                                                               cGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAg                                                 <  1:1221878/143‑1 (MQ=255)
                                                                                                 gTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCa                                               >  1:31796/1‑143 (MQ=255)
                                                                                                  ttcAGTTGTTGAAATACAGCGTGTTGTTCTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCTAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTTTTATCTTTCGTTTTATCACGAGCaa                                              <  2:1122316/141‑1 (MQ=255)
                                                                                                    cAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATa                                            >  1:508870/1‑143 (MQ=255)
                                                                                                          tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGGCCATTTTTTCATTTAACGTCCCCCGTGTTCCCTTTCGTTTTATCACACGCAAACAgctgc                                      >  2:789576/1‑142 (MQ=255)
                                                                                                          tAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctgg                                      >  1:620117/1‑143 (MQ=255)
                                                                                                            gAAATACAGCGTGTTGTTGTCATGGGCGTGAGAGGCGCTCATGCCCAGTGACAGAGCACACGCAGCGGGCAGGTGGCGTGACCATTTTTTCAATTACGAAGCCATGAGTTATCTTGCGATTTATCTCGGGGAGAAAgctggca                                    >  2:505258/1‑142 (MQ=255)
                                                                                                            gAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAgctggct                                    >  2:24910/1‑143 (MQ=255)
                                                                                                                 aCAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGGGCACCCGCAGCGAACAGGTGGCGGAACCATTTTTTCATTTAGCGCCCCGGGAGTAACCTTTCGCTGTATCACGAGCAATAAGCTGGCCGGCag                               >  2:302067/1‑142 (MQ=255)
                                                                                                                   aGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATa                             >  2:261452/1‑143 (MQ=255)
                                                                                                                     cGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCGGGTGGCGGGACCATTTTTTCATTTAACGTCCCCGGTGTTACCTTTCGTTTTATCACGCGCACTACGCTGCCCGACAATGAc                           >  2:579260/1‑143 (MQ=255)
                                                                                                                     cGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAAc                           <  1:607392/143‑1 (MQ=255)
                                                                                                                         ttgttgcCATCGCCGTGAGCGGCGCCCATGCCCAGTGCCAAAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATc                       <  1:800005/143‑1 (MQ=255)
                                                                                                                            ttgtCATCGGCGTGAGCGGGGCTGATGCCCGGTGCCAGAGCAGCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCAcc                    <  1:252319/143‑1 (MQ=255)
                                                                                                                             tgtCACGGGCGGGGGCTGCGTTCAGGCCAACTGCCAGAGAACCCGACGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTCTCGTTATATCTAGAGCAATAAGCTGGCTGGCAATAAACATCACCa                   <  2:516239/143‑1 (MQ=255)
                                                                                                                               tgttCGGCGTGGGGGGCGCTCGTGCCCAGTGTCAGAGCGTCGGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAATCTGGCTGGCAATAACCATCACCAGc                 <  2:1004808/140‑1 (MQ=255)
                                                                                                                                 aTCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGGAGGTGGCGTGACCATTTTTTCAGTTAACGGCCGCTGGGTTACCTTTCGTTTTATCACGAGGAATAAACTCGCTAGCTGGAACCATCACCAGCgc               >  2:607940/1‑142 (MQ=255)
                                                                                                                                  tcggcgtcggcgGCTCCCACGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCAGTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGAc              <  1:92507/143‑1 (MQ=255)
                                                                                                                                  tCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTACCGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGAc              <  1:95115/143‑1 (MQ=255)
                                                                                                                                   cGGCGTGAGCGGCGCTCATGCCCAGGGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACa             <  1:460549/143‑1 (MQ=255)
                                                                                                                                    ggCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAg            <  2:792213/143‑1 (MQ=255)
                                                                                                                                      catGAGCGCCGCTCACGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCGCTGTATTACCTTTCGTTTTATCACGAGCAATACGCTGGCTGGCAATAACCAGCACCAGCGACAGCa          >  1:389946/3‑143 (MQ=255)
                                                                                                                                             ggcgCTCGTGCCCAGTGCCAGCGCACCCGCCGCGAGCAGACGGCGGGACCATTTTTTCATTTAACGCCCCCTGTGTCCCCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTaa   >  1:270164/1‑143 (MQ=255)
                                                                                                                                              gcgcCCATGCCCACTGCCTGGGCTCCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACTTCCCCTGTGTTACCTTTCGCTTTATCACCAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAat  <  2:675145/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GGCGTACATGGTTTCGTTCGACTCGTAAGTCGAATAGATAACCTTAATACCGGTTTCTTTGGTGAACTGTTCAAGCAGTCCTGGCGGCACGTACTCGGTCCAGTTGTAGAAATACAGCGTGTTGTTGTCATCGGCGTGAGCGGCGCTCATGCCCAGTGCCAGAGCACCCGCCGCGAGCAGGTGGCGTGACCATTTTTTCATTTAACGTCCCCTGTGTTACCTTTCGTTTTATCACGAGCAATAAGCTGGCTGGCAATAACCATCACCAGCGACAGCACCAGTAAT  >  NC_000913/1182629‑1182913

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