Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,001,304 (G)7→8 coding (960/2259 nt) xdhA → putative xanthine dehydrogenase molybdenum‑binding subunit XdhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,001,2971.G96.8% 101.2 / NA 32A318G (GCG→GGG) xdhAputative xanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (18/12);  minor base A (1/0);  total (19/13)
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

TCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTT  >  NC_000913/3001159‑3001427
                                                                                                                                           |                                                                                                                                   
gcgtgctccGTGCCCCCCGCGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCGGTCTAACACCGGCGGTTATGCATCTCGCGGGCGCTCCATCGCTTCTGCg                                                                                                                                    <  2:104622/133‑1 (MQ=255)
   cAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggg                                                                                                                                 >  2:157845/1‑140 (MQ=255)
     aCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTAATGCATCCCACGGGCACTCCATCGCTTCTGCgggggg                                                                                                                              >  1:304343/1‑140 (MQ=255)
                ccACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTGATGCATCTTACGAGCAGTCCATCGCTTCTGCAGGGGGGGGATAAAGTc                                                                                                                    >  2:321742/1‑140 (MQ=255)
                   cGCTTTTACCATTGACCGGCAACTGGGCGTGTACCGCGAGGGAACATTCAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                 <  2:316187/140‑1 (MQ=255)
                     cTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTa                                                                                                               <  2:373181/140‑1 (MQ=255)
                        ttACCATTGACGGGCTAATGGGAGTGAACAGCGACGGTACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCt                                                                                                            <  1:313917/140‑1 (MQ=255)
                                           gggCGTGAACCGCGACCAAACATTGAAAGGTCATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTAc                                                                                         <  2:304343/140‑1 (MQ=255)
                                                 gAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTAc                                                                                   <  2:17974/140‑1 (MQ=255)
                                                             aaCATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGAcc                                                                       <  1:636183/140‑1 (MQ=255)
                                                                          tataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAcc                                                          >  2:359514/1‑140 (MQ=255)
                                                                          tataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAcc                                                          >  2:146290/1‑140 (MQ=255)
                                                                           ataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCt                                                         >  2:514121/1‑140 (MQ=255)
                                                                            taGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTc                                                        >  2:139832/1‑140 (MQ=255)
                                                                            taGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGACTACGCTAACAATTAAAAGACTTGCTATACCAACCTc                                                        >  1:238695/1‑140 (MQ=255)
                                                                              gTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTccc                                                      >  2:17620/1‑140 (MQ=255)
                                                                                  ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCg                                                  >  1:21739/1‑140 (MQ=255)
                                                                                  ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTCCCTATGCTCACATTTCAAAGATCTGCGATACCAACCTCCCCTCg                                                  >  2:280330/1‑140 (MQ=255)
                                                                                   gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTGGTGCCTACACTTGCAGGTCAAAGACCTGCTATACCAACCTCCCCTCgg                                                 >  2:206358/1‑140 (MQ=255)
                                                                                     tGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCt                                               <  1:37635/140‑1 (MQ=255)
                                                                                         cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGATCTGGTAGACCAACGTCCCCTCGGCTGGTg                                           >  1:182871/1‑140 (MQ=255)
                                                                                         cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTg                                           <  1:359514/140‑1 (MQ=255)
                                                                                           gTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCg                                         <  2:88552/140‑1 (MQ=255)
                                                                                           gTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCg                                         <  1:466092/140‑1 (MQ=255)
                                                                                              tAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                      >  2:118289/1‑140 (MQ=255)
                                                                                                      gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTa                              >  1:482179/1‑140 (MQ=255)
                                                                                                      gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTAAACCAACCTCCCCGCGGCTAGTGCGATGCGTGGTTa                              >  2:309145/1‑140 (MQ=255)
                                                                                                      gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCAATACCAAACTCACCCCGGCTGGTGCGATGCGTGGTTa                              >  1:227306/1‑140 (MQ=255)
                                                                                                          ttATGCATCTCACGGTCACTTCTTCTCTTGTGCGGGGCGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGgc                          <  1:205055/140‑1 (MQ=255)
                                                                                                                aTCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCca                    <  2:58031/140‑1 (MQ=255)
                                                                                                                      cGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGATTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTc              >  1:475775/1‑140 (MQ=255)
                                                                                                                       gggCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCg             >  2:602768/1‑140 (MQ=255)
                                                                                                                            cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTAttt        >  2:397082/1‑140 (MQ=255)
                                                                                                                             tCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATCCGTGGTTATGGCGTGCCACAAGTCGTATTTg       >  2:206545/1‑140 (MQ=255)
                                                                                                                                  cGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt  <  2:21739/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
TCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTT  >  NC_000913/3001159‑3001427

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: