Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 911,132 A→G intergenic (‑83/+50) poxB ← / ← hcr pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding/HCP oxidoreductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913911,1320AG98.3% 210.3 / ‑6.8 60intergenic (‑83/+50)poxB/hcrpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding/HCP oxidoreductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (30/29);  minor base C (0/1);  total (30/30)
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

CACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCA  >  NC_000913/910996‑911264
                                                                                                                                          |                                                                                                                                    
cACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                                    <  1:1062813/139‑1 (MQ=255)
                       tGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                                    <  1:42261/116‑1 (MQ=255)
                       tGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                                    >  2:42261/1‑116 (MQ=255)
                            atatAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAat                                                                                                        <  1:1681106/139‑1 (MQ=255)
                                     gCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAGATAGCCCGATGAAAGGAATATCATCggg                                                                                               <  2:1164189/139‑1 (MQ=255)
                                            tttgtttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                                    >  2:1345180/1‑95 (MQ=255)
                                                tttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGAAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                                    >  1:247734/1‑91 (MQ=255)
                                                          ctcCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGAtaaggggaaaa                                                                          >  2:1737823/1‑129 (MQ=255)
                                                                  ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                                             <  1:756401/112‑1 (MQ=255)
                                                                  ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                                             >  2:756401/1‑112 (MQ=255)
                                                                   ctaatgttaaTAACGTAAAAAATTTAATTCCTTCGACCGAGGGGGAAAGGAATGGAAGAAAAAGAGGCCCCAACCGTTAAATTACCCCATTAAAGGAATTTCCACCGGCATAAAGCCATTATGCGAGAACCAAATcccc                                                                <  1:1737823/130‑1 (MQ=255)
                                                                    tGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATcccc                                                                >  1:677317/1‑139 (MQ=255)
                                                                      aaTGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                                    <  1:1739497/70‑1 (MQ=255)
                                                                       aTGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTg                                                             >  1:71172/1‑139 (MQ=255)
                                                                          tgATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCgg                                                          <  1:409069/139‑1 (MQ=255)
                                                                           gATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCt                                                                                                                                    <  2:1040696/64‑1 (MQ=255)
                                                                              aaCGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACg                                                                                                                                >  2:1440174/1‑65 (MQ=255)
                                                                              aaCGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACg                                                                                                                                <  1:1440174/65‑1 (MQ=255)
                                                                                  gTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAgg                                                  >  1:1438346/1‑139 (MQ=255)
                                                                                    aaCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                                                >  2:590327/1‑139 (MQ=255)
                                                                                    aaCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAg                                                <  2:71172/139‑1 (MQ=255)
                                                                                     aCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTaaa                                                                                                                          <  1:293677/64‑1 (MQ=255)
                                                                                     aCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTaaa                                                                                                                          >  2:293677/1‑64 (MQ=255)
                                                                                        aGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAgg                                            >  1:790849/1‑139 (MQ=255)
                                                                                           ttAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                                             <  1:642981/87‑1 (MQ=255)
                                                                                           ttAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCa                                                                                             >  2:642981/1‑87 (MQ=255)
                                                                                               ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                   >  1:1841635/1‑77 (MQ=255)
                                                                                               ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                   <  2:1841635/77‑1 (MQ=255)
                                                                                                 cATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                   <  1:22396/75‑1 (MQ=255)
                                                                                                 cATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                   >  2:22396/1‑75 (MQ=255)
                                                                                                   tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGc                                                                                                                      <  2:372260/54‑1 (MQ=255)
                                                                                                   tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGc                                                                                                                      >  1:372260/1‑54 (MQ=255)
                                                                                                    cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCg                                                                                                                   >  1:597511/1‑56 (MQ=255)
                                                                                                    cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCg                                                                                                                   <  2:597511/56‑1 (MQ=255)
                                                                                                    cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                              >  1:31279/1‑61 (MQ=255)
                                                                                                    cTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                              <  2:31279/61‑1 (MQ=255)
                                                                                                        cGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTc                            >  2:436156/1‑139 (MQ=255)
                                                                                                         ggAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAAt                                                                    <  1:547121/98‑1 (MQ=255)
                                                                                                         ggAGGGGGAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAAt                                                                    >  2:547121/1‑98 (MQ=255)
                                                                                                             gggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                   <  1:585820/64‑1 (MQ=255)
                                                                                                             gggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATatcat                                                                                                   >  2:585820/1‑64 (MQ=255)
                                                                                                              ggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgag                                                                            >  1:235149/1‑86 (MQ=255)
                                                                                                              ggggAAGGGATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgag                                                                            <  2:235149/86‑1 (MQ=255)
                                                                                                                      ggATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGc                >  2:462013/1‑139 (MQ=255)
                                                                                                                       gATGGGAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCg               >  2:1302533/1‑139 (MQ=255)
                                                                                                                          gggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgaga                                                                           <  2:1412210/76‑1 (MQ=255)
                                                                                                                          gggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCgaga                                                                           >  1:1412210/1‑76 (MQ=255)
                                                                                                                           ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                              <  2:260324/40‑1 (MQ=255)
                                                                                                                           ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGaa                                                                                                              >  1:260324/1‑40 (MQ=255)
                                                                                                                           ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAata                                                                                                       <  2:546545/47‑1 (MQ=255)
                                                                                                                           ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAata                                                                                                       >  1:546545/1‑47 (MQ=255)
                                                                                                                           ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGc                                                     >  2:345099/1‑97 (MQ=255)
                                                                                                                           ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGc                                                     <  1:345099/97‑1 (MQ=255)
                                                                                                                           ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCgg           >  1:213212/1‑139 (MQ=255)
                                                                                                                           ggAGAGAAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGt          <  1:927807/139‑1 (MQ=255)
                                                                                                                            gagagaAAGGAGGCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGt          <  2:790849/139‑1 (MQ=255)
                                                                                                                                   aggaggCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCa                                         <  2:990664/101‑1 (MQ=255)
                                                                                                                                   aggaggCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCa                                         >  1:990664/1‑101 (MQ=255)
                                                                                                                                    ggaggCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCa  <  2:1438346/139‑1 (MQ=255)
                                                                                                                                    ggaggCGCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCa  <  1:462013/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                    
CACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCA  >  NC_000913/910996‑911264

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