Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,772,240 A→G 6.7% K390K (AAA→AAG ydiM → putative MFS transporter, membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,772,2400AG 6.7% 187.2 / 3.4 60K390K (AAA→AAGydiMputative MFS transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base A (25/31);  new base G (2/2);  total (27/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.42e-01

GCGGGCAGTATTGCGACCTTTACTATTCCGTTGATTACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATCGCGATCAAAAACAAACTATGACATGCA  >  NC_000913/1772094‑1772383
                                                                                                                                                  |                                                                                                                                               
gCGGGCAGTATTGCGACCTTTACTATTCCGTTGATTACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAAc                                                                                                                                               <  1:463551/149‑1 (MQ=255)
      aGTATTGCGACCTTTACTATTCCGTTGATTACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCa                                                                                                                                         >  2:334008/1‑149 (MQ=255)
                  tttACTATTCCGTTGATTACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTaaa                                                                                                                             >  2:89887/1‑149 (MQ=255)
                     aCTATTCCGTTGATTACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGa                                                                                                                          <  1:172042/149‑1 (MQ=255)
                      cTATTCCGTTGATTACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGaa                                                                                                                         <  2:531/149‑1 (MQ=255)
                       tATTCCGTTGATTACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGaaa                                                                                                                        >  3:86064/1‑149 (MQ=255)
                       tATTCCGTTGATTACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGaaa                                                                                                                        >  2:307566/1‑149 (MQ=255)
                        aTTCCGTTGATTACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGaaaa                                                                                                                       <  2:357554/149‑1 (MQ=255)
                                  ttACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCg                                                                                                             >  4:118053/1‑149 (MQ=255)
                                   tACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCgg                                                                                                            <  3:20021/149‑1 (MQ=255)
                                      ggCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCgggtg                                                                                                         <  2:340504/149‑1 (MQ=255)
                                        cTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCgggtggg                                                                                                       <  2:363785/149‑1 (MQ=255)
                                            tCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAAt                                                                                                   >  4:194911/1‑149 (MQ=255)
                                              tGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATGACTCGCTAAAGGAAAATGTCGCTCCgggt                                                                                                          <  4:148954/140‑1 (MQ=255)
                                              tGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATGACTCGCTAAAGGAAAATGTCGCTCCgggt                                                                                                          >  3:148954/1‑140 (MQ=255)
                                                  ccAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCTGGGTGGGTAATGCAata                                                                                             >  4:493754/1‑149 (MQ=255)
                                                  ccAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAata                                                                                             >  1:549401/1‑149 (MQ=255)
                                                   cAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAatat                                                                                            >  3:93053/1‑149 (MQ=255)
                                                          gTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTc                                                                                     <  4:364591/149‑1 (MQ=255)
                                                            aTTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAg                                                                                   <  3:415259/149‑1 (MQ=255)
                                                            aTTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAg                                                                                   <  2:549401/149‑1 (MQ=255)
                                                               gCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTa                                                                                                                               >  2:12900/1‑102 (MQ=255)
                                                               gCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTa                                                                                                                               <  1:12900/102‑1 (MQ=255)
                                                                ccGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGaaaa                                                                                                                       <  4:165649/109‑1 (MQ=255)
                                                                ccGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGaaaa                                                                                                                       >  3:165649/1‑109 (MQ=255)
                                                                 cGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAGACGCGGCAt                                                                                                                                        <  1:392638/91‑1 (MQ=255)
                                                                 cGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAGACGCGGCAt                                                                                                                                        >  2:392638/1‑91 (MQ=255)
                                                                    tagtatGTGGTTCGATACTGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTGAAGCAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCa                                                                           >  4:392155/4‑149 (MQ=255)
                                                                        atGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTGAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGAt                                                                       <  3:392155/149‑1 (MQ=255)
                                                                        atGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGAt                                                                       <  1:334008/149‑1 (MQ=255)
                                                                        atGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGAt                                                                       <  1:593727/149‑1 (MQ=255)
                                                                          gtgGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCt                                                                     >  1:78685/1‑149 (MQ=255)
                                                                            ggTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAg                                                                                                                            <  4:27200/92‑1 (MQ=255)
                                                                            ggTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAg                                                                                                                            >  3:27200/1‑92 (MQ=255)
                                                                             gTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTAc                                                                  <  1:75937/149‑1 (MQ=255)
                                                                              ttCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACg                                                                 <  2:40492/149‑1 (MQ=255)
                                                                                cGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGa                                                               >  1:446670/1‑149 (MQ=255)
                                                                                       gCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGataaataa                                                        <  3:27942/149‑1 (MQ=255)
                                                                                               tGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCt                                                <  1:89887/149‑1 (MQ=255)
                                                                                                 ccATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTgc                                              >  3:420177/1‑149 (MQ=255)
                                                                                                   aTCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTgcgc                                            >  4:519548/1‑149 (MQ=255)
                                                                                                    tCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGct                                           >  2:82734/1‑149 (MQ=255)
                                                                                                       gTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGctaac                                        >  1:336808/1‑149 (MQ=255)
                                                                                                          ttttACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGctaactaa                                     <  4:93053/149‑1 (MQ=255)
                                                                                                          ttttACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGctaactaa                                     <  3:102867/149‑1 (MQ=255)
                                                                                                           tttACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAg                                    <  4:420177/149‑1 (MQ=255)
                                                                                                            ttACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAgg                                   >  2:264187/1‑149 (MQ=255)
                                                                                                                  gCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATc                             <  4:484385/149‑1 (MQ=255)
                                                                                                                      tGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATCGCGa                         >  1:350785/1‑149 (MQ=255)
                                                                                                                        tttATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATCGCGATc                       >  1:179131/1‑149 (MQ=255)
                                                                                                                            tCGGCTTACGCAGCCGCAAAAAGACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGa                                                               <  2:183238/105‑1 (MQ=255)
                                                                                                                            tCGGCTTACGCAGCCGCAAAAAGACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGa                                                               >  1:183238/1‑105 (MQ=255)
                                                                                                                            tCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATCGCGATCaaaa                   <  1:264187/149‑1 (MQ=255)
                                                                                                                               gCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATCGCGATCaaaaa                  >  2:329169/1‑147 (MQ=255)
                                                                                                                               gCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATCGCGATCaaaaa                  <  1:329169/147‑1 (MQ=255)
                                                                                                                                cTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATCGCGATCAAaaacaa               <  2:479745/149‑1 (MQ=255)
                                                                                                                                 ttACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTCTTCACGTCATGCAAGATCTTACGGATAAATACCTCTTTCTGCGCTAACTAAGGAAAATCGCGATCAAaaccaaa              <  2:480576/149‑1 (MQ=255)
                                                                                                                                            cAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAat                                                                                              >  4:118268/1‑58 (MQ=255)
                                                                                                                                            cAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAat                                                                                              <  3:118268/58‑1 (MQ=255)
                                                                                                                                            cAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATCGCGATCAAAAACAAACTATGACATGc   <  4:516375/149‑1 (MQ=255)
                                                                                                                                             aaaaaaaCGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATCGCGATCAAAAACAAACTATGACATGCa  >  1:304374/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCGGGCAGTATTGCGACCTTTACTATTCCGTTGATTACGGCTCATCTGTCCCAAAGAAGTATTGCCGATATTATGTGGTTCGATACCGCCATCGCTGCCATCGGTTTTTTACTGGCACTGTTTATCGGCTTACGCAGCCGCAAAAAAACGCGGCATCACTCGCTAAAGGAAAATGTCGCTCCGGGTGGGTAATGCAATATTCTTTTCAGGTCATGCAAGATCTTACGGATAAATAACTCTTTCTGCGCTAACTAAGGAAAATCGCGATCAAAAACAAACTATGACATGCA  >  NC_000913/1772094‑1772383

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: