Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,195,626 A→G intergenic (+30/+464) icdA → / ← ymfD isocitrate dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,195,6260AG97.4% 230.0 / ‑2.1 78intergenic (+30/+464)icdA/ymfDisocitrate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (33/43);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

GGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGT  >  NC_000913/1195536‑1195718
                                                                                          |                                                                                            
ggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGc                                                                                            <  1:1323655/93‑1 (MQ=255)
 gatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCg                                                                                           <  2:2963422/93‑1 (MQ=255)
 gatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCg                                                                                           >  2:3213159/1‑93 (MQ=255)
  atggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGt                                                                                          <  2:1795086/93‑1 (MQ=255)
  atggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGt                                                                                          <  1:2320224/93‑1 (MQ=255)
   tggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTa                                                                                         >  1:3929200/1‑93 (MQ=255)
       gcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGc                                                                                     >  1:3535823/1‑93 (MQ=255)
        cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCt                                                                                    >  2:1863034/1‑93 (MQ=255)
         tAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTc                                                                                   >  1:2355588/1‑93 (MQ=255)
         tAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTc                                                                                   >  1:4316761/1‑93 (MQ=255)
          aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTcc                                                                                  >  1:542236/1‑93 (MQ=255)
          aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTcc                                                                                  <  1:2262929/93‑1 (MQ=255)
           aaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTccc                                                                                 <  2:2051382/93‑1 (MQ=255)
            aCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCg                                                                                >  2:3349552/1‑93 (MQ=255)
             ctgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCgt                                                                               >  2:1850189/1‑93 (MQ=255)
             ctgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCgt                                                                               >  1:2969722/1‑93 (MQ=255)
              tgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCgtt                                                                              <  2:2347628/93‑1 (MQ=255)
                  gAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGttt                                                                          <  2:2272752/93‑1 (MQ=255)
                   aaaTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGtttt                                                                         <  1:3178790/93‑1 (MQ=255)
                    aaTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGttttt                                                                        >  2:2335005/1‑93 (MQ=255)
                      tGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTg                                                                      >  1:1236494/1‑93 (MQ=255)
                      tGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTg                                                                      >  2:4418174/1‑93 (MQ=255)
                         tCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTa                                                                   <  2:4256356/93‑1 (MQ=255)
                         tCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTa                                                                   <  2:2312373/93‑1 (MQ=255)
                           agagTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAgg                                                                 >  2:3457139/1‑93 (MQ=255)
                            gagTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGgc                                                                <  1:2217777/93‑1 (MQ=255)
                            gagTTTGGTGACGCGATCAACGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGgc                                                                >  2:4564027/1‑93 (MQ=255)
                             agTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGgct                                                               >  1:1452058/1‑93 (MQ=255)
                              gTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGgctg                                                              >  2:2484188/1‑93 (MQ=255)
                                  ggTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTaa                                                          <  1:1775309/93‑1 (MQ=255)
                                   gTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAAc                                                         <  1:3761296/93‑1 (MQ=255)
                                   gTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAAc                                                         <  1:3918994/93‑1 (MQ=255)
                                   gTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAAc                                                         <  1:4263628/93‑1 (MQ=255)
                                      aCGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGt                                                      <  2:1815935/93‑1 (MQ=255)
                                      aCGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGt                                                      >  1:3083173/1‑93 (MQ=255)
                                      aCGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGt                                                      <  1:1951712/93‑1 (MQ=255)
                                          gATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATc                                                  >  1:3171298/1‑93 (MQ=255)
                                           atcatcGACAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCa                                                 >  2:4568208/1‑93 (MQ=255)
                                            tcatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCaa                                                <  2:633624/93‑1 (MQ=255)
                                            tcatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCaa                                                <  2:2707779/93‑1 (MQ=255)
                                               tcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAt                                             <  2:2276322/93‑1 (MQ=255)
                                               tcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAt                                             >  1:4588169/1‑93 (MQ=255)
                                                 gAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAttt                                           >  1:1293454/1‑93 (MQ=255)
                                                  aaaaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAtttt                                          >  2:216874/1‑93 (MQ=255)
                                                   aaaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTa                                         <  1:910914/93‑1 (MQ=255)
                                                   aaaCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTa                                         <  2:2124297/93‑1 (MQ=255)
                                                     aCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATc                                       <  2:3995269/93‑1 (MQ=255)
                                                     aCATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATc                                       >  2:2378782/1‑93 (MQ=255)
                                                       aTGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaa                                     >  2:3869411/1‑93 (MQ=255)
                                                       aTGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaa                                     <  1:3755745/93‑1 (MQ=255)
                                                        tGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaa                                    <  1:2164692/93‑1 (MQ=255)
                                                        tGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaa                                    <  2:1987217/93‑1 (MQ=255)
                                                        tGTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaa                                    >  1:4110110/1‑93 (MQ=255)
                                                         gTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaaa                                   <  2:4227012/93‑1 (MQ=255)
                                                         gTAATGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaaa                                   >  2:2384519/1‑93 (MQ=255)
                                                          tAATGCCGTAGTTGGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaaaa                                  <  1:988271/93‑1 (MQ=255)
                                                            aTGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaaaaaa                                <  2:2524335/93‑1 (MQ=255)
                                                             tGCCGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAg                               <  1:910423/93‑1 (MQ=255)
                                                              ggcGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGt                              <  1:1547525/91‑1 (MQ=255)
                                                                cGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTa                            <  1:2183783/93‑1 (MQ=255)
                                                                cGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTa                            >  2:2409661/1‑93 (MQ=255)
                                                                cGTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTa                            <  1:3452440/93‑1 (MQ=255)
                                                                 gTAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTAt                           >  1:2521315/1‑93 (MQ=255)
                                                                  tAGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATc                          <  1:730607/93‑1 (MQ=255)
                                                                   aGTTTGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCa                         <  2:4519/93‑1 (MQ=255)
                                                                     tttGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCaaa                       <  1:946844/93‑1 (MQ=255)
                                                                       tGTTAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAc                     <  1:4419809/93‑1 (MQ=255)
                                                                          tAAATTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAcccc                  >  2:2017735/1‑93 (MQ=255)
                                                                           aaaTTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCt                 <  2:4014459/93‑1 (MQ=255)
                                                                           aaaTTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCt                 <  1:2476563/93‑1 (MQ=255)
                                                                            aaTTTATTAACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTc                >  1:3278379/1‑93 (MQ=255)
                                                                                tattaACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTa            >  1:4660191/1‑93 (MQ=255)
                                                                                tattaACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTa            >  1:3943698/1‑93 (MQ=255)
                                                                                  ttaACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAACAAGTTATCAAAACCCCTCGGTAGt          >  1:4426506/1‑93 (MQ=255)
                                                                                  ttaACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGt          >  1:1429704/1‑93 (MQ=255)
                                                                                   taACGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGtt         >  1:2056676/1‑93 (MQ=255)
                                                                                    aaCGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGttt        <  2:2164693/93‑1 (MQ=255)
                                                                                     aCGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGtttt       <  2:3886210/93‑1 (MQ=255)
                                                                                     aCGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGtttt       <  1:3684045/93‑1 (MQ=255)
                                                                                      cGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTg      <  2:2932262/93‑1 (MQ=255)
                                                                                      cGGGGGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTg      <  2:2786894/93‑1 (MQ=255)
                                                                                       ggggGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgg     >  2:1181612/1‑93 (MQ=255)
                                                                                       ggggGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgg     <  1:1013907/93‑1 (MQ=255)
                                                                                          aGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGt  >  1:4412521/1‑93 (MQ=255)
                                                                                          |                                                                                            
GGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGT  >  NC_000913/1195536‑1195718

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑