Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,109,187 A→G intergenic (‑32/‑366) speC ← / → yqgA ornithine decarboxylase, biosynthetic/DUF554 domain‑containing protein YqgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,109,1870AG97.3% 125.6 / NA 37intergenic (‑32/‑366)speC/yqgAornithine decarboxylase, biosynthetic/DUF554 domain‑containing protein YqgA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (20/16);  minor base C (0/1);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCGTAAAATCAGTATCTCCCAACGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATACCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACTGTTGTCCTGCAGTGAGCGTGCGT  >  NC_000913/3109060‑3109311
                                                                                                                               |                                                                                                                            
gTCCGTAAAATCAGTATCTCCCAACGCCACCACGCGACGATGAGCAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGa                                                                                                                         <  2:470475/133‑1 (MQ=255)
  ccGTAAAATCAGTATCTCCCAACGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACa                                                                                                                       <  2:113171/133‑1 (MQ=255)
    gTAAAATCAGTATCTCCCAACGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAg                                                                                                                     <  1:47133/133‑1 (MQ=255)
     tAAAATCAGTATCTCCCAACGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAgg                                                                                                                    >  2:475336/1‑133 (MQ=255)
         aTCAGTACCTCCCAACGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGa                                                                                                                <  2:595132/133‑1 (MQ=255)
                  tcCCAACGCCCCCACCCGACGATGCGAATAAATTCGGGATATCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGc                                                                                                       <  2:182827/133‑1 (MQ=255)
                      aaCGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGGCAAGGGGGAAAACGCTGCCCCa                                                                                                   >  1:241529/1‑133 (MQ=255)
                      aaCGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCa                                                                                                   >  1:544762/1‑133 (MQ=255)
                        cGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGCATGCCGGCCAAGATGGAAAACCCTGCCCCATg                                                                                                 >  2:255986/1‑133 (MQ=255)
                        cGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATg                                                                                                 <  2:24004/133‑1 (MQ=255)
                        cGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATg                                                                                                 >  1:632921/1‑133 (MQ=255)
                            agcacGCGACTATGAGAAGAAAGGTGGGCTACTAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGat                                                                                             <  2:315380/131‑1 (MQ=255)
                             ccacGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGata                                                                                            >  2:543029/1‑133 (MQ=255)
                                cgcgACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACTCTGCCCCATCAGATATgg                                                                                         >  1:264530/1‑133 (MQ=255)
                                         gagaAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTg                                                                                <  1:571429/133‑1 (MQ=255)
                                                         tACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGAcc                                                                >  2:577586/1‑133 (MQ=255)
                                                            cAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGcc                                                             >  2:609179/1‑133 (MQ=255)
                                                             aGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCa                                                            <  2:544762/133‑1 (MQ=255)
                                                                    tactGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCa                                                     >  2:592918/1‑133 (MQ=255)
                                                                         gcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCGTGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTaaa                                                >  2:143013/1‑133 (MQ=255)
                                                                         gcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTaaa                                                >  1:90957/1‑133 (MQ=255)
                                                                         gcggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAACTCCGACCACCAGACATCAGTaaa                                                >  2:233086/1‑133 (MQ=255)
                                                                           ggcAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAg                                              >  1:16043/1‑133 (MQ=255)
                                                                              aaTATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCAACAGATACTGCCGCCAGTCATCCGTAAAAGCAg                                           >  2:305808/1‑133 (MQ=255)
                                                                                     aTTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCGGACGGCCGGACATCAGTAAAAGCAGAAACGct                                    >  2:194857/1‑133 (MQ=255)
                                                                                       tGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGctct                                  <  1:592918/133‑1 (MQ=255)
                                                                                               tAAGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTAc                          >  2:550914/1‑133 (MQ=255)
                                                                                                 aGCGCAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACtg                        <  2:16043/133‑1 (MQ=255)
                                                                                                   cgcAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGAATTTACtgtt                      >  1:76343/1‑133 (MQ=255)
                                                                                                     cAAACCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACtgttgt                    >  2:351300/1‑133 (MQ=255)
                                                                                                      aaaCCCGTTTCGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACTGTTGTc                   <  1:577586/133‑1 (MQ=255)
                                                                                                       aaCCCGTTTCGGGGACGTAAGAATCCCGGACACGGTGGAAAGCCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACTGTTGTcc                  <  2:195945/133‑1 (MQ=255)
                                                                                                                cGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACTGTTGTCCTGCAGTGAg         <  2:105207/133‑1 (MQ=255)
                                                                                                                cGGGGAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACTGTTGTCCTGCAGTGAg         <  1:475336/133‑1 (MQ=255)
                                                                                                                 ggggAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACTGTTGTCCTGCAGTGAgc        <  2:90957/133‑1 (MQ=255)
                                                                                                                 ggggAAGTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACTGTTGTCCTGCAGTGAgc        <  1:609179/133‑1 (MQ=255)
                                                                                                                       gTAAGAATGCCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACTTTTGTCCAGCAGTGAgcgtgcgt  >  1:231590/1‑133 (MQ=255)
                                                                                                                               |                                                                                                                            
GTCCGTAAAATCAGTATCTCCCAACGCCACCACGCGACGATGAGAAGAAAGTCGGGATACCAGTTCACTACTGGCGGCAATATTCATTGATTTCATAAGCGCAAACCCGTTTCGGGGAAGTAAGAATACCGGACAAGGTGGAAAACCCTGCCCCATGAGATATGGGTCAAACTGGTCACCAGCTCCGACCGCCAGACATCAGTAAAAGCAGAAACGCTCTGATTTTACTGTTGTCCTGCAGTGAGCGTGCGT  >  NC_000913/3109060‑3109311

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