Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,839,078 G→A 7.3% R111H (CGT→CAT)  ynjD → putative ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,839,0780GA 7.3% 161.7 / 2.9 55R111H (CGT→CAT) ynjDputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base A (2/2);  total (29/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 = 5.60e-01

TGAGCAACGGATTGACATCCTACCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAGCGTTACTCCTGGATGA  >  NC_000913/1838938‑1839222
                                                                                                                                            |                                                                                                                                                
tGAGCAACGGATTGACATCCTACCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCg                                                                                                                                          <  1:13842/149‑1 (MQ=255)
 gAGCAACGGATTGACATCCTACCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGt                                                                                                                                         >  2:135700/1‑149 (MQ=255)
  aGCAACGGATTGACATCCTACCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTGATGCCGt                                                                                                                                         <  7:26675/148‑1 (MQ=255)
    cAACGGATTGACATCCTACCCACCGCACAGCGTCAGATTGGCAATCAATTTCAGGATTCACTGTTATTTGACAAGATAACTGTCGGGCAAAATTTACTGCTGGCGCTACCGGAGACACTTAAAGGGAATGCCCGACGAAATGCCGTGaa                                                                                                                                      >  1:90415/1‑149 (MQ=255)
     aaCGGATTGACATCCTACCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAat                                                                                                                                     >  6:82526/1‑149 (MQ=255)
          aTTGACATCCTACCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAatgatg                                                                                                                                 >  7:28037/1‑148 (MQ=255)
                aTCCTACCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACtt                                                                                                                            >  3:219236/1‑147 (MQ=255)
                 tCCTACCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAg                                                                                                                         >  5:239480/1‑149 (MQ=255)
                    tACCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGt                                                                                                                      >  4:34994/1‑149 (MQ=255)
                     aCCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGtt                                                                                                                     >  8:29565/1‑149 (MQ=255)
                         aCCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAgg                                                                                                                 >  5:26891/1‑149 (MQ=255)
                           cGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGcc                                                                                                               <  7:29565/149‑1 (MQ=255)
                              acaGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCtt                                                                                                             >  8:137299/1‑148 (MQ=255)
                               caGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGa                                                                                                           <  6:144728/149‑1 (MQ=255)
                                aGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAACGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAg                                                                                                          <  3:212822/149‑1 (MQ=255)
                                       gATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCtt                                                                                                   <  6:26891/149‑1 (MQ=255)
                                       gATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTAGAGCGTTCAGGCCTTGAGGGAGCtt                                                                                                   >  8:45269/1‑149 (MQ=255)
                                         ttGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCGAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTc                                                                                                 >  8:93696/1‑149 (MQ=255)
                                            gCATTCTTTTTCAGGATGCACTTTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTTGCGCTACCGGCGACACTTAATGGGTATGCCCGACGTAATGACGTGACTGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCa                                                                                               <  5:73898/148‑1 (MQ=255)
                                              aTTCTTTTTCAGGATGCACTGTTATTTGCCCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCCGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCa                                                                                            <  1:137148/149‑1 (MQ=255)
                                              aTTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCAt                                                                                              >  3:171291/1‑147 (MQ=255)
                                              aTTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACATAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCAt                                                                                              >  7:132609/1‑147 (MQ=255)
                                              aTTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACATAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCAt                                                                                              >  7:88279/1‑147 (MQ=255)
                                               ttCTTTTTCAGGATGCACTGTTATTTGACGAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTGATAGGGAATGCCCGACGTAATGCCGTGAATGATGCACGTGCGCGTTCAGGGCTTGAGGGAGCTGTCCATCAg                                                                                           >  4:123257/1‑149 (MQ=255)
                                               ttCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAg                                                                                           >  7:213723/1‑149 (MQ=255)
                                               ttCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAg                                                                                           >  8:124297/1‑149 (MQ=255)
                                                        aGGCTGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGc                                                                                   <  7:93634/148‑1 (MQ=255)
                                                         ggATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCa                                                                                                                   >  1:66998/1‑115 (MQ=255)
                                                         ggATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCa                                                                                                                   <  2:66998/115‑1 (MQ=255)
                                                                 tGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCt                                                                         >  7:213331/1‑149 (MQ=255)
                                                                     aTTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTc                                                                                                 >  2:98563/1‑121 (MQ=255)
                                                                     aTTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTc                                                                                                 <  1:98563/121‑1 (MQ=255)
                                                                     aTTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCg                                                                     >  6:200229/1‑149 (MQ=255)
                                                                        tGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACAGTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGt                                                                   <  3:123257/148‑1 (MQ=255)
                                                                           ccAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCagc                                                               <  8:161618/149‑1 (MQ=255)
                                                                            cAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCagc                                                               >  7:126329/1‑148 (MQ=255)
                                                                                      tCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGt                                                     <  3:34994/148‑1 (MQ=255)
                                                                                          gCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGctc                                                <  8:213331/149‑1 (MQ=255)
                                                                                            aaaaTTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTcta                                              <  7:74710/149‑1 (MQ=255)
                                                                                            aaaaTTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTcta                                              >  1:10056/1‑149 (MQ=255)
                                                                                                            gcTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGACCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCg                              >  3:109113/1‑149 (MQ=255)
                                                                                                             cTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGc                             <  3:224590/149‑1 (MQ=255)
                                                                                                             cTACCGGCGACACTTAAAGGGAATGCCCGACATAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGc                             <  8:88279/149‑1 (MQ=255)
                                                                                                             cTACCGGCGACACTTAAAGGGAATGCCCGACATAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGc                             <  8:132609/149‑1 (MQ=255)
                                                                                                                  ggCGACACTTAAAGGGAATGCCGGAGGTAATGCCGTGAATGTTGCAGTTGAGCGTTCATGCCTTGAGGGAGCTTTCCATCAGGATCCTGCGACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCccaac                        <  2:90415/149‑1 (MQ=255)
                                                                                                                        acTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAg                  <  4:171291/149‑1 (MQ=255)
                                                                                                                            aaaGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAGCGt               <  1:135700/148‑1 (MQ=255)
                                                                                                                             aaGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAGCGTTa             <  7:124297/149‑1 (MQ=255)
                                                                                                                              aGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAGCGTTAc            <  4:219236/149‑1 (MQ=255)
                                                                                                                               gggAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAGCGTTACt           <  2:10056/149‑1 (MQ=255)
                                                                                                                                  aaTGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAGCGTTACTCCt        >  7:148288/1‑149 (MQ=255)
                                                                                                                                  aaTGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGCCGTCAGCGCTCCCGCGGTTCGCAACTACCCCGCCATAACCCCCCAACAAAAGCCTTACTccat       >  7:211163/1‑147 (MQ=255)
                                                                                                                                      cccGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAGCGTTACTCCTGGAt    >  3:187306/1‑149 (MQ=255)
                                                                                                                                       ccGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAGCGTTACTCCTGGAt    >  3:189265/1‑148 (MQ=255)
                                                                                                                                        cGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAGCGTTACTCCTGGATGa  <  8:213723/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                
TGAGCAACGGATTGACATCCTACCCACCGCACAGCGTCAGATTGGCATTCTTTTTCAGGATGCACTGTTATTTGACCAGTTCAGTGTCGGGCAAAATTTACTGCTGGCGCTACCGGCGACACTTAAAGGGAATGCCCGACGTAATGCCGTGAATGATGCACTTGAGCGTTCAGGCCTTGAGGGAGCTTTCCATCAGGATCCTGCCACTTTGTCTGGCGGTCAGCGAGCGCGCGTTGCTCTACTACGCGCCCTTCTCGCCCAACCAAAAGCGTTACTCCTGGATGA  >  NC_000913/1838938‑1839222

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