Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,000,899 A→G E185E (GAA→GAG xdhA → xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,000,8990AG98.0% 176.5 / ‑4.0 50E185E (GAA→GAGxdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (24/25);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

AGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAACCATTTG  >  NC_000913/3000752‑3001041
                                                                                                                                                   |                                                                                                                                              
agaagaCGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGc                                                                                                                                             <  1:505040/151‑1 (MQ=255)
       gCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACAt                                                                                                                                      >  1:908054/1‑151 (MQ=255)
         tGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACa                                                                                                                                       >  1:416793/1‑148 (MQ=255)
          gCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATc                                                                                                                                     <  1:728453/149‑1 (MQ=255)
            aCCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTg                                                                                                                                 >  1:618397/1‑151 (MQ=255)
              cAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGc                                                                                                                               >  1:761980/1‑151 (MQ=255)
                 tCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGtg                                                                                                                            >  1:616656/1‑151 (MQ=255)
                  ccATAACGGTGGCAATTTACTGAAACAAAGCACGATTTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGtgg                                                                                                                           <  1:681089/151‑1 (MQ=255)
                      aaCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGtggatgg                                                                                                                       <  1:190911/151‑1 (MQ=255)
                       aCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGtggatgg                                                                                                                       >  1:611605/1‑150 (MQ=255)
                         ggTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATggagg                                                                                                                    >  1:234589/1‑151 (MQ=255)
                            ggCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATg                                                                                                                 <  2:676819/151‑1 (MQ=255)
                            ggCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATg                                                                                                                 <  2:769883/151‑1 (MQ=255)
                                  ttACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTcg                                                                                                            <  1:34773/150‑1 (MQ=255)
                                  ttACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTca                                                                                                            <  2:631387/150‑2 (MQ=255)
                                   tACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTcgc                                                                                                           <  1:483041/150‑1 (MQ=255)
                                       gAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAAtt                                                                                                      <  1:710397/151‑1 (MQ=255)
                                       gAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAAt                                                                                                       <  2:783518/150‑1 (MQ=255)
                                        aaaCCAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTa                                                                                                     <  1:489664/151‑1 (MQ=255)
                                         aaCAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTAc                                                                                                    <  2:640625/151‑1 (MQ=255)
                                                 aCGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGtt                                                                                             >  1:537225/1‑150 (MQ=255)
                                                      gTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGc                                                                                       >  1:374285/1‑151 (MQ=255)
                                                           cGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCa                                                                                  <  1:334940/151‑1 (MQ=255)
                                                            gggTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCAccc                                                                                   >  1:45718/1‑149 (MQ=255)
                                                             ggTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGa                                                                                <  2:143213/151‑1 (MQ=255)
                                                              gTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGAt                                                                               <  1:597248/151‑1 (MQ=255)
                                                                 aTGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATccc                                                                            <  2:548256/151‑1 (MQ=255)
                                                                      caacaaACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGcaca                                                                       <  2:175026/151‑1 (MQ=255)
                                                                           aacaaTCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATTCCGCACATTGtt                                                                  <  1:586327/151‑1 (MQ=255)
                                                                            acaaTCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTc                                                                 >  2:544871/1‑151 (MQ=255)
                                                                             caaTCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTcg                                                                >  1:710364/1‑151 (MQ=255)
                                                                                          cgACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGtggttggt                                                   >  1:256642/1‑151 (MQ=255)
                                                                                              tACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAgg                                               >  1:560579/1‑151 (MQ=255)
                                                                                                   ggTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTg                                          >  2:842594/1‑151 (MQ=255)
                                                                                                    gTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTg                                          >  1:595714/1‑150 (MQ=255)
                                                                                                     tACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTgg                                         >  1:115696/1‑150 (MQ=255)
                                                                                                         ggggCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATAtt                                    >  1:561020/1‑151 (MQ=255)
                                                                                                              aCTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTg                               >  1:73002/1‑151 (MQ=255)
                                                                                                               cTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTg                               >  2:766799/1‑150 (MQ=255)
                                                                                                                 aTCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTc                            <  1:596508/151‑1 (MQ=255)
                                                                                                                   cAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCa                           >  2:732743/1‑150 (MQ=255)
                                                                                                                        ccccGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGt                      >  2:94401/1‑150 (MQ=255)
                                                                                                                         cccGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTAc                    >  1:749/1‑151 (MQ=255)
                                                                                                                          ccGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACg                   <  1:357367/151‑1 (MQ=255)
                                                                                                                           cGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTGCGa                  <  2:82376/151‑1 (MQ=255)
                                                                                                                                     cATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAAc        <  1:842661/151‑1 (MQ=255)
                                                                                                                                     cATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAAc        >  1:269927/1‑151 (MQ=255)
                                                                                                                                      aTTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGa                  <  2:383338/140‑1 (MQ=255)
                                                                                                                                      aTTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGa                  >  1:383364/1‑140 (MQ=255)
                                                                                                                                           cacaTGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAACCATTTg  <  1:666014/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAACCATTTG  >  NC_000913/3000752‑3001041

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