Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,717,779 G→T 12.5% intergenic (+150/‑45) cti → / → PP_2377 esterified fatty acid cis/trans isomerase/acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,717,7790GT12.5% 88.1 / 2.3 32intergenic (+150/‑45)cti/PP_2377esterified fatty acid cis/trans isomerase/acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (15/13);  new base T (2/2);  total (17/15)
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.67e-01

GAGTACCTTGCGGTCGGAACGGGTGGGTGGTCTTGATCGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCTGCCACCATTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCAT  >  NC_002947/2717650‑2717924
                                                                                                                                  |                                                                                                                                                 
gAGTACCTTGCGGTCGGAACGGGTGGGTGGTCTTGATCGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACgg                                                                                                                                       <  1:141051/142‑1 (MQ=255)
gAGTACCTTGCGGTCGGAACGGGTGGGTGGTCTTGATCGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACgg                                                                                                                                       >  2:141051/1‑142 (MQ=255)
  gTACCTTGCGGTCGGAACGGGTGGGTGGTCTTGATCGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGc                                                                                                                             <  1:669552/150‑1 (MQ=255)
       ttGCGGTCGGAACGGGTGGGTGGTCTTGATCGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTTCACTCGATACGGGCCTCAATGCACCGc                                                                                                                        >  2:455801/1‑150 (MQ=255)
       ttGCGGTCGGAACGGGTGGGTGGTCTTGATCGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTTCACTCGATACGGGCCTCAATGCACCGc                                                                                                                        >  2:456030/1‑150 (MQ=255)
          cGGTCGGAACGGGTGGGTGGTCTTGATCGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCaa                                                                                                                     <  1:484866/150‑1 (MQ=255)
                      gtgggtggTCTTGATCGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACACCTTTCTGGAATAGGCTGTCGCTGTCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTAtt                                                                                                         >  1:693605/1‑150 (MQ=25)
                                  gATCGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGAATTCCGTCCTGTAtt                                                                                             >  2:232037/1‑150 (MQ=255)
                                     cGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTTCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGc                                                                                          <  1:455801/150‑1 (MQ=255)
                                     cGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTTCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGc                                                                                          <  1:456030/150‑1 (MQ=255)
                                            aaaGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTTTTTGCTTCAGGc                                                                                   >  1:409061/1‑150 (MQ=255)
                                                   gctgctTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAgg                                                                                                              >  1:337261/1‑116 (MQ=255)
                                                   gctgctTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAgg                                                                                                              <  2:337261/116‑1 (MQ=255)
                                                   gctgctTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCgg                                                                            <  1:68972/150‑1 (MQ=255)
                                                     tgctTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGc                                                                          >  2:398738/1‑150 (MQ=255)
                                                     tgctTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGc                                                                          >  1:672128/1‑150 (MQ=255)
                                                                ccTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCt                                                               >  2:58110/1‑150 (MQ=255)
                                                                  ttCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTgg                                                             >  2:448974/1‑150 (MQ=255)
                                                                            aaCCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCt                                                               >  1:314665/1‑138 (MQ=255)
                                                                            aaCCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCt                                                               <  2:314665/138‑1 (MQ=255)
                                                                            aaCCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCt                                                   <  1:148565/150‑1 (MQ=255)
                                                                            aaCCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCt                                                   <  2:672128/150‑1 (MQ=255)
                                                                                cGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCTGcca                                               <  1:58110/150‑1 (MQ=255)
                                                                                 gCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCTGccac                                              <  1:232037/150‑1 (MQ=255)
                                                                                      tACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGc                                                                                          >  1:627244/1‑101 (MQ=255)
                                                                                      tACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGc                                                                                          <  2:627244/101‑1 (MQ=255)
                                                                                         aacaGGTACGGCGCTGCTCTCGAGGCTGGCGCTGGCCTCTGCACTCGATACGAGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCTGCAACCATTTAAt                                      >  2:199896/1‑150 (MQ=255)
                                                                                              gTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCTGCCACCATTTc                                        >  1:744178/1‑143 (MQ=255)
                                                                                              gTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCTGCCACCATTTc                                        <  2:744178/143‑1 (MQ=255)
                                                                                                      ctgctCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCTGCCACCATTTCATGCAGATCttcttc                         <  1:448974/150‑1 (MQ=255)
                                                                                                              aaGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGAACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCTGCCACCATTTCATGCAGATCTTCTTCGACTTccc                 >  1:240117/1‑149 (MQ=255)
                                                                                                                              ctctGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCTGCCACCATTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCAt  >  1:329879/1‑150 (MQ=255)
                                                                                                                                  |                                                                                                                                                 
GAGTACCTTGCGGTCGGAACGGGTGGGTGGTCTTGATCGATTGCAAAGGGGGCTGCTTTGCAGCCCTTCGCGGGTGAACCCGCTCCTACAACAGGTACGGCGCTGCTCTCAAGGCTGGCGCTGGCCTCTGCACTCGATACGGGCCTCAATGCACCGCCAAGCGCAGGGTATTCCATGCTGTATTCGCTTCAGGCACTGCGGGCGTTCGCCGCCTGGGTGGTGGTCTGCCACCATTTCATGCAGATCTTCTTCGACTTCCACGCCACCGGCCCCAT  >  NC_002947/2717650‑2717924

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