Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 656,095 A→C intergenic (‑127/‑462) dcuC ← / → pagP anaerobic C4‑dicarboxylate transporter DcuC/Lipid IVA palmitoyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913656,0950AC97.1% 121.7 / NA 35intergenic (‑127/‑462)dcuC/pagPanaerobic C4‑dicarboxylate transporter DcuC/Lipid IVA palmitoyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (13/21);  minor base G (1/0);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAATTTTTCCTGTCTCCAGGCCCCAAAGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGAAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATCACATAATAATGTGTGCAAATTGATGGTGGATTGATGCAAATTTGTTAATTGCAATGGTGTTAGC  >  NC_000913/655966‑656215
                                                                                                                                 |                                                                                                                        
cATAATTTTTCCTGTCTCCAGGCCCCAAAGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGt                                                                                                                       >  1:396946/1‑133 (MQ=255)
   aaTTTTTCCTGTCTCCAGGCCCCAAAGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGc                                                                                                                    <  2:107564/133‑1 (MQ=255)
     tttttCCTGTCTCCAGGCCCCAAAGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTg                                                                                                                  >  1:108979/1‑133 (MQ=255)
     tttttCCTGTCTCCAGGCCCCAAAGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTg                                                                                                                  <  1:327147/133‑1 (MQ=255)
     tttttCCTGTCTCCAGGCCCCAAAGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTg                                                                                                                  <  1:356470/133‑1 (MQ=255)
     tttttCCTGTATCCAAGCCCCAAAATAAGTAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTg                                                                                                                  <  2:257456/133‑1 (MQ=255)
         tCCTGTCTCCAGGCCCCAAAGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGtttt                                                                                                              <  2:82814/133‑1 (MQ=255)
              tctcCAGGCCCCAAAGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATAc                                                                                                         <  1:609814/133‑1 (MQ=255)
               ctcCAGGCCCCAAAGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACa                                                                                                        >  1:137346/1‑133 (MQ=255)
                          aaaGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACaaa                                                                                             >  1:347168/1‑133 (MQ=255)
                          aaaGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACaaa                                                                                             >  1:433218/1‑133 (MQ=255)
                          aaaGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACaaa                                                                                             >  2:84344/1‑133 (MQ=255)
                          aaaGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACaaa                                                                                             >  2:577929/1‑133 (MQ=255)
                          aaaGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACaaa                                                                                             >  2:565261/1‑133 (MQ=255)
                          aaaGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACaaa                                                                                             >  2:499105/1‑133 (MQ=255)
                                      aaaaaaTTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGGAGTTGCTGGTTTAATACATAGATAACAAAACGCATATTTAg                                                                                 >  2:195590/1‑133 (MQ=255)
                                          aaTTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCgga                                                                             <  1:499105/133‑2 (MQ=255)
                                          aaTTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGt                                                                             <  2:137346/133‑1 (MQ=255)
                                                 aaaGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAAt                                                                      <  2:146493/133‑1 (MQ=255)
                                                   aGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTa                                                                    <  1:409464/133‑1 (MQ=255)
                                                   aGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTa                                                                    <  2:73033/133‑1 (MQ=255)
                                                   aGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTa                                                                    <  2:529548/133‑1 (MQ=255)
                                                   aGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTa                                                                    <  1:33248/133‑1 (MQ=255)
                                                   aGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTa                                                                    <  1:577929/133‑1 (MQ=255)
                                                     cTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATc                                                                  <  2:396946/133‑1 (MQ=255)
                                                     cTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATc                                                                  <  2:433218/133‑1 (MQ=255)
                                                     cTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATc                                                                  <  1:565261/133‑1 (MQ=255)
                                                       tAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATcac                                                                <  2:347168/133‑1 (MQ=255)
                                                          ggAAAAAATATGCCCTATAAATTGTCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATCACata                                                             <  2:264443/133‑1 (MQ=255)
                                                          ggAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATCACata                                                             <  2:347115/133‑1 (MQ=255)
                                                               aaaTATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATCACATaataat                                                        >  1:311617/1‑133 (MQ=255)
                                                                     gCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCTGTTGGCAATTATCACATAATAATGTGTGc                                                  >  1:603039/1‑133 (MQ=255)
                                                                                            cTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATCACATAATAATGTGTGCAAATTGATGGTGGATTGATGCa                            >  2:29812/1‑132 (MQ=255)
                                                                                            cTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATCACATAATAATGTGTGCAAATTGATGGTGGATTGATGCa                            <  1:29812/132‑1 (MQ=255)
                                                                                                                     aTTAAGACATTGCAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATCACATAATAATGTGTGCAAATTGATGGTGGATTGATGCAAATTTGTTAATTGCAATGGTGTTAGc  >  1:150237/1‑133 (MQ=255)
                                                                                                                                 |                                                                                                                        
CATAATTTTTCCTGTCTCCAGGCCCCAAAGTAAATAATAAAAAATTCCTAAAGCTTAAGGAAAAAATATGCCCAATAAATTGGCGATGAATGCTGATTAAAATCAAGAAAAACTGCCATTAAGACATTGAAGTTGCTGTTTTTATACATAGATAACAAAACGCATATTTAGCGGTTGGCAATTATCACATAATAATGTGTGCAAATTGATGGTGGATTGATGCAAATTTGTTAATTGCAATGGTGTTAGC  >  NC_000913/655966‑656215

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: