Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,162,765 C→A 12.5% L247L (CTG→CTT)  pssA ← CDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,162,7650CA12.5% 89.5 / 3.3 32L247L (CTG→CTT) pssACDP‑diacylglycerol‑‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base C (13/15);  new base A (2/2);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GCGGTGCGGTCACCGACAATCAGCTCCACCCCAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATGCGCCGCAGCCGTGCGCGCAAACGCCGAATGTCACCGCGCAAACTGCGG  >  NC_002947/4162625‑4162909
                                                                                                                                             |                                                                                                                                                
gcggtgcggtCACCGACAATCAGCTCCACCCCAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCgg                                                                                                                                         >  2:378951/1‑150 (MQ=255)
    tgcggtCACCGACAATCAGCTCCACCCCAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCa                                                                                                                                     <  2:517607/150‑1 (MQ=255)
         tCACCGACAATCAGCTCCACCCCAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCa                                                                                                                                     >  1:162291/1‑145 (MQ=255)
         tCACCGACAATCAGCTCCACCCCAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCa                                                                                                                                     <  2:162291/145‑1 (MQ=255)
         tCACCGACAATCAGCTCCACCCCAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGa                                                                                                                                <  2:729584/150‑1 (MQ=255)
             cGACAATCAGCTCCACCCCACCCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGc                                                                                                                            <  1:729936/150‑1 (MQ=255)
                   tCAGCTCCACCCCAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTg                                                                                                                                  >  2:638792/1‑138 (MQ=255)
                   tCAGCTCCACCCCAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTg                                                                                                                                  <  1:638792/138‑1 (MQ=255)
                          cACCCCAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGAGCGCCACTGACAGCAAGGCGCACAGCCCCCGTTTCAGTGGATTGCCACGGCCCACAcc                                                                                                               >  2:235906/1‑150 (MQ=255)
                              ccAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTGCCAGCGGG‑TGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCcagc                                                                                                           <  1:235906/150‑1 (MQ=255)
                              ccAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCcagc                                                                                                           >  1:249022/1‑150 (MQ=255)
                                     cccGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCgg                                                                                                                                         <  2:469834/113‑1 (MQ=255)
                                     cccGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCgg                                                                                                                                         >  1:469834/1‑113 (MQ=255)
                                     cccGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGa                                                                                                    >  1:445627/1‑150 (MQ=255)
                                         cgcCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCAc                                                                                                <  2:445627/150‑1 (MQ=255)
                                                      tCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAGATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGc                                                                                   <  2:249022/150‑1 (MQ=255)
                                                        aaGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTg                                                                                                                                  >  2:720443/1‑101 (MQ=255)
                                                        aaGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTg                                                                                                                                  <  1:720443/101‑1 (MQ=255)
                                                           ttCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCcag                                                                                                            <  2:518229/120‑1 (MQ=255)
                                                           ttCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCcag                                                                                                            >  1:518229/1‑120 (MQ=255)
                                                                         ccTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTa                                                                <  1:155383/150‑1 (MQ=255)
                                                                           tG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTAcc                                                              >  2:396842/1‑150 (MQ=255)
                                                                                    gTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGc                                                      >  2:204367/1‑150 (MQ=255)
                                                                                                                   gtgcgCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATGCGCCGCAGCCGTGCGCGCAAACGCCGa                       >  1:297175/1‑150 (MQ=255)
                                                                                                                          ccACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATGCGCCGCAGCCGTGCGCGCAAACGCCGAATGTCAc                <  1:396842/150‑1 (MQ=255)
                                                                                                                              tGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATGCGCCGCAGCCGTgcgc                                  >  2:118802/1‑128 (MQ=255)
                                                                                                                              tGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATGCGCCGCAGCCGTgcgc                                  <  1:118802/128‑1 (MQ=255)
                                                                                                                                   gCAAGGCGCAAAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATgc                                                 <  2:458904/108‑1 (MQ=255)
                                                                                                                                   gCAAGGCGCAAAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATgc                                                 <  2:458674/108‑1 (MQ=255)
                                                                                                                                   gCAAGGCGCAAAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATgc                                                 >  1:458904/1‑108 (MQ=255)
                                                                                                                                   gCAAGGCGCAAAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATgc                                                 >  1:458674/1‑108 (MQ=255)
                                                                                                                                        gcgcACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATGCGCCGCAGCCGTGCGCGCAAACGCCGAATGTCACCGCGCAAACTGCgg  <  2:308564/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GCGGTGCGGTCACCGACAATCAGCTCCACCCCAACCCCCCGCGCCAGGGCATGGTCAAGTTCGCGCATCACGACCTG‑CGGCGGGTTGAAATACGGCGTACTGATGATGATCTGCGTGCGCGCCACTGCCAGCAAGGCGCACAGCGCCCGGTTCAGTGGATTGCCACGGCCCACACCCAGCAGCGGGATCACCTGCAGGCCCTGCCCCGTCACGCTGGCCGGTACCTCGTAGGCCATGCGCCGCAGCCGTGCGCGCAAACGCCGAATGTCACCGCGCAAACTGCGG  >  NC_002947/4162625‑4162909

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