Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 759,544 A→G 9.7% D280G (GAC→GGC)  sucA → 2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913759,5440AG 9.7% 112.3 / 3.4 47D280G (GAC→GGC) sucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base A (18/23);  new base G (2/3);  total (21/26)
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.59e-01

TTAATCCCGATGCTTAAAGAGATGATCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTA  >  NC_000913/759408‑759692
                                                                                                                                        |                                                                                                                                                    
ttAATCCCGATGCTTAAAGAGATGATCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCgg                                                                                                                                          <  1:256568/149‑1 (MQ=255)
  aaTCCCGATGCTCAATGAGATGATCCGCGACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACAGCAAGACTTGTTCGACGAGTTccccttta                                                                                                                                        <  2:785289/149‑5 (MQ=255)
   aTCCCGATGCTTAAAGAGATGATCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGGCGGTaa                                                                                                                                       >  2:731489/1‑149 (MQ=255)
   aTCCCGATGCTTAAAGAGATGATCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTaa                                                                                                                                       >  2:646999/1‑149 (MQ=255)
   aTCCCGATCCTGCAAGAGGCGCGCGGCCGCGCTGGCAACTGCCGCGCCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGGCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCggcca                                                                                                                                       <  2:653364/149‑4 (MQ=255)
    tCCCGATGCTTAAAGAGATGATCCGCCACGCTGGCAACAGCGGCACCCTCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTaaa                                                                                                                                      <  1:237419/149‑1 (MQ=255)
    tCCCGATGCTTAAAGAGATGATCCGCCACGCTGGCAACAACGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAACAAACCGCAAGACTTGTTCGCCGAGTTCGCCGGTaaa                                                                                                                                      >  1:359753/1‑149 (MQ=255)
         aTGCTCAAACACATGAGCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCCGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATaa                                                                                                                                 <  1:135110/149‑1 (MQ=255)
             ttAAAGAGATGATCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTGAAAGACCGCAAGGCTTGTTCGACGAGTTCGCCGGTAAACATAAAGaa                                                                                                                             >  2:641872/1‑149 (MQ=255)
               aaaGAGATGATCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAaca                                                                                                                           <  2:373644/149‑1 (MQ=255)
                gatggatgatCCGCCACGCTGGCAGAAGCGGCACCCGGGAAGTGGTTCTCGGGATGGCGCACCGTGGTAGTCTGATCGTGCTGGTGAACGTGCGGGGTAAAAAACCGGATGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAacac                                                                                                                          <  2:262006/145‑1 (MQ=255)
                    gatgatCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGACCGTGCTGGTTAAAAAACCGCATGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCg                                                                                                                      <  1:123332/149‑1 (MQ=255)
                      tgatCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGc                                                                                                                    >  1:711585/1‑149 (MQ=255)
                        atCCGCCACGCTGGCAACAGCGGCAACCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCAc                                                                                                                  >  1:443126/1‑149 (MQ=255)
                         tcccgcACGGTGGCGACGGCGGCACCCGCGAAGTGGTGTTCGGGATGGCGCACCGTGGTCGTCGGAACGTGCTGGTGAACGTGCTGGGTAAAAGACCGCAAGACTTGTTCGACGAGGTCGCCGGTAAACATAAAGAACACCTCGGCACg                                                                                                                 <  2:371581/144‑1 (MQ=255)
                         tCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCCCg                                                                                                                 >  2:303010/1‑149 (MQ=255)
                         tCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACg                                                                                                                 <  1:231947/149‑1 (MQ=255)
                         tCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACg                                                                                                                 <  1:79197/149‑1 (MQ=255)
                            gCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGGCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGGACACCTCGGCACGGGt                                                                                                              >  1:93772/1‑149 (MQ=255)
                             ccACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTg                                                                                                             >  1:852976/1‑149 (MQ=255)
                                            cGGCAGCCGCGAAGTGGTGCTCGGGATGGCGCAGCGTGGTCGCCTGAACGTGCTGGTGAACGCGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACcaca                                                                                              <  2:253988/149‑1 (MQ=255)
                                             ggCACCCGCGAAGTGGTTCCCGGGATGGCGCCCCGCGGTCGCCTGAACGCGCTGGTGCACGTGCTGGGTAAAAAACCGCAAGACTTGTGCGACGAGTTCGCCGGGATACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACAt                                                                                             <  2:807433/149‑1 (MQ=255)
                                              gCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATg                                                                                            >  2:300735/1‑149 (MQ=255)
                                               cACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATgg                                                                                           >  2:99239/1‑149 (MQ=255)
                                               cACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATgg                                                                                           <  1:184491/149‑1 (MQ=255)
                                                   cgcgAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCtt                                                                                       <  1:722464/149‑1 (MQ=255)
                                                            gTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGGCGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGa                                                                              >  1:181472/1‑149 (MQ=255)
                                                            gTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGGCGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGa                                                                              >  1:181940/1‑149 (MQ=255)
                                                             ttCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGAc                                                                             >  1:77493/1‑149 (MQ=255)
                                                                  gggATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCa                                                                        <  1:457182/149‑1 (MQ=255)
                                                                            cgcGTGGTCGTCTCGACGTGCTGGTGAACGGGCTGGGGAAAAAACCGCAAGACTTGTCCGACGAGTTTGCCGGTAAACATAAAGAACCCCGCGGCACGGGTGACGTGAAATACCACATGGGCTCCTCGTCTGACTTCCAGACCGATggc                                                              <  2:43753/147‑1 (MQ=255)
                                                                                          aacgtgctggtgaacgtgctggGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACct                                                >  1:742387/1‑149 (MQ=255)
                                                                                          aacgtgctggtgaacgtgctggGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCAAATGTGCTTCCCGTATGACTTCCAGGCCGATGGCGGCCTTGTGCAcgt                                                >  2:128093/1‑147 (MQ=255)
                                                                                             gtgctgggtaacgtgctggGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCCGACCGATGGCGGCCTGGTGCACctggc                                             <  2:403834/149‑1 (MQ=255)
                                                                                                         gtgctggGTAAAAAACCGGACGACTTGTTCGGCGAGTTCGCCGGTAGACGTAGAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTaa                                 <  2:181940/149‑1 (MQ=255)
                                                                                                         gtgctggGTAAAAAACCGCAAGACTTGTTCGGCGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTaa                                 <  2:181472/149‑1 (MQ=255)
                                                                                                          tgctggGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAAc                                >  1:77562/1‑149 (MQ=255)
                                                                                                              gggTAAAAAACCGCAAGACTGGTGCGAGGGGTTCGCCGGTAAACATAAGGGCCACCTAGGCACGGGTGACGGGGAATACCACAGGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGt                            <  2:93772/149‑1 (MQ=255)
                                                                                                               ggTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGtc                           <  2:706659/149‑1 (MQ=255)
                                                                                                                   aaaaaCCGCAAACCGTGGTCGACGGGGTCGCCGGTACACATAAAGAACCCCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCAc                       <  2:663860/149‑1 (MQ=255)
                                                                                                                    aaaaCCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACCCCTCGGCACGGGGGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCAcc                      >  2:142885/1‑149 (MQ=255)
                                                                                                                     aaaCCGCAAGACTTGGTCGACGAGTTCGCCGGGAAACATAAAGAACACCTCGGCACGGGTGCGGGGCAATACCACATGGGCAACTCGTCTGACTTCCAGACCGATGGCGGCGTGAGGCACCTGACGCTGGGGGTTAACGCGtctgccgt                     >  2:445004/1‑144 (MQ=255)
                                                                                                                          gCATGACTTGTTCGGCGAGTTCGCCGGTAAACATAAAGATCCCCTCGGCACGGGTGACGTGAAATACCACATGGGCTGCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTgaga                <  2:197643/149‑1 (MQ=255)
                                                                                                                          gCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTgaga                <  2:711585/149‑1 (MQ=255)
                                                                                                                              gACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGt            <  1:300735/149‑1 (MQ=255)
                                                                                                                               aCTTGTTCGACGAGTTCGCCGGTAAACATGAAGAACACCTCGGCACGGGTGACGTGAAATACCTCATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTa           >  1:557445/1‑149 (MQ=255)
                                                                                                                                 ttGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAg         >  1:346108/1‑149 (MQ=255)
                                                                                                                                        acgaGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGGCGTGAAATACCACATGAACCACAAGGTTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTa  >  1:377128/1‑149 (MQ=255)
                                                                                                                                        acgaGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTa  <  1:646999/149‑1 (MQ=255)
                                                                                                                                        acgaGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTa  <  1:731489/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
TTAATCCCGATGCTTAAAGAGATGATCCGCCACGCTGGCAACAGCGGCACCCGCGAAGTGGTTCTCGGGATGGCGCACCGTGGTCGTCTGAACGTGCTGGTGAACGTGCTGGGTAAAAAACCGCAAGACTTGTTCGACGAGTTCGCCGGTAAACATAAAGAACACCTCGGCACGGGTGACGTGAAATACCACATGGGCTTCTCGTCTGACTTCCAGACCGATGGCGGCCTGGTGCACCTGGCGCTGGCGTTTAACCCGTCTCACCTTGAGATTGTAAGCCCGGTA  >  NC_000913/759408‑759692

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