Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 954,356 C→A 25.5% L37L (CTG→CTT focA ← formate channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913954,3560CA25.5% 69.4 / 2.2 36L37L (CTG→CTTfocAformate channel
Reads supporting (aligned to +/- strand):  ref base C (10/16);  new base A (2/8);  total (12/24)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCC  >  NC_000913/954126‑954548
                                                                                                                                                                                                                                            |                                                                                                                                                                                                  
cAGTTTTTCGCCAACTGACCCCAGGTTATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGa                                                                                                                                                                                             <  2:129106/238‑1 (MQ=255)
cAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGa                                                                                                                                                                                             <  2:520653/238‑1 (MQ=255)
        cGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAAGGCGATTGAGAGGAAAACACCGGCGGGAATCTCCAGATAGaaaaagataa                                                                                                                                                                                    >  1:603136/1‑232 (MQ=255)
                    ctcggTGTTCCGCCCCCTCCCCTTAGAAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCTCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATTCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACGCCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTcc                                                                                                                                                                       <  1:104260/236‑2 (MQ=255)
                    cctgggtagggcctcACTCGCTTTAGAAACAACAAATAACGCGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCCGCCCCGGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGGGCCAGTGGTTGCTGTGATATAGAATACGAATGCGATTGATATGAAATCTCCCGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCg                                                                                                                                                                       <  1:732586/225‑1 (MQ=255)
                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCg                                                                                                                                                                       <  2:689889/239‑1 (MQ=255)
                    ccAGGGGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCg                                                                                                                                                                       <  2:814709/239‑1 (MQ=255)
                    ccAAGTTATGCCCCCTCTCGCCTTAGCAACAACAAGTAACACGGTGGAAGTAAAGAGATCCGCTCCGCAGACAACACAAATAATCAGCCCCATAGAGAAGCAAATTCCGACAACCAGTTTTGCCATGCCGAAGGGCAGTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACCCCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCg                                                                                                                                                                       <  1:151319/239‑1 (MQ=255)
                        gggttGCGCCCCCTCCCCTAAGAAAAAAAAATCAACACGGGGGAAATAAAGAGATCGGCTCCGCAGAAAACCCAAAGAACTAGCCCCAAGAAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGGCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACGCCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGc                                                                                                                                                                   <  1:701664/235‑1 (MQ=255)
                             gcgcCCACTCGCCTTAGCAACAACAATCAACACGGTGGCAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGAGGTTGCTGGGATATAGAAGACGAATGCGATTGAGATGCAAACACCGGCGGTGATCCCCAGAAAGAAACACTTaaagggtggttttggtgtgttta                                                                                                                                                               >  1:164802/1‑218 (MQ=255)
                                                 caacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGATTGCGGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGTCAATTGCCCGAAAGACAGCCTCAAGCGGAGGTTTCACTGTTTTTTAGAACACACCCTCTTCGGc                                                                                                                                          >  2:492339/1‑238 (MQ=255)
                                                    caaTCAACGCGGGGGAAGAAAAGAGATCGGCCCCGAAGACAACAAAAAAAATCAGCCCAAGAGAAAAGCAAATGCCGCCAACCAGTTTTGCCATGCAGAAGGGCATTGTGCCTGTGCCAGGGGTGGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGTTTTATAGACACCCGCCTCTTCGGCCACt                                                                                                                                      <  2:603136/239‑1 (MQ=255)
                                                    caaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCg                                                                                                                                                                       <  2:168199/207‑1 (MQ=255)
                                                       tCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGGAAATCCCAAGTAAGAAGGCTTTAAGGGGATGTTCGTTGGTTTATAGGCACCCCCCCCTTGCGCCACTTTg                                                                                                                                   >  2:354358/1‑239 (MQ=255)
                                                       acaccacgggggagaaaaaGAGGTCGGCCCCGAAAACAACACAAACAACCACCCCAAGGGAGAAGAAAATGCCGACACCCATTTTTGCCAGGCCGAAGCGATTGGGGCCTGTGCCAGTGGTTGGTGTGATAGAAGAGAAAAATGCGATTGAGAAGAAAACCCCGGGGGTAATCGCCAGAAAGAAATTCTTAAGCGGTTGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTg                                                                                                                                   <  1:623523/223‑1 (MQ=255)
                                                               gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGCGATGAAAACACCGGCGGGAAATCCCCAAATAAAAAGTTTTAAGGGGGGATTTATTGGTTTTATATAAACCCGCCCTTTAGGCCCATTTGGcatttc                                                                                                                            >  1:104355/1‑234 (MQ=255)
                                                               gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATAGCCAGATAGAAAGGCTTAAACGGGTGTTTTGTTGGTTTTTAGACACCCCCCCCTTTGGCCCCTTTTGCCATTgc                                                                                                                           >  2:66297/1‑239 (MQ=255)
                                                                   aaGTAAAGAGGCCGGCCCCCAAAACACTAAAAAAAACCACCCCAGGGGAGAACAAAATGCCGCCAACAGGTTTTCCCAGGCCGAAGGGGATGGGGCTGGTGCCGGGGGTTGTGGTAAAATGGAAAACGAATGGGATTGAGATAAAACCCCGGGCGGTAATCGCAAGATAGAAAGTTTAAAGCGGAGTTTTCGTTGCTTTAAAAACCCCCGCCTCTTCGGCCACTTTGGCCATTgcagca                                                                                                                       <  1:57672/239‑1 (MQ=255)
                                                                      tAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGGGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGAGGTAAGAGCACGATAGAAAGTCTTTAGACGATGGTTTGTTGCTTTATTGAACACCCCCTCTTCGGCCACCTTGGGCCGTGTAGCAGGa                                                                                                                    >  1:189378/1‑239 (MQ=255)
                                                                             aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTGATCGCCCATGCGAAAGGCTTTAAGGGGTGTTTGGTTGCTTTATAGAACCCCACCCCCTCGGGCACTTTGGGCCTTTCAGCAGGAGGGaaaa                                                                                                             >  1:118397/1‑239 (MQ=255)
                                                                             aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATTGAAAAGCTTAAGCGGATGTTTCGTTGCTTTATAGAAACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                             >  2:329878/1‑239 (MQ=255)
                                                                              tCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAAAACCGGCGGTAATCGCCCGAAAGAAAGACTTAAGCGGGAGGTTTCGTGCGTTTAAGACACCCGCCCCTTCGGCCCCTTTGGGCCTTGTACCAGGAAAGaaaaa                                                                                                            >  1:380114/1‑238 (MQ=255)
                                                                                            aacaCAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGa                                                                                                                                                                                             >  2:116476/1‑146 (MQ=255)
                                                                                            aacaCAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGa                                                                                                                                                                                             <  1:116476/146‑1 (MQ=255)
                                                                                                                   ggaaacaaactcgccccaCCAGTTTTCGCCAGGCAGGGGGGCTTTGTGTGGGCGCCAGGGGTTTGTGTAGTATAGAAGTGAAATGCGTGTGGAGTAAAACCACGGGCGGCAATCGCCAGAGAAAATTTCTAAAGGGTGTGTTTCGTGGTTTAAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctc                                                                          <  1:774844/220‑1 (MQ=255)
                                                                                                                                     ccggTTTGGCCAGGCCGAAGGGATTGTGGCCGGTCCCAGGGGTTGCTGTAATATAGAAGACAAAGGCGTTTGAGATGAAAACCCCGGGGGTAATCGCAAGAGAAAAAGTCTTAAGCGGATGTTTCTTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgacga                                                     <  2:853610/236‑1 (MQ=255)
                                                                                                                                                cATGCGGAGGGGCATTGTGCTTTGTCCGGTGTTTGGTATAATATAGAAGAAGAATGCGATTGAGATGAAAACCCCGGCGGTAATCGCCAGATAGAAAGTCTTAGGCGTATTTTTCGTTGCTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                           <  2:93523/239‑1 (MQ=255)
                                                                                                                                                    gcgatgggatttttgccctgggcgAGTGGTTGTGTGGAAAGAAAAGAGGAATGAGATGAGAATAAAACCCCGGGCGGTAATCGCCAGATAAAAATCCTTAAGCGGTTGTTTCGTGGCTTAAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                           <  1:106563/220‑1 (MQ=255)
                                                                                                                                                    gcGAGGGGCTTTGTGCCTGTGCCAGGGTTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCACAAGATGAAAATTCTTAAGCGATTTTTTCGTTGTTTTATAACCACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAg                                        <  2:118397/237‑1 (MQ=255)
                                                                                                                                                                     ctgtgccAGTGTTTGCTGGTAAATAGAAGACAAAGGCGTTTGAGATAAAACCCCCGGCGGTAACCGCAAGATAGAAAGTCTAAGGCGGATGTTTCGTTGCTTTAGAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAAt                       <  2:380114/239‑1 (MQ=255)
                                                                                                                                                                     ctgggcaAGGGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAACCCCCGGCGGTAACCGCAAGATAGAAAGTTTAAAGCGGATGTTTCGTTGCTTAAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAAt                       <  2:104355/239‑1 (MQ=255)
                                                                                                                                                                           cAGGGGTTGGTGTAATAAAAAAGGCGATTGCGATGGAAATGAAAACCCCGGCGGTAACCGCAAGATAGAAAGTCTAAAGGGGTTGTTGCGTTGCTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                  <  2:116694/238‑1 (MQ=255)
                                                                                                                                                                                   gCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTTGCCATTGCAGCAGGAAGTAAAAAATCAAAAGGGTTTGCAGCTTTTCCACTAACTCTCTCTTTTTTAAGTCGGGGACGGGATACTAACAAACCATTATAGCTgtcccttatacacacctccc          >  2:236766/1‑219 (MQ=255)
                                                                                                                                                                                    gtgtgATAAAGAAGAGAAAGGCGATTGAGATGAAAACCCCGGGGGTAATCGCCAGATAGAAAGTCTTAAGCGGTTTTTCCGTTGCTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCatat        <  2:741800/238‑1 (MQ=255)
                                                                                                                                                                                           taaagaagaCGAATGCGATTGAGATAAAACCCCGTGGGTTAACCGCAAGATAGAAAGTCTTAAGGGGTTGTTTCGTGTCTTTATAAACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc  <  2:405251/236‑1 (MQ=255)
                                                                                                                                                                                           taaaaaagaCAAAGGCGATTGATAGAAAAACCCCGGCGGTAACGGCAAGATAGAAATTCTTAAGCGGTTGTTTCTTTGCTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTACAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc  <  2:586636/234‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                  
CAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCC  >  NC_000913/954126‑954548

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