Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,022,455 G→A E380E (GAG→GAA xdhD → fused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,022,4550GA97.3% 126.7 / NA 37E380E (GAG→GAAxdhDfused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (15/21);  total (15/22)
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

CAAAAATATCGCCGCCAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAGCGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCTTTACTCGTATGGACGGC  >  NC_000913/3022377‑3022527
                                                                              |                                                                        
cAAAAATATCGCCGCCAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACg                                                                        <  1:523282/81‑1 (MQ=255)
  aaaaTATCGCCGCCAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGaa                                                                      >  2:523316/1‑81 (MQ=255)
    aaTATCGCCGCCAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAAcc                                                                    <  2:523283/81‑1 (MQ=255)
     atatCGCCGCCAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCt                                                                   <  2:523284/81‑1 (MQ=255)
      tatCGCCGCCAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCta                                                                  <  1:523285/81‑1 (MQ=255)
        tCGCCGCCAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAGCGAACCtat                                                                 <  2:523286/80‑1 (MQ=255)
        tCGCCGCCAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCtata                                                                <  1:523288/81‑1 (MQ=255)
        tCGCCGCCAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCtat                                                                 <  1:523287/80‑1 (MQ=255)
           ccgccAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACt                                                             <  2:523289/81‑1 (MQ=255)
              ccAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCaa                                                          <  2:523290/81‑1 (MQ=255)
              ccAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCaa                                                          >  1:523337/1‑81 (MQ=255)
                aGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACg                                                        <  1:523291/81‑1 (MQ=255)
                      cATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAgc                                                   <  1:523292/80‑1 (MQ=255)
                      cATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAgc                                                   <  1:523293/80‑1 (MQ=255)
                       aTGGTCATATTGGCGATATGTACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAgcg                                                  >  2:523306/1‑80 (MQ=255)
                          gTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCgcag                                               >  2:523336/1‑80 (MQ=255)
                               atTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGc                                         >  2:523319/1‑81 (MQ=255)
                                  ggCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGccc                                       >  1:523325/1‑80 (MQ=255)
                                   gCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCg                                      <  2:523294/80‑1 (MQ=255)
                                        atGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGa                                 >  2:523328/1‑80 (MQ=255)
                                          ggACAAAGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGaaa                               >  1:523341/1‑80 (MQ=255)
                                               aaGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACata                         <  1:523295/81‑1 (MQ=255)
                                                aGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACatat                        <  1:523297/81‑1 (MQ=255)
                                                aGGCTTTGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACatat                        <  1:523296/81‑1 (MQ=255)
                                                   cTTTGCCGATGACGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTg                     >  2:523348/1‑81 (MQ=255)
                                                    tttGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGc                    >  1:523331/1‑81 (MQ=255)
                                                     ttGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCt                   >  2:523335/1‑81 (MQ=255)
                                                     ttGCCGATGCCGATGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCt                   >  2:523347/1‑81 (MQ=255)
                                                      tgccgatgccgaTGTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCt                   <  2:523298/80‑1 (MQ=255)
                                                        ccgatgccgatgTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCttt                 <  2:523299/80‑1 (MQ=255)
                                                        ccgatgccgatgTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCttt                 <  1:523300/80‑1 (MQ=255)
                                                          gatgccgatgTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCTTTACt              >  1:523333/1‑81 (MQ=255)
                                                            tgccgatgTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATTTCTGCTTTACTCg            <  2:523301/81‑1 (MQ=255)
                                                            tgccgatgTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCTTTACTCg            >  1:523343/1‑81 (MQ=255)
                                                                 atgTGATCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCTTTACTCGTATgg       <  2:523302/81‑1 (MQ=255)
                                                                      aTCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCTTTACTCGTATGGACGGc  >  1:523361/1‑81 (MQ=255)
                                                                      aTCATTGAACGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCTTTACTCGTATGGACGGc  <  1:523303/81‑1 (MQ=255)
                                                                              |                                                                        
CAAAAATATCGCCGCCAGTATTCATGGTCATATTGGCGATATGGACAAAGGCTTTGCCGATGCCGATGTGATCATTGAGCGAACCTATAACTCAACGCAAGCGCAGCAGTGCCCGACTGAAACACATATCTGCTTTACTCGTATGGACGGC  >  NC_000913/3022377‑3022527

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