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

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: ‑