Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 17,527 G→C 7.3% S13S (TCG→TCC)  nhaA → sodium‑proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,5270GC 7.3% 238.3 / 8.2 68S13S (TCG→TCC) nhaAsodium‑proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base C (3/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATG  >  NC_000913/17403‑17673
                                                                                                                            |                                                                                                                                                  
ttCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGc                                                                                                                             >  2:651901/1‑148 (MQ=255)
       acTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCg                                                                                                                                                    <  4:940787/118‑1 (MQ=255)
       acTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCg                                                                                                                                                    >  3:940787/1‑118 (MQ=255)
          ataATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTAt                                                                                                                                   <  2:540295/132‑1 (MQ=255)
          ataATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTAt                                                                                                                                   >  1:540295/1‑132 (MQ=255)
                 gATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCg                                                                                                                >  4:298120/1‑144 (MQ=255)
                 gATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCg                                                                                                                <  3:298120/144‑1 (MQ=255)
                                 cGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAgcgg                                                                                            <  1:651901/148‑1 (MQ=255)
                                            aaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCg                                                                                                                <  1:1033772/117‑1 (MQ=255)
                                            aaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCg                                                                                                                >  2:1033772/1‑117 (MQ=255)
                                                            cTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAg                                                                                               <  1:186144/118‑1 (MQ=255)
                                                            cTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAg                                                                                               >  2:186144/1‑118 (MQ=255)
                                                              aTTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTAt                                                                                                                                   >  2:961504/1‑80 (MQ=255)
                                                              aTTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTAt                                                                                                                                   <  1:961504/80‑1 (MQ=255)
                                                               ttCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGc                                                                                                                                              <  4:188836/68‑1 (MQ=255)
                                                               ttCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGc                                                                                                                                              >  3:188836/1‑68 (MQ=255)
                                                                    cTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCattat                                                                                                                                         >  1:1070465/1‑68 (MQ=255)
                                                                    cTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCattat                                                                                                                                         <  2:1070465/68‑1 (MQ=255)
                                                                    cTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAg                                                                                   >  2:821678/1‑122 (MQ=255)
                                                                    cTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAg                                                                                   <  1:821678/122‑1 (MQ=255)
                                                                      gAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGc                                                                                                                                              >  1:775413/1‑61 (MQ=255)
                                                                      gAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGc                                                                                                                                              <  2:775413/61‑1 (MQ=255)
                                                                      gAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCatt                                                                                                                                           >  4:1105444/1‑64 (MQ=255)
                                                                      gAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCatt                                                                                                                                           <  3:1105444/64‑1 (MQ=255)
                                                                               aTAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGa                                                                               >  2:489355/1‑115 (MQ=255)
                                                                               aTAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGa                                                                               <  1:489355/115‑1 (MQ=255)
                                                                               aTAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCAGTATTCTTATCATTGCCGCTATCCTg                                                                                                                   >  2:910381/1‑79 (MQ=255)
                                                                               aTAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGAGAGGCATTATTCTTATCATTGCCGCTATCCTg                                                                                                                   <  1:910381/79‑1 (MQ=255)
                                                                                 aaaaaGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTAt                                                                                                        <  3:849931/88‑1 (MQ=255)
                                                                                 aaaaaGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTAt                                                                                                        >  4:849931/1‑88 (MQ=255)
                                                                                   aaaGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATgg                                                                                                      <  4:66660/88‑1 (MQ=255)
                                                                                   aaaGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATgg                                                                                                      >  3:66660/1‑88 (MQ=255)
                                                                                   aaaGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGt                                                                           >  1:1097237/1‑115 (MQ=255)
                                                                                   aaaGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGt                                                                           <  2:1097237/115‑1 (MQ=255)
                                                                                    aaGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgat                                                                                                              >  2:694298/1‑79 (MQ=255)
                                                                                    aaGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgat                                                                                                              <  1:694298/79‑1 (MQ=255)
                                                                                        gAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGa                                                                    <  3:583032/117‑1 (MQ=255)
                                                                                        gAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGa                                                                    >  4:583032/1‑117 (MQ=255)
                                                                                            cATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAgcgg                                                                                            >  4:18722/1‑89 (MQ=255)
                                                                                            cATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAgcgg                                                                                            <  3:18722/89‑1 (MQ=255)
                                                                                            cATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGggttggtt                                 >  3:491105/1‑148 (MQ=255)
                                                                                              tCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgatgat                                                                                                           <  4:1107580/72‑1 (MQ=255)
                                                                                              tCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgatgat                                                                                                           >  3:1107580/1‑72 (MQ=255)
                                                                                               cTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCAc                              <  2:512128/148‑1 (MQ=255)
                                                                                               cTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCAc                              >  1:512128/1‑148 (MQ=255)
                                                                                                   aTCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTgg                                                            >  4:258900/1‑114 (MQ=255)
                                                                                                   aTCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTgg                                                            <  3:258900/114‑1 (MQ=255)
                                                                                                    tCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTgg                                                            <  3:903914/113‑1 (MQ=255)
                                                                                                    tCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTgg                                                            >  4:903914/1‑113 (MQ=255)
                                                                                                     cGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAg                                              <  1:351299/126‑1 (MQ=255)
                                                                                                     cGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAg                                              >  2:351299/1‑126 (MQ=255)
                                                                                                       aTTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTa                                                                                                         >  1:1074065/1‑65 (MQ=255)
                                                                                                       aTTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTa                                                                                                         <  2:1074065/65‑1 (MQ=255)
                                                                                                        ttCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATca                     >  3:1095097/1‑148 (MQ=255)
                                                                                                            ttAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACg                                                                     >  1:71904/1‑96 (MQ=255)
                                                                                                            ttAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACg                                                                     <  2:71904/96‑1 (MQ=255)
                                                                                                            ttAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACg                                                       >  1:578008/1‑110 (MQ=255)
                                                                                                            ttAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACg                                                       <  2:578008/110‑1 (MQ=255)
                                                                                                            ttAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaacaa                 >  1:554064/1‑148 (MQ=255)
                                                                                                                 aataaTGCCTCCCGAGGCATTATTCTTGTCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTgg                                                            >  3:978570/5‑100 (MQ=255)
                                                                                                                  gTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACAt           <  2:554064/148‑1 (MQ=255)
                                                                                                                  ataaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGaaa                        >  2:584185/4‑135 (MQ=255)
                                                                                                                  ataaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGaaa                        <  1:584185/132‑1 (MQ=255)
                                                                                                                   taaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATg                                                                             >  3:429691/3‑81 (MQ=255)
                                                                                                                   taaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATg                                                                             <  4:429691/79‑1 (MQ=255)
                                                                                                                       gCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGt      >  4:926493/1‑148 (MQ=255)
                                                                                                                          tCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCa                               >  2:1010324/1‑120 (MQ=255)
                                                                                                                          tCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCa                               <  1:1010324/120‑1 (MQ=255)
                                                                                                                          tCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTAt   <  4:491105/148‑1 (MQ=255)
                                                                                                                           cGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATg  <  4:1095097/148‑1 (MQ=255)
                                                                                                                            |                                                                                                                                                  
TTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATG  >  NC_000913/17403‑17673

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