Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,005,996 T→A 7.1% I429L (ATA→TTA)  U712_04925 ← Na(+)/H(+) antiporter NhaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,005,9960TA 7.1% 173.5 / 4.6 56I429L (ATA→TTA) U712_04925Na(+)/H(+) antiporter NhaC
Reads supporting (aligned to +/- strand):  ref base T (26/26);  new base A (2/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

GCTGCTAGCATTTCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATGTTTCGTTTGTCATATAAAGATT  >  CP006881/1005859‑1006133
                                                                                                                                         |                                                                                                                                         
gctgctAGCATTTCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTAt                                                                                                                                           >  6:324496/1‑138 (MQ=255)
       gCATTTCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTTTGACTgg                                                                                                                                     <  7:80349/137‑1 (MQ=255)
       gCATTTCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTgg                                                                                                                                     <  7:321748/137‑1 (MQ=255)
       gCATTTCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTgg                                                                                                                                     <  2:293635/137‑1 (MQ=255)
       gCATTTAATTAATGAACGTAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTgg                                                                                                                                     <  2:162132/137‑1 (MQ=255)
          tttCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACt                                                                                                                                       <  8:156110/132‑1 (MQ=255)
          tttCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACt                                                                                                                                       >  7:156110/1‑132 (MQ=255)
          tttCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAacac                                                                                                                                >  7:30664/1‑139 (MQ=255)
          tttCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAAACCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTATATGACTGGAaccc                                                                                                                                >  4:57080/1‑137 (MQ=255)
            tCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACAcc                                                                                                                               >  4:100565/1‑138 (MQ=255)
              aTAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTAAGACTGGAACACCAAg                                                                                                                            >  5:66887/1‑139 (MQ=255)
              aTAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTAAGACTGGAACACCAAg                                                                                                                            >  5:61315/1‑139 (MQ=255)
                  aTGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTgcg                                                                                                                        <  1:271007/139‑1 (MQ=255)
                     aaCGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCt                                                                                                                      <  7:121555/138‑1 (MQ=255)
                     aaCGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCt                                                                                                                      <  5:53280/138‑1 (MQ=255)
                     aaCGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCt                                                                                                                      <  5:360850/138‑1 (MQ=255)
                     aaCGAAGAAAAGAGCACGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTATTCTATGACTGGAACACCAAGTGCGCt                                                                                                                      <  8:369075/138‑1 (MQ=255)
                       cGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGc                                                                                                                   >  6:151902/1‑139 (MQ=255)
                            aaaaGAGCCCGCTTTTTACAGCGGGGCTCTTTTTTATTTTTTACAAATACAAAAAAAAAGGAGATAAATGGCGAATAATAAAGGAATGACGTAAATGGGATTTTGTCTATGCCTGGAAAAACAAATGAGCTCCACATa                                                                                                               >  2:271007/1‑138 (MQ=255)
                                   cccGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGcgcc                                                                                                       >  4:364257/1‑139 (MQ=255)
                                          tttACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGa                                                                                                  <  8:36845/137‑1 (MQ=255)
                                               aGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCa                                                                                           <  5:324496/139‑1 (MQ=255)
                                                           tttATTTTTTAACAAATCCAATCAAAATGGATATCATTGGCGATATATAGAGTAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAAt                                                                                 <  8:173635/137‑1 (MQ=255)
                                                             tATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTAAGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATggg                                                                              <  6:66887/138‑1 (MQ=255)
                                                             tATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTAAGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATggg                                                                              <  6:61315/138‑1 (MQ=255)
                                                                        aaaTCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGtt                                                                  <  3:100565/139‑1 (MQ=255)
                                                                           tCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTATCATTTAATGTTATGCTTACGAACCAAGGAATCAATGGGTTAATTAATGttact                                                               <  8:2637/139‑4 (MQ=255)
                                                                                tCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTc                                                          <  1:205503/139‑1 (MQ=255)
                                                                                tCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTc                                                          <  5:214962/139‑1 (MQ=255)
                                                                                tCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCGCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTc                                                          <  8:107363/139‑1 (MQ=255)
                                                                                tCAAAATGCAGATCATTGGAGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATCCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTc                                                          <  6:144445/139‑1 (MQ=255)
                                                                                                tGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGt                                           >  5:358258/1‑138 (MQ=255)
                                                                                                             agGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATgg                                                                               <  4:316228/89‑1 (MQ=255)
                                                                                                             agGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATgg                                                                               >  3:316228/1‑89 (MQ=255)
                                                                                                             agGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCt                             >  2:370111/1‑139 (MQ=255)
                                                                                                             agGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCt                             >  1:186436/1‑139 (MQ=255)
                                                                                                             agGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCt                             >  7:46470/1‑139 (MQ=255)
                                                                                                             agGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCt                             >  4:340150/1‑139 (MQ=255)
                                                                                                               gaagaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCt                                                                                                                      <  7:280267/48‑1 (MQ=255)
                                                                                                               gaagaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCt                                                                                                                      >  8:280267/1‑48 (MQ=255)
                                                                                                               gaagaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTg                            >  8:110002/1‑138 (MQ=255)
                                                                                                                aagaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCt                                                                                                                      <  3:166572/47‑1 (MQ=255)
                                                                                                                aagaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCt                                                                                                                      >  4:166572/1‑47 (MQ=255)
                                                                                                                  gaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGCCGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGCATg                        >  2:244787/1‑139 (MQ=255)
                                                                                                                  gaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATggg                                                                              >  2:116662/1‑85 (MQ=255)
                                                                                                                  gaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATggg                                                                              <  1:116662/85‑1 (MQ=255)
                                                                                                                  gaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATg                        >  2:236789/1‑139 (MQ=255)
                                                                                                                  gaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATg                        >  8:251784/1‑139 (MQ=255)
                                                                                                                  gaaCGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATg                        >  3:313556/1‑139 (MQ=255)
                                                                                                                         aaTGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATgg                                                                               <  7:206672/77‑1 (MQ=255)
                                                                                                                         aaTGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATgg                                                                               >  8:206672/1‑77 (MQ=255)
                                                                                                                                tGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATGTTTCGTTTGTCata          >  2:333823/1‑139 (MQ=255)
                                                                                                                                    gTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATGTTTCGTTTGTCATATaa       <  7:251784/138‑1 (MQ=255)
                                                                                                                                       tATGACTGGAGCACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATGTTTCGTTTGTCATATAAAGaa   >  5:239989/1‑138 (MQ=255)
                                                                                                                                       tATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGCCGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATGTTTCGTTTGTCATATAAAGAt   >  6:256740/1‑139 (MQ=255)
                                                                                                                                       tATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCa                                                                                           >  7:220940/1‑51 (MQ=255)
                                                                                                                                       tATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCa                                                                                           <  8:220940/51‑1 (MQ=255)
                                                                                                                                         tGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATGTTTCGTTTGTCATATAAAGAtt  >  1:19358/1‑138 (MQ=255)
                                                                                                                                         tGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATGTTTCGTTTGTCATATAAAGAt   <  4:313556/137‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                         
GCTGCTAGCATTTCATAAATGAACGAAGAAAAGAGCCCGCTTTTTACAGCGGGGCTCTTTTTATTTTTTAACAAATCCAATCAAAATGGAGATCATTGGCGATATATAGAGGAAGAACGTAAATGGGATGTAGTCTATGACTGGAACACCAAGTGCGCTCGCCATAAATGCGCCGCTGACGCCCCAAGGAATCAATGGGTTAATTAATGTTCCGCCGTCCTCTAATGATCTTGTCAGAAATTTTCGCTGGATGTTTCGTTTGTCATATAAAGATT  >  CP006881/1005859‑1006133

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑