Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,072,718 A→G V49A (GTC→GCC)  rutB ← ureidoacrylate amidohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,072,7180AG96.9% 97.2 / NA 32V49A (GTC→GCC) rutBureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (14/17);  minor base T (0/1);  total (14/18)
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

GGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCG  >  NC_000913/1072581‑1072859
                                                                                                                                         |                                                                                                                                             
gggCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTg                                                                                                                                    <  2:16444/149‑1 (MQ=255)
gggCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGt                                                                                                                                  <  1:299259/151‑1 (MQ=255)
 ggCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTGGCAATGGCCcggggggggg                                                                                                                                   >  2:412199/1‑142 (MQ=255)
     ttCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAAGTTGGCAATGGCCgggggggggg                                                                                                                                   >  2:252681/1‑138 (MQ=255)
        gATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACa                                                                                                                           <  2:805133/150‑1 (MQ=255)
             aTGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAAc                                                                                                                     <  2:519172/151‑1 (MQ=255)
                   aTTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTTATTTTTGCAATGGCCGGGGGGGTGGGTGGGAAATCAAAacc                                                                                                                   >  1:746335/1‑147 (MQ=255)
                      cGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCgag                                                                                                             <  2:410914/150‑1 (MQ=255)
                                       cAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTTGCCATGGCCGGGCGAGTGGTTGAGACATCAAACCCCGCGAGATCTAAGTAgccgcc                                                                                              >  2:190778/1‑148 (MQ=255)
                                        aGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATTTTGGCAATGGGCGGGGGGGTGGGTGAGACATCAAACCCCGCGAGATCTAAGTAGCCGCCggg                                                                                           >  2:783637/1‑150 (MQ=255)
                                                cATGCTGTTCACCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATa                                                                                  <  1:56824/151‑1 (MQ=255)
                                                cATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGGATGTTGGCAATGGCCGGGCGCGGGGGTGGGACATCAAACCCCGCGGGGGCTAAATAGGCGCCTgggggg                                                                                       >  1:149807/1‑142 (MQ=255)
                                                           tCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAAAGGTGGCAACGGCGGGGGGGGAGGGTGAGGCATCAAACCCAGGGAGATCTTAGGAGGCGCCCTGCGGGGCATAAGGGGTTTg                                                                         >  1:764627/1‑149 (MQ=255)
                                                              cAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCatat                                                                    <  2:356236/151‑1 (MQ=255)
                                                                 ccATTTTGAAACAAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTg                                                                         <  2:363784/143‑1 (MQ=255)
                                                                   aTTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATccacc                                                               <  1:783581/151‑1 (MQ=255)
                                                                       tGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGGGGGTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCCGGCGTGGCATAAGCGTTTTTCATATccacc                                                               >  2:461177/1‑147 (MQ=255)
                                                                             cAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGGGGGGGGTTGAGACATCAAAACCGGCGAGATCTAAGTAACCGCCTGGCGTGGGATAAGCGGTTTGCATATccacc                                                               >  2:730024/1‑141 (MQ=255)
                                                                              aGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAgcgc                                                    <  2:559844/151‑1 (MQ=255)
                                                                              aGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAgcg                                                     <  1:461141/150‑1 (MQ=255)
                                                                                atgatCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTGGCAATGGGCGGGCGGGGGGGTGGGACATCAAAACCCGCGGGATCTAAGGAGCCGCCCGGCGGGGGATAAAGGGTTTTCATATCCACCACGATCAgcgccc                                                  >  2:535024/1‑149 (MQ=255)
                                                                                   atCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACttt                                               <  2:577385/151‑1 (MQ=255)
                                                                                     cAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTtg                                              <  1:815315/150‑1 (MQ=255)
                                                                                      aGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTTATGTTGGCCAAGGCCGGGCGGGGGGGTGAGACATCAAACCCGGCGGGATCTAAGTAGGCGCCTGGCGGGGCATAAGAGGTTTGCATATCCACCACGATCAGGGCACttt                                               >  2:181436/1‑148 (MQ=255)
                                                                                           cccTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGa                                       <  1:343830/151‑1 (MQ=255)
                                                                                              tGCCGCTCGCGCTGCGGTCACGGCGGTTTGGATGTTGGCAATGGCCGGGCGAGCGGTTGAGACACCAAACCCGGCGAGGTCTAAGTAGCCGCCTTGCGTGGCCTAAGGGTTTTTCATATCCACCACGATCAGCGCCCTTTTCTGCGGACCg                                    >  2:566976/1‑151 (MQ=255)
                                                                                                  gCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAgg                                <  2:425227/151‑1 (MQ=255)
                                                                                                       cgcTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAAt                            <  2:193350/150‑1 (MQ=255)
                                                                                                               tCACGGTGGTTTCAATGTTGGCAATGTCCGGGCGAGTGGTTGAGACATCAAACCCGGCCAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCTCCACGACCAGCGCCCTTTGCTGCGGATCGAAGGTAATGTCTTcc                     <  1:295588/149‑1 (MQ=255)
                                                                                                                       gTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCgg           >  1:685606/1‑151 (MQ=255)
                                                                                                                               gTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGgtc   <  2:162563/151‑1 (MQ=255)
                                                                                                                                ttGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGgtcg  >  2:688159/1‑151 (MQ=255)
                                                                                                                                         |                                                                                                                                             
GGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCG  >  NC_000913/1072581‑1072859

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