Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 327,990 T→C E248E (GAA→GAG)  betB ← betaine aldehyde dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913327,9900TC95.3% 141.2 / ‑1.3 43E248E (GAA→GAG) betBbetaine aldehyde dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (22/19);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATA  >  NC_000913/327848‑328132
                                                                                                                                              |                                                                                                                                              
ccGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAg                                                                                                                                        <  1:42423/151‑1 (MQ=255)
        aaGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCCAGATCGAGACCCGCATCATCGAAAAAGATCAGCGGTGATTTACCGCCCCGTTCCATGGTCACCTCTTTCAGGGAagag                                                                                                                                 >  1:9190/1‑150 (MQ=255)
                            acacacCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACt                                                                                                            <  2:317528/151‑1 (MQ=255)
                             cacacCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACtt                                                                                                           <  2:388808/151‑1 (MQ=255)
                                acCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACttttt                                                                                                        <  1:281578/151‑1 (MQ=255)
                                  cTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTg                                                                                                      >  1:82840/1‑151 (MQ=255)
                                  cTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGGGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACtttttt                                                                                                       >  1:397623/1‑150 (MQ=255)
                                   tGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGc                                                                                                     >  2:16578/1‑151 (MQ=255)
                                     aCCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGcc                                                                                                    >  1:592441/1‑150 (MQ=255)
                                         gAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTgg                                                                                               <  2:132064/151‑1 (MQ=255)
                                          aGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGc                                                                                              <  2:453674/151‑1 (MQ=255)
                                                agaagTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACAcc                                                                                        <  2:727966/151‑1 (MQ=255)
                                                  aagTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACAccgc                                                                                      >  2:377446/1‑151 (MQ=255)
                                                   agTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACAccgcc                                                                                     <  1:640763/151‑1 (MQ=255)
                                                       tGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGt                                                                                  >  1:96571/1‑150 (MQ=255)
                                                             catGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATg                                                                             <  1:698785/149‑1 (MQ=255)
                                                             catGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGac                                                                           >  2:298543/1‑151 (MQ=255)
                                                               tGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGacac                                                                         >  2:187107/1‑151 (MQ=255)
                                                                        cggcggCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAAt                                                                >  1:235618/1‑151 (MQ=255)
                                                                         ggcggcGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAAt                                                                >  1:747948/1‑150 (MQ=255)
                                                                            ggcgAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGccc                                                            <  1:463537/151‑1 (MQ=255)
                                                                             gcgAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCg                                                           <  2:84470/151‑1 (MQ=255)
                                                                             gcgAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCg                                                           >  1:542059/1‑151 (MQ=255)
                                                                              cgagatcgagatCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCgg                                                          >  1:13808/1‑151 (MQ=255)
                                                                               gagatcgagatcCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCgg                                                          >  1:520836/1‑150 (MQ=255)
                                                                                      agaccCGCATCATCGAAAACGATCAGCGGTGATTTCCCGCCCAGTTCCATGGTCCCCCCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCg                                                   <  1:547615/150‑1 (MQ=255)
                                                                                            gCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGat                                            <  2:560566/151‑1 (MQ=255)
                                                                                             catcatCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGata                                           <  1:282157/151‑1 (MQ=255)
                                                                                             catcatCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGat                                            >  2:335366/1‑150 (MQ=255)
                                                                                                 atcGAAAACGATCAGCGGTGATTTCCCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGc                                       <  1:716836/151‑1 (MQ=255)
                                                                                                   cGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGccc                                     >  1:540422/1‑151 (MQ=255)
                                                                                                         cGATCAGCGGTGATTTACCGCTCAGTTCCATGTTCACTTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCTCTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGtctc                               <  2:559018/151‑1 (MQ=255)
                                                                                                             cAGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCcgc                            <  2:521739/150‑1 (MQ=255)
                                                                                                              aGCGGTGATTTACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCgcgc                          >  1:705283/1‑151 (MQ=255)
                                                                                                                        tACCGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCgg                >  2:597808/1‑151 (MQ=255)
                                                                                                                          ccGCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGCCCGCGGATTTAGCCATCTCTTTTTGGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCTGTCTCCGCGCCCACGCCCGGc               <  1:383149/150‑1 (MQ=255)
                                                                                                                            gCCCAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAac            >  2:97720/1‑151 (MQ=255)
                                                                                                                             cccAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAac            >  2:491285/1‑150 (MQ=255)
                                                                                                                              ccAGTTCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAacac          <  2:565592/151‑1 (MQ=255)
                                                                                                                                   tCCATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTa      >  1:81680/1‑150 (MQ=255)
                                                                                                                                    ccATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTaa     >  2:372387/1‑150 (MQ=255)
                                                                                                                                    ccATGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTaa     <  2:748023/150‑1 (MQ=255)
                                                                                                                                      aTGGTCACCTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATa  <  2:58112/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CCGGAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATA  >  NC_000913/327848‑328132

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