Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,088,095 T→A intergenic (+235/‑189) metK → / → galP S‑adenosylmethionine synthetase/D‑galactose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,088,0950TA96.5% 173.6 / NA 58intergenic (+235/‑189)metK/galPS‑adenosylmethionine synthetase/D‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (21/35);  minor base C (2/0);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACT  >  NC_000913/3087955‑3088237
                                                                                                                                            |                                                                                                                                              
aaTAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGaaa                                                                                                                                              >  2:669724/1‑143 (MQ=255)
aaTAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGaaa                                                                                                                                              >  1:858362/1‑143 (MQ=255)
 aTAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGaaaa                                                                                                                                             <  1:137770/143‑1 (MQ=255)
   aaCATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGt                                                                                                                                           <  2:635032/143‑1 (MQ=255)
    aCATCAATCTCCATGATCACGGTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTc                                                                                                                                          <  1:349743/143‑1 (MQ=255)
      atcatTCTTCCTGATCACGTTTCACCGCAGAGTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTGGTTACATTTCTTTTCAGAAAAGTCtt                                                                                                                                        >  1:753923/1‑143 (MQ=255)
         atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCTACCACTACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAAt                                                                                                                                     <  2:573494/143‑1 (MQ=255)
         atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAAt                                                                                                                                     >  1:464329/1‑143 (MQ=255)
         atTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAAACAAAAATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAACGTCTTAAt                                                                                                                                     >  2:64087/1‑143 (MQ=255)
             ttcCTGACCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAATATTTGTAATAGTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCa                                                                                                                                 <  1:67702/143‑1 (MQ=255)
             ttcCGGATCACTTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCa                                                                                                                                 <  2:739376/143‑1 (MQ=255)
              gccTGATCACGTTTCACCGCAGATTATCATCACACTTGCAACCTATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAg                                                                                                                                <  2:674028/142‑1 (MQ=255)
                 tGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATa                                                                                                                             >  1:754963/1‑143 (MQ=255)
                 tGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATACTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATa                                                                                                                             >  2:661848/1‑143 (MQ=255)
                 tGATCACGTTTCACCGCAAATTATCATCCCAACTGAAACCGATTACACCAACCACAACAGCCAAAGCTTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCCGAAAAGTCTTCATTGCAAATa                                                                                                                             >  2:20750/1‑143 (MQ=255)
                  gATCACGCTTCACCGCAGAGGAGCAGCACAACTGGAACCGAGTACAGCAACCACAACAGACAAAGATGTGTAATATTTTCATAGTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATaa                                                                                                                            <  1:668890/143‑1 (MQ=255)
                        gTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGCAAAGTCTTAATTGCAGATAACAGCGt                                                                                                                      >  1:711967/1‑143 (MQ=255)
                        gTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAGGTCTTAATTGCAGATAACAGCGt                                                                                                                      >  1:771949/1‑143 (MQ=255)
                               cGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCCGAAAAGTCTTAATTGCAGATAACACCGTCTAATcc                                                                                                               >  2:859244/1‑142 (MQ=255)
                                  aGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCGGATACCAGCGTTTAACCTatg                                                                                                            >  2:453888/1‑143 (MQ=255)
                                   gATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTatga                                                                                                           <  2:271907/143‑1 (MQ=255)
                                   aaTTATCCTCACGACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTatga                                                                                                           <  2:845847/142‑1 (MQ=255)
                                     ttATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTatgatg                                                                                                         >  1:612874/1‑143 (MQ=255)
                                       atcatcACAACTGAAACCGATTACACCAACCCCAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATgatgat                                                                                                       <  1:382986/143‑1 (MQ=255)
                                           tcaCAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACCGCGTTTATTCTATGCAGataaca                                                                                                   >  2:452340/1‑140 (MQ=255)
                                                aCTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCaa                                                                                              <  1:49684/143‑1 (MQ=255)
                                                 cTGAAACCGATTACACCACCCACGACAGACAAAGATGTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAat                                                                                             <  1:348155/143‑1 (MQ=255)
                                                 cTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAat                                                                                             <  1:823431/143‑1 (MQ=255)
                                                    aaaCCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAattat                                                                                          >  2:645936/1‑143 (MQ=255)
                                                     aaCCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTGCATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAattatt                                                                                         <  1:612041/143‑1 (MQ=255)
                                                        cGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCAAATTATTATTCGGTTTTCACAGTTGTTACAAAAAAAAAAAAAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTc                                                                                      <  2:183347/143‑1 (MQ=255)
                                                         gATTACACCAACCACAACAGACGAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCa                                                                                     <  1:547701/143‑1 (MQ=255)
                                                            tACGCCAACGGCCACAGACAAAGATCGGTAATCTTTTCATATTATTATTCGGTTTGCACAGTTGTTACAGTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGc                                                                                  <  1:814204/143‑1 (MQ=255)
                                                              cacCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCAc                                                                                <  1:196777/143‑1 (MQ=255)
                                                                      acaacaGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCa                                                                        <  1:824018/143‑1 (MQ=255)
                                                                       caacaGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCAt                                                                       <  1:614586/143‑1 (MQ=255)
                                                                       caacaGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATCTCTTTTCAGAAAAGTCTTAATTGCTGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCAt                                                                       <  2:474796/143‑1 (MQ=255)
                                                                        aacaGACAAAGATTTGTAATATTTTCATCTGATGATTCGGTTTTCACAGTTGTCCCATCACGGTGCTGCAATGCCCTAAGTGCAGGTAGCAGCGTTTAATCGATGATGATATAACTCAATTATTTTCATGCACTTAAATCATa                                                                      >  2:931557/1‑143 (MQ=255)
                                                                         acaGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATaa                                                                     <  2:171232/143‑1 (MQ=255)
                                                                             actaaGATTTTGAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTaa                                                                 <  2:436539/140‑1 (MQ=255)
                                                                             aCAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTGATTACAGCTAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTaa                                                                 <  1:359881/143‑1 (MQ=255)
                                                                                   aTTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATaaa                                                           <  1:866986/143‑1 (MQ=255)
                                                                                   aTTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATaaa                                                           <  1:674530/143‑1 (MQ=255)
                                                                                    tttGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAAt                                                          <  2:372429/143‑1 (MQ=255)
                                                                                    tttGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAAt                                                          <  2:754963/143‑1 (MQ=255)
                                                                                        tAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTa                                                      <  2:222055/143‑1 (MQ=255)
                                                                                         aaTATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAg                                                     <  1:669724/143‑1 (MQ=255)
                                                                                         aaTATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGAAGATATAACTCAATTATTTTCATTCACTTAAATCATAACTAAGATAAATGTTAg                                                     >  2:333484/1‑143 (MQ=255)
                                                                                              tttCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATACATGTTAGTGTaa                                                >  1:123702/1‑143 (MQ=255)
                                                                                                  atatTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGa                                            <  2:953508/143‑1 (MQ=255)
                                                                                                     ttattattCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTa                                         >  2:343348/1‑143 (MQ=255)
                                                                                                      tattattCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAACTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTac                                        >  2:905760/1‑143 (MQ=255)
                                                                                                                   tCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGAttt                           <  1:93820/143‑1 (MQ=255)
                                                                                                                            ttACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACAt                  >  2:302122/1‑143 (MQ=255)
                                                                                                                              aCATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCt                >  2:164349/1‑143 (MQ=255)
                                                                                                                               cATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGGTGTGTTTTGCTCTActtctt               >  1:279012/1‑143 (MQ=255)
                                                                                                                                aTTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCttt              <  1:86660/143‑1 (MQ=255)
                                                                                                                                 tttCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCtttt             <  1:310373/143‑1 (MQ=255)
                                                                                                                                 tttCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCtttt             <  1:181187/143‑1 (MQ=255)
                                                                                                                                     ttttCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTAcg         >  1:28816/1‑143 (MQ=255)
                                                                                                                                         cAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATGACACTGATGTGATTTGCTTCACATCTTTTTAcgtcgg     >  1:732423/1‑142 (MQ=255)
                                                                                                                                            tAAACTCTTCATTTCAGATAACATCTTTTAATCTATTATTATATAACTCATTTATTTTCATGCACTTAAATCATAATTAATATAAATCTTAGTGTATTCTATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACt  <  2:476380/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
AATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACT  >  NC_000913/3087955‑3088237

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