Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,041,411 G→A G13S (GGT→AGT)  zinT → metal‑binding protein ZinT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,041,4110GA92.9% 40.8 / NA 14G13S (GGT→AGT) zinTmetal‑binding protein ZinT
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (5/8);  minor base C (1/0);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAA  >  NC_000913/2041222‑2041655
                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
ccTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCtt                                                                                                                                                                                         >  1:235454/1‑251 (MQ=255)
                 aatTTATTTTATTTAAATAAGATAGATAAGTAGAACTGAGAAAGCCATGATCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTAGCGATTAGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGAGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACaa                                                                                                                                                                         >  2:385434/1‑250 (MQ=255)
                    ttattttattTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACaaaaa                                                                                                                                                                      <  1:381880/250‑1 (MQ=255)
                           atttCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTaa                                                                                                                                                              <  2:215742/251‑1 (MQ=255)
                                            aaGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGatgatg                                                                                                                                             <  1:231365/251‑1 (MQ=255)
                                            aaGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATGAACTCTGGAGGAAACTGTTTTGGCTATTCGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCGCCTGCCTTTTCGCATGGTCATCACTCACACGTCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGatgatg                                                                                                                                             <  1:21561/251‑1 (MQ=255)
                                                                                                   tatGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTACGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTCACAGAGGTCGACCAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGttt                                                                                       >  1:80617/1‑250 (MQ=255)
                                                                                                                          atcatATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAAc                                                               >  1:212536/1‑251 (MQ=255)
                                                                                                                               ataCATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACtt                                                             >  1:343079/1‑248 (MQ=255)
                                                                                                                               ataCATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTGAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAAGGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAAGCCGTTTATCCTTTACTGCAAAGACGCAAACtt                                                             <  2:343079/248‑1 (MQ=255)
                                                                                                                                   aTTAACTCTGGAGGAAACTGTTTTGGCGGTTCGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCGTGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGTATGGAGTCTGGCAATCCGTTTATCTTTTACTGCAAAGTGGCAAAATTGGCCCCg                                                      >  2:306446/1‑251 (MQ=255)
                                                                                                                                                      gTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCgg                                    <  1:326880/250‑1 (MQ=255)
                                                                                                                                                                       ttACAAACTGGCTGTTGCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGAtaaaactaaa                   <  2:101515/250‑1 (MQ=255)
                                                                                                                                                                                       gCTTTAAGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAATTAAATCATTTGCTGAAATTaa  <  2:235454/251‑1 (MQ=255)
                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
CCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATATAACATTACACATCATATACATTAACTCTGGAGGAAACTGTTTTGGCGATTCGTCTTTACAAACTGGCTGTTGCTTTAGGTGTCTTTATTGTTAGCGCTCCTGCCTTTTCGCATGGTCATCACTCACACGGCAAACCCTTAACAGAGGTCGAACAAAAAGCTGCTAATGGTGTTTTTGATGATGCCAATGTACAAAACCGAACGCTCAGTGACTGGGATGGAGTCTGGCAATCCGTTTATCCTTTACTGCAAAGTGGCAAACTTGACCCCGTCTTTCAGAAGAAAGCGGATGCAGATAAAACTAAAACATTTGCTGAAATTAA  >  NC_000913/2041222‑2041655

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