Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,062 T→A 11.4% intergenic (+29/+199) glxK → / ← allE glycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,0620TA11.4% 136.6 / 6.3 52intergenic (+29/+199)glxK/allEglycerate 2‑kinase 2/(S)‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base A (3/3);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTA  >  NC_000913/542918‑543207
                                                                                                                                                |                                                                                                                                                 
ttttGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGc                                                                                                                                              <  1:989934/150‑1 (MQ=255)
   tGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACAACCCTTACCTTTAAAGCGGATGCGATATAAATCGCAtcc                                                                                                                                          <  2:1427059/150‑3 (MQ=255)
     ccTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAGATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATATAAATCGCAtccg                                                                                                                                         <  2:1337008/150‑4 (MQ=255)
        cGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGt                                                                                                                                      <  1:190854/150‑1 (MQ=255)
                aCCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACt                                                                                                                              <  2:19974/150‑1 (MQ=255)
                 ccTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTc                                                                                                                             >  1:134525/1‑150 (MQ=255)
                       gCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCg                                                                                                                       >  1:1041226/1‑150 (MQ=255)
                        ccGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCgg                                                                                                                      >  1:61886/1‑150 (MQ=255)
                         cGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGa                                                                                                                     <  2:578876/150‑1 (MQ=255)
                            tgTGCTCGCCTGCGGTGATACCAATCTCTTCAACAGCGCGCGAAATATTGGCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGc                                                                                                                  <  1:807137/149‑1 (MQ=255)
                            aGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGc                                                                                                                  <  1:1122219/150‑1 (MQ=255)
                               gCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGc                                                                                                               <  1:1032070/150‑1 (MQ=255)
                                cTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCa                                                                                                              <  1:1197430/150‑1 (MQ=255)
                                  cGCTAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCGATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTTGCTTAAAAGCGGTTGCGATTAAAATCGCATAAGAGTGCAGTTCTGATGCGGGATGCGGCATg                                                                                                            <  2:1394174/150‑1 (MQ=255)
                                    ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAg                                                                                                          >  1:406395/1‑150 (MQ=255)
                                    ccAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAg                                                                                                          >  1:769796/1‑150 (MQ=255)
                                      aGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTa                                                                                                        >  2:174944/1‑150 (MQ=255)
                                      aGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGCATCCGGCATGAGTa                                                                                                        <  1:148123/150‑1 (MQ=255)
                                       gCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAc                                                                                                       <  1:153862/150‑1 (MQ=255)
                                        cGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAcc                                                                                                      <  2:1368790/150‑1 (MQ=255)
                                           tGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCata                                                                                                   >  2:450182/1‑150 (MQ=255)
                                             aaaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtc                                                                                                 <  1:575706/150‑1 (MQ=255)
                                              aaCCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcc                                                                                                >  1:1217745/1‑150 (MQ=255)
                                                ccAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATAtcctt                                                                                              >  1:261647/1‑150 (MQ=255)
                                                    tctcTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTg                                                                                          <  2:916692/150‑1 (MQ=255)
                                                       ctTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGaaa                                                                                       >  2:1352955/1‑150 (MQ=255)
                                                             catcgcgcgAAATATTGACTGCGCCATTAAAATAGTTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATcgcg                                                                                 <  2:331843/147‑1 (MQ=255)
                                                                        atatTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtata                                                                      <  1:1352955/150‑1 (MQ=255)
                                                                         aatTGCCTGCGCCATTAAAATAGGTCATGGAATTAAAAACTAACCCTTACCTTTAAATCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCtatat                                                                     <  2:772502/149‑1 (MQ=255)
                                                                            tGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATAtt                                                                  >  1:1217939/1‑150 (MQ=255)
                                                                                  cgcCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCagag                                                            <  1:234391/150‑1 (MQ=255)
                                                                                     cATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                         >  2:1175834/1‑150 (MQ=255)
                                                                                     cATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATc                                                         >  2:295492/1‑150 (MQ=255)
                                                                                      aTTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATTCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTAGCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCa                                                        <  2:1041226/150‑1 (MQ=255)
                                                                                         aaaaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGt                                                     <  2:1217745/150‑1 (MQ=255)
                                                                                         aaaaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGt                                                     <  2:134525/150‑1 (MQ=255)
                                                                                          aaaTAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTa                                                    >  2:1230984/1‑150 (MQ=255)
                                                                                               ggTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCt                                               >  2:295889/1‑150 (MQ=255)
                                                                                                    gggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaa                                          >  2:186516/1‑150 (MQ=255)
                                                                                                     ggAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaag                                         <  1:1002712/150‑1 (MQ=255)
                                                                                                       aaTTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATaagcg                                       >  2:254698/1‑150 (MQ=255)
                                                                                                          tAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTATGCCTGATaagcgaag                                    >  2:284955/1‑150 (MQ=255)
                                                                                                             aaaCTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAgcgc                                 >  2:1297772/1‑150 (MQ=255)
                                                                                                                cTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATc                              >  1:24372/1‑150 (MQ=255)
                                                                                                                    cccTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGc                          <  2:406395/150‑1 (MQ=255)
                                                                                                                           cTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTa                   >  1:1243424/1‑150 (MQ=255)
                                                                                                                              tAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaa                >  2:1292299/1‑150 (MQ=255)
                                                                                                                                 ctcttaTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaa             >  1:621289/6‑150 (MQ=255)
                                                                                                                                   cttaTGCGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACaaaaaaa           >  1:51003/4‑150 (MQ=255)
                                                                                                                                         cGATATAAATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACg     >  1:491838/1‑150 (MQ=255)
                                                                                                                                            tttATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCTCATCAGGCAATGTTACAAAAAAAGCCACGGta  <  2:327919/150‑1 (MQ=255)
                                                                                                                                            tttATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGta  <  1:294609/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTTTGCCTCGTCTGGCACCTTTAGCCGAAGTGCTCGCCAGCGGTGAAACCAATCTCTTCAACAGCGCGCGAAATATTGCCTGCGCCATTAAAATAGGTCAGGGAATTAAAAACTAACCCTTACCTTTAAAGCGGATGCGATTTATATCGCATAAGAGTGCAGTACTCATGCCGGATGCGGCATGAGTACCATATCCTTCCTGAAAATCGCGCAAATTCTATATATTGCAGAGATCATGTAGGCCTGATAAGCGAAGCGCATCAGGCAATGTTACAAAAAAAGCCACGGTA  >  NC_000913/542918‑543207

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