Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,995,567 A→G T635A (ACC→GCC)  xdhA → xanthine dehydrogenase molybdenum‑binding subunit XdhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,995,5670AG96.8% 192.1 / ‑2.0 62T635A (ACC→GCC) xdhAxanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (29/31);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAAT  >  NZ_CP009273/2995480‑2995649
                                                                                       |                                                                                  
gAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGcc                                                                                  >  2:488641/1‑90 (MQ=255)
     cTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCaa                                                                             >  1:1562646/1‑90 (MQ=255)
         aTCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGc                                                                         <  2:1024037/90‑1 (MQ=255)
          tCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCa                                                                        <  2:1225139/90‑1 (MQ=255)
               accaccACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTc                                                                   <  2:769553/90‑1 (MQ=255)
                 caccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCaa                                                                 <  2:1467762/90‑1 (MQ=255)
                 caccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCaa                                                                 <  2:811247/90‑1 (MQ=255)
                 caccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCaa                                                                 <  2:833704/90‑1 (MQ=255)
                   ccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACg                                                               >  2:179330/1‑90 (MQ=255)
                   ccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACg                                                               <  1:422274/90‑1 (MQ=255)
                        aaCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          <  2:346046/90‑1 (MQ=255)
                        aaCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          <  1:203295/90‑1 (MQ=255)
                        aaCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          <  1:910876/90‑1 (MQ=255)
                              ccGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCa                                                    <  1:292793/90‑1 (MQ=255)
                                  cGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGc                                                <  1:179330/90‑1 (MQ=255)
                                   gTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCa                                               <  2:611571/90‑1 (MQ=255)
                                    tttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCat                                              <  1:1376534/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  2:867358/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  1:1005660/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  2:126010/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  1:780144/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  1:404050/90‑1 (MQ=255)
                                     ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  2:1177242/90‑1 (MQ=255)
                                           gTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTa                                       >  2:416739/1‑90 (MQ=255)
                                           gTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTa                                       >  2:262343/1‑90 (MQ=255)
                                               cTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATcc                                   <  2:38242/90‑1 (MQ=255)
                                                tttGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCa                                  >  1:146034/1‑90 (MQ=255)
                                                tttGTTGATCTGACGATCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCa                                  >  1:65229/1‑90 (MQ=255)
                                                 ttgttgATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCAc                                 >  2:1430/1‑90 (MQ=255)
                                                 ttgttgATCTGACGGTCGATATTGCGCTGTGCAAAGCCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCAc                                 >  2:590196/1‑90 (MQ=255)
                                                      gATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCActgctg                            >  2:45721/1‑90 (MQ=255)
                                                        tCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          <  2:678386/58‑1 (MQ=255)
                                                        tCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          >  1:678386/1‑58 (MQ=255)
                                                        tCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  1:342914/71‑1 (MQ=255)
                                                        tCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             >  2:342914/1‑71 (MQ=255)
                                                        tCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGc                          >  2:621256/1‑90 (MQ=255)
                                                           gACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGa                       >  2:1383327/1‑90 (MQ=255)
                                                            aCGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGaa                      <  2:270662/90‑1 (MQ=255)
                                                             cGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAg                     >  2:1023541/1‑90 (MQ=255)
                                                              ggTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAgg                    <  1:874985/90‑1 (MQ=255)
                                                                    tatTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGt              <  1:45721/90‑1 (MQ=255)
                                                                        gcgcTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTacac          >  2:283047/1‑90 (MQ=255)
                                                                           gtgtgCCAAGTCGCCATCGACCGCACCCCCAACGTTCATGATTCAGGCCATCTTCTTAATCCACTGCTGGCAGAAGATCAGGTACAcggc       <  1:821688/89‑1 (MQ=255)
                                                                           cTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggc       <  1:748857/90‑1 (MQ=255)
                                                                           cTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggc       <  2:1405450/90‑1 (MQ=255)
                                                                           cTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggc       <  2:333362/90‑1 (MQ=255)
                                                                             gtgCAAAGTCGCCATCGACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  1:1418594/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          <  2:742028/37‑1 (MQ=38)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                          >  1:742028/1‑37 (MQ=38)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             >  1:856412/1‑50 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                             <  2:856412/50‑1 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  2:860500/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  1:108864/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  2:480016/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  2:374905/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  2:239398/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  2:1486175/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  1:590968/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  1:250761/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTCCTGGCAGAAGGTCAGGTACAcggcgg     >  2:108497/1‑90 (MQ=255)
                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACCTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg     >  1:655827/1‑90 (MQ=255)
                                                                                cAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCggaat  <  2:1364884/90‑1 (MQ=255)
                                                                                       |                                                                                  
GAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAAT  >  NZ_CP009273/2995480‑2995649

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