Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,725,845 A→C 8.7% L773L (CTT→CTG)  plsB ← glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,725,8450AC 8.7% 138.1 / 5.5 46L773L (CTT→CTG) plsBglycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (18/24);  new base C (2/2);  total (20/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 = 8.16e-01

CCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGA  >  NC_002947/1725698‑1725990
                                                                                                                                                   |                                                                                                                                                 
cccGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGaa                                                                                                                                                  <  5:159548/149‑1 (MQ=255)
cccGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGaa                                                                                                                                                  <  5:255661/149‑1 (MQ=255)
      gtcACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAg                                                                                                                                             <  2:121301/148‑1 (MQ=255)
         aCGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTTCAGCAAggcgt                                                                                                                                         <  4:78135/149‑4 (MQ=255)
         aCGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCC‑‑CGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtg                                                                                                                                        >  3:36472/1‑148 (MQ=255)
         aCGGTGCAGGGCCACCGGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCt                                                                                                                                         <  8:168031/149‑1 (MQ=255)
              gCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaatg                                                                                                                                    <  3:207612/149‑1 (MQ=255)
                aGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGACGACAGTACTCGCTTGGCCACGCCCTAGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCtgaatg                                                                                                                                    >  2:158528/1‑147 (MQ=255)
                      aCCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCg                                                                                                                            <  7:64647/149‑1 (MQ=255)
                       ccTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCgg                                                                                                                           <  1:97295/149‑1 (MQ=255)
                             ggATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACa                                                                                                                       >  8:172570/1‑147 (MQ=255)
                                tCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTc                                                                                                                  <  3:241852/149‑1 (MQ=255)
                                   aCAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTg                                                                                                               <  3:172378/149‑1 (MQ=255)
                                       gcgcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCg                                                                                                            >  4:195167/1‑148 (MQ=255)
                                         gcAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCgaag                                                                                                         <  7:172570/149‑1 (MQ=255)
                                             cTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACt                                                                                                     >  3:220521/1‑149 (MQ=255)
                                             cTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACt                                                                                                     >  3:91067/1‑149 (MQ=255)
                                                  ccGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                <  5:8019/149‑1 (MQ=255)
                                                       aGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                <  3:189306/144‑1 (MQ=255)
                                                       aGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGt                                                                                                >  4:189306/1‑144 (MQ=255)
                                                        gTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGtt                                                                                           >  2:134959/1‑148 (MQ=255)
                                                          aCTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGa                                                                                                                          <  6:230603/115‑1 (MQ=255)
                                                          aCTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGa                                                                                                                          >  5:230603/1‑115 (MQ=255)
                                                            tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGc                                                                                      >  2:181756/1‑149 (MQ=255)
                                                            tCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGc                                                                                      >  1:255009/1‑149 (MQ=255)
                                                              gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCa                                                                                         <  1:203273/144‑1 (MQ=255)
                                                              gCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCa                                                                                         >  2:203273/1‑144 (MQ=255)
                                                                   gCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAg                                                                               <  4:220521/149‑1 (MQ=255)
                                                                        gCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATc                                                                           >  7:47237/1‑148 (MQ=255)
                                                                          ccTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGAc                                                                        <  7:146666/149‑1 (MQ=255)
                                                                              ggCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAg                                                                      <  1:181756/147‑1 (MQ=255)
                                                                                   gCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTgg                                                               <  2:255009/149‑1 (MQ=255)
                                                                                   gCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTgg                                                               <  3:195167/149‑1 (MQ=255)
                                                                                        ccGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCcat                                                          <  1:158528/149‑1 (MQ=255)
                                                                                               tGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACa                                                                       >  8:42154/1‑129 (MQ=255)
                                                                                               tGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACa                                                                       <  7:42154/129‑1 (MQ=255)
                                                                                                   gTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGca                                                 >  7:7491/1‑147 (MQ=255)
                                                                                                      aTAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAgg                                             <  8:47237/148‑1 (MQ=255)
                                                                                                             aGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAGCTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCa                                      <  2:252746/148‑1 (MQ=255)
                                                                                                                    cttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCg                              <  3:129985/149‑1 (MQ=255)
                                                                                                                     ttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAgc                                    >  3:59116/1‑142 (MQ=255)
                                                                                                                     ttgcGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAgc                                    <  4:59116/142‑1 (MQ=255)
                                                                                                                              ggCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGc                    <  4:91067/149‑1 (MQ=255)
                                                                                                                                      tcagaccTTGCTTCAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGATTCGGCGCTCAGGCAGTTCTGc            >  5:184095/6‑149 (MQ=255)
                                                                                                                                         gaccTTGCTTCAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTg             <  3:257762/143‑1 (MQ=255)
                                                                                                                                         gaccTTGCTTCAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTg             >  4:257762/3‑145 (MQ=255)
                                                                                                                                             ttGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGt     >  4:92781/1‑149 (MQ=255)
                                                                                                                                                 tGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGa  >  1:229306/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCCGTCGTCACGGTGCAGGGCCACCTGGCGGATCGACAGGCGCAGCTCGGCCGACAGTACTCGCTTGGCCACGCCCTCGGCCAGCTCGCCGAGCTTGTCGTGATAACTCAGCTTGCCTTGCGCGTCGGCGTGCAGCACGCCTTGCTGAAGCAAGGTCTGAATGAAGTGGCGGAACAGCGTCTTGTCGAAGAACTCCGGTGCGTTCAGGCCATGCAGGATCGACAGGCGCTGGGCCATCATCACGCACAGGTCTTCCAGCGCTTCGGCGCTCAGGCTGTTCTGCCCGCTGTTGA  >  NC_002947/1725698‑1725990

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