Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,549,091 C→T 7.7% P564L (CCG→CTG)  fdnG → formate dehydrogenase‑N, alpha subunit, nitrate‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,549,0910CT 7.7% 169.2 / 4.5 52P564L (CCG→CTG) fdnGformate dehydrogenase‑N, alpha subunit, nitrate‑inducible
Reads supporting (aligned to +/- strand):  ref base C (27/21);  new base T (2/2);  total (29/23)
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.50e-01

AACAACTGGGGCTATGACTGGCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTATTCAGACTG  >  NC_000913/1548961‑1549228
                                                                                                                                  |                                                                                                                                         
aacaacTGGGGCTATGACTGGCTGCCGAAGTGGGACCAGACCAACGACGTCATCAAGAATTTCAAAATGATGGATGAAGGAAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACaaaaa                                                                                                                                  >  2:1151342/1‑140 (MQ=255)
        gggCTATGACTGGCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAgtgg                                                                                                                          >  1:860607/1‑140 (MQ=255)
          gCTATGACTGGCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACACAAACAAAgtggtg                                                                                                                        >  2:641649/1‑140 (MQ=255)
          gCTATGACTGGCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAgtggtg                                                                                                                        >  2:641524/1‑140 (MQ=255)
            tATGACTGGCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAg                                                                                                                      <  2:500372/140‑1 (MQ=255)
             aTGACTGGCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGc                                                                                                                     >  1:295823/1‑140 (MQ=255)
                    gCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCTGTTGCGTCCTTCCTGGACAAAAACAAAGTGGTGAGCTGCCTGa                                                                                                              >  1:1153672/1‑140 (MQ=255)
                    gCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCTGTTGCGTCCTTCCTGGACAAAAACAAAGTGGTGAGCTGCCTGa                                                                                                              >  1:1153728/1‑140 (MQ=255)
                    gCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGa                                                                                                              >  1:333201/1‑140 (MQ=255)
                    gCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGa                                                                                                              >  2:485886/1‑140 (MQ=255)
                     cTGCCTAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACAAGATGGATGAAGGCAAAGTCACCGGTTATTTCTGACAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGat                                                                                                             <  2:1099044/140‑2 (MQ=255)
                     cTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAg                                                                                                             <  1:456924/140‑1 (MQ=255)
                             gTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGa                                                                                                     >  1:212261/1‑140 (MQ=255)
                               gggACCAGACCTACGACGTCATCAAGTATTTCATCATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAg                                                                                                   <  2:1134427/140‑1 (MQ=255)
                               gggACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAg                                                                                                   <  2:1050795/140‑1 (MQ=255)
                                     aGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAATGTGGTGAGCTGCCTGAGCAAGCTGTAGTACAtg                                                                                             <  1:282186/140‑1 (MQ=255)
                                     aGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACAtg                                                                                             <  1:282312/140‑1 (MQ=255)
                                      gACCTACTACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAATTCAACGGTTATTTCTGACAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAATTACAtgg                                                                                            >  1:395581/1‑140 (MQ=255)
                                                   atcaAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTg                                                                                                      >  2:375465/1‑117 (MQ=255)
                                                   atcaAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTg                                                                                                      >  2:375403/1‑117 (MQ=255)
                                                   atcaAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTg                                                                                                      <  1:375403/117‑1 (MQ=255)
                                                   atcaAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTg                                                                                                      <  1:375465/117‑1 (MQ=255)
                                                   atcaAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATcc                                                                               >  1:48797/1‑140 (MQ=255)
                                                        gTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTgg                                                                          <  2:324129/140‑1 (MQ=255)
                                                                                 aaaGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCa                                                 <  1:1125566/140‑1 (MQ=255)
                                                                                  aaGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACCTGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCGTCTGGCAg                                                >  2:482770/1‑140 (MQ=255)
                                                                                  aaGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAg                                                >  2:309390/1‑140 (MQ=255)
                                                                                  aaGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGACCTTCCCGGACCAAAACAAACTGCTGAGCTGCCTGAGCAAGCTGAAGTACATGGGGGATATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAg                                                >  2:483002/1‑140 (MQ=39)
                                                                                       aCCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCATAACCa                                           >  2:86122/1‑140 (MQ=255)
                                                                                        ccGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGAGAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCAc                                          >  1:366949/1‑140 (MQ=255)
                                                                                         cGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACg                                         <  2:137001/140‑1 (MQ=255)
                                                                                         cGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACg                                         >  2:720294/1‑140 (MQ=255)
                                                                                                        gggCTTTAACCCTGTTGCGTCCTTCCTGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATTGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATg                          <  2:1153672/140‑1 (MQ=255)
                                                                                                        gggCTTTAACCCTGTTGCGTCCTTCCTGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATTGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATg                          <  2:1153728/140‑1 (MQ=255)
                                                                                                               aaCCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATcc                   >  2:801717/1‑140 (MQ=255)
                                                                                                                   cGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGTTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCg               >  2:59259/1‑140 (MQ=255)
                                                                                                                    ggTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGt              >  1:493796/1‑140 (MQ=255)
                                                                                                                    ggTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGt              >  1:848510/1‑140 (MQ=255)
                                                                                                                     gTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGTGTCGAACGATGTCGATCCGGCggc             >  1:805881/1‑138 (MQ=255)
                                                                                                                     gTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGaa                               >  2:1013078/1‑122 (MQ=255)
                                                                                                                     gTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGaa                               <  1:1013078/122‑1 (MQ=255)
                                                                                                                     gTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTc             >  1:171912/1‑140 (MQ=255)
                                                                                                                     gTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTc             >  1:804113/1‑140 (MQ=255)
                                                                                                                         cGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTAtt         <  1:29973/140‑1 (MQ=255)
                                                                                                                             cgtCCCGGACAAACACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTATTCAGa     <  2:307466/138‑1 (MQ=255)
                                                                                                                             cTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTATTCAGa     <  2:1116182/140‑1 (MQ=255)
                                                                                                                             cTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTATTCAGa     <  1:269764/140‑1 (MQ=255)
                                                                                                                             cTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTATTCAGa     <  2:504497/140‑1 (MQ=255)
                                                                                                                             cTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTATTCAGa     <  2:104406/140‑1 (MQ=255)
                                                                                                                             cTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTATTCAGa     <  2:307758/140‑1 (MQ=255)
                                                                                                                             cTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACAATGTCGATCCGGCGTCTATTCAGa     <  2:403570/140‑1 (MQ=255)
                                                                                                                               tCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTATTCAGACt   >  2:1181656/1‑140 (MQ=255)
                                                                                                                                cccGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTATTCAGACTg  >  1:1167962/1‑140 (MQ=255)
                                                                                                                                  |                                                                                                                                         
AACAACTGGGGCTATGACTGGCTGCCGAAGTGGGACCAGACCTACGACGTCATCAAGTATTTCAACATGATGGATGAAGGCAAAGTCACCGGTTATTTCTGCCAGGGCTTTAACCCGGTTGCGTCCTTCCCGGACAAAAACAAAGTGGTGAGCTGCCTGAGCAAGCTGAAGTACATGGTGGTTATCGATCCGCTGGTGACTGAAACCTCTACCTTCTGGCAGAACCACGGTGAGTCGAACGATGTCGATCCGGCGTCTATTCAGACTG  >  NC_000913/1548961‑1549228

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: