Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,361,393 C→T R92C (CGT→TGT)  puuD → gamma‑glutamyl‑gamma‑aminobutyrate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,361,3930CT97.0% 101.6 / NA 33R92C (CGT→TGT) puuDgamma‑glutamyl‑gamma‑aminobutyrate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/15);  minor base A (0/1);  total (17/16)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGA  >  NC_000913/1361257‑1361529
                                                                                                                                        |                                                                                                                                        
acaTGCGCTGGCGGAACCGTCATTACTTGAATAACTTTTACCGAAACTGGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGtgtg                                                                                                                                       <  2:342457/140‑1 (MQ=255)
acaTGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAATGGCGATGAGCCTGACGCCGATCCCGGGtgtg                                                                                                                                       <  2:790716/140‑1 (MQ=255)
acaTGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGtgtg                                                                                                                                       <  1:196429/140‑1 (MQ=255)
acaTGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGtgtg                                                                                                                                       <  1:810097/140‑1 (MQ=255)
        tGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTg                                                                                                                               >  1:492290/1‑140 (MQ=255)
          gCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGCCGATGAGGCTGACGCCTATCCAGGGTGTGATCTTCTGAg                                                                                                                             >  2:259445/1‑140 (MQ=255)
              aaCCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATg                                                                                                                         <  2:492290/140‑1 (MQ=255)
                 cGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCg                                                                                                                      <  2:569041/140‑1 (MQ=255)
                 cGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCg                                                                                                                      <  1:598430/140‑1 (MQ=255)
                                       gCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGaaa                                                                                                >  1:581883/1‑140 (MQ=255)
                                                                   ggTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCg                                                                    <  1:473057/140‑1 (MQ=255)
                                                                       gTCCCAGCAATGTGCAGCCGCACCTTTATTGTGAAAAGGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGt                                                                <  1:259445/140‑1 (MQ=255)
                                                                        tCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCTATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGtt                                                               >  2:775237/1‑140 (MQ=255)
                                                                          ccAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTa                                                             >  2:643949/1‑140 (MQ=255)
                                                                               aaTGTGCAGCCGCACCTATATGCTGAAAACGGCGATGAGCCTGACGCCCATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGa                                                        <  2:333104/140‑1 (MQ=255)
                                                                                aTGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGaa                                                       <  1:744267/140‑1 (MQ=255)
                                                                                   tgCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACtg                                                    >  2:197260/1‑140 (MQ=255)
                                                                                   tgCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCAACCCCATTATCGTCATCTGCCGAGGTTTACAAGAACtg                                                    >  1:322712/1‑140 (MQ=255)
                                                                                     cAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACtggt                                                  >  2:465630/1‑140 (MQ=255)
                                                                                      aGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCTATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACtggtg                                                 <  1:775237/140‑1 (MQ=255)
                                                                                                   tGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCggtg                                    <  2:581883/140‑1 (MQ=255)
                                                                                                      tGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGtt                                 >  2:724671/1‑140 (MQ=255)
                                                                                                      tGAAAACGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGAGGtt                                 >  1:627424/1‑140 (MQ=255)
                                                                                                        aaaaCGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTGGCCATCTGGCGGAGTTAACAAGAACAGGTGGTGGCAACCGGTGGTTCg                               >  2:227931/1‑140 (MQ=255)
                                                                                                        aaaaCGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTATAAGCACTGGTGGTAGCAACCGGTGGTTCg                               >  1:271030/1‑140 (MQ=255)
                                                                                                        aaaaCGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCgg                                      <  1:774792/133‑1 (MQ=255)
                                                                                                        aaaaCGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCgg                                      >  2:774792/1‑133 (MQ=255)
                                                                                                        aaaaCGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCgg                                      >  2:96010/1‑133 (MQ=255)
                                                                                                        aaaaCGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCgg                                      <  1:96010/133‑1 (MQ=255)
                                                                                                        aaaaCGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCg                               >  2:611115/1‑140 (MQ=255)
                                                                                                        aaaaCGGCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCg                               >  2:787033/1‑140 (MQ=255)
                                                                                                             ggCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCa                          >  2:38076/1‑140 (MQ=255)
                                                                                                              gCGATGAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCAt                         <  1:643949/140‑1 (MQ=255)
                                                                                                                  tGAGCCTGACGCCGACCCCGGGAGTTATCCTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCa                     <  1:227931/140‑1 (MQ=255)
                                                                                                                   gAGCCTGACGCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCaa                    <  2:122799/140‑1 (MQ=255)
                                                                                                                            gCCGATCCCGGGTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGa           <  1:611115/140‑1 (MQ=255)
                                                                                                                                     gggTGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGa  <  1:38076/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
ACATGCGCTGGCGGAACCGTCATTACTTGAACAACTTTTGCCGAAACTCGATGGCATTTATCTTCCTGGTAGTCCCAGCAATGTGCAGCCGCACCTATATGGTGAAAACGGCGATGAGCCTGACGCCGATCCCGGGCGTGATCTTCTGAGCATGGCGATAATTAATGCCGCACTCGAAAGGCGCATCCCCATTTTCGCCATCTGCCGGGGTTTACAAGAACTGGTGGTGGCAACCGGTGGTTCGTTGCATCGCAAGCTGTGCGAACAGCCTGA  >  NC_000913/1361257‑1361529

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: