Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 323,355 C→T P340S (CCA→TCA)  ykgF → ferridoxin‑like LutB family protein, putative electron transport chain YkgEFG component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913323,3550CT98.2% 179.6 / NA 57P340S (CCA→TCA) ykgFferridoxin‑like LutB family protein, putative electron transport chain YkgEFG component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/27);  minor base G (0/1);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGCGTC  >  NC_000913/323221‑323490
                                                                                                                                      |                                                                                                                                       
aTAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAAt                                                                                                                                     <  2:1528280/139‑1 (MQ=255)
aTAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTTCCAGGGTCAAtt                                                                                                                                    >  1:42564/1‑139 (MQ=255)
    cGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGt                                                                                                                                 <  1:1742224/139‑1 (MQ=255)
     gggCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTTCCAGGGTCAATTGGTGc                                                                                                                               <  2:42564/139‑1 (MQ=255)
       gCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGt                                                                                                                                 <  2:872711/136‑1 (MQ=255)
       gCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGt                                                                                                                                 >  1:872711/1‑136 (MQ=255)
                 gTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAttt                                                                                                                       <  1:808382/136‑1 (MQ=255)
                 gTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAttt                                                                                                                       >  2:808382/1‑136 (MQ=255)
                    cTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAGTTGGTGCGGTGATTTCTCCGc                                                                                                                 >  1:1019369/1‑139 (MQ=255)
                    cTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGc                                                                                                                 <  1:1747161/139‑1 (MQ=255)
                       gCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTAc                                                                                                              >  2:616886/1‑139 (MQ=255)
                           ctGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                                                          <  2:279191/139‑1 (MQ=255)
                            tGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTggc                                                                                                         >  1:275525/1‑139 (MQ=255)
                                   cGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCtata                                                                                                  >  1:564247/1‑139 (MQ=255)
                                        tgtgCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGACATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagat                                                                                             >  1:272526/1‑139 (MQ=255)
                                                                ggCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGc                                                                     <  2:422047/139‑1 (MQ=255)
                                                                ggCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGc                                                                     <  2:564247/139‑1 (MQ=255)
                                                                          aaTACTTGTCCGGCATATCGCCATATTGGCGGACATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGca                                                           <  2:272526/139‑1 (MQ=255)
                                                                            tACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAt                                                                                                                         >  2:694148/1‑75 (MQ=255)
                                                                            tACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAt                                                                                                                         <  1:694148/75‑1 (MQ=255)
                                                                                tGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCg                                                                                                                  <  2:1577746/78‑1 (MQ=255)
                                                                                tGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCg                                                                                                                  >  1:1577746/1‑78 (MQ=255)
                                                                                tGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCt                                                                       <  2:239524/121‑1 (MQ=255)
                                                                                tGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCt                                                                       >  1:239524/1‑121 (MQ=255)
                                                                                tGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCtt                                                     >  2:103941/1‑139 (MQ=255)
                                                                                     ggCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGAc                                                <  1:616886/139‑1 (MQ=255)
                                                                                      gCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAt                                                                                                                         >  1:499913/1‑65 (MQ=255)
                                                                                      gCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAt                                                                                                                         <  2:499913/65‑1 (MQ=255)
                                                                                         tatCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                                                   <  1:531881/100‑1 (MQ=255)
                                                                                         tatCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                                                   >  2:531881/1‑100 (MQ=255)
                                                                                                atatTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTcc                                     <  2:275525/139‑1 (MQ=255)
                                                                                                 tatTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTc                                      >  1:32738/1‑137 (MQ=255)
                                                                                                 tatTGGCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTc                                      <  2:32738/137‑1 (MQ=255)
                                                                                                      gCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAt                                                                                                                         >  1:485838/1‑49 (MQ=255)
                                                                                                      gCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGAt                                                                                                                         <  2:485838/49‑1 (MQ=255)
                                                                                                      gCGGTCATGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCg                               >  1:1638294/1‑139 (MQ=255)
                                                                                                            aTGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCtt                                                                 <  1:1182411/99‑1 (MQ=255)
                                                                                                            aTGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCtt                                                                 >  2:1182411/1‑99 (MQ=255)
                                                                                                            aTGGATATGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTcc                         >  1:1657012/1‑139 (MQ=255)
                                                                                                                atatGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                                                   >  1:19775/1‑77 (MQ=255)
                                                                                                                atatGGCTCTATTTATCCAGGGGCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                                                   <  2:19775/77‑1 (MQ=255)
                                                                                                                   tGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                                                          <  1:1192403/51‑1 (MQ=255)
                                                                                                                   tGGCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTgg                                                                                                          >  2:1192403/1‑51 (MQ=255)
                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagat                                                                                    <  1:1105366/72‑1 (MQ=255)
                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagat                                                                                    >  2:1105366/1‑72 (MQ=255)
                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCt                                                                       <  2:1787393/85‑1 (MQ=255)
                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCt                                                                       >  1:1787393/1‑85 (MQ=255)
                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGca                                                           >  1:182924/1‑97 (MQ=255)
                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGca                                                           <  2:182924/97‑1 (MQ=255)
                                                                                                                    ggCTCTATTTATCCAGGGTCAATTGGTACGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCaa                 >  1:1347637/1‑139 (MQ=255)
                                                                                                                      ctctATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCaaaa               >  1:564347/1‑139 (MQ=255)
                                                                                                                         tATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                         <  1:564291/62‑1 (MQ=255)
                                                                                                                         tATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                                                         >  2:564291/1‑62 (MQ=255)
                                                                                                                         tATTTATCCAGGGTCAATTGGTGCGGTGATTTCTCCGCTACTTGACGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTg            >  1:1334034/1‑139 (MQ=255)
                                                                                                                                   gggTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGCg    <  2:469467/137‑1 (MQ=255)
                                                                                                                                   gggTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGCGTc  <  2:1657012/139‑1 (MQ=255)
                                                                                                                                   gggTCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGCGTc  >  1:790170/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                       
ATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGGTGCGTATTCCGCTGTCAAAACTGATTTTGCGTC  >  NC_000913/323221‑323490

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: