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,7180AG97.8% 141.2 / ‑2.7 45coding (146/693 nt)rutBureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/22);  minor base . (1/0);  total (23/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 = 8.85e-01

GTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGTCATC  >  NC_000913/1072574‑1072864
                                                                                                                                                 |                                                                                                                                                  
gTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGc                                                                                                                                              <  2:356593/151‑1 (MQ=255)
 ttttCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTGGCAATTGCCggg                                                                                                                                               >  1:502095/1‑149 (MQ=255)
   ttCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTGGGAATGGGCggggg                                                                                                                                             >  2:719264/1‑147 (MQ=255)
            ttCGATTTATGAAAAGTCGGTGAGCCGGGCCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGTTCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTgag                                                                                                                                   <  1:606387/150‑1 (MQ=255)
             tCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTgaga                                                                                                                                  <  2:403458/150‑1 (MQ=255)
                    aTGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCaaa                                                                                                                           <  1:412233/150‑1 (MQ=255)
                     tGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGGATGTTGGGAATGGCCGGGCGGGGGGGTGAGACATCaaa                                                                                                                           >  1:130147/1‑149 (MQ=255)
                           ttCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTTATTTTTGCAATGGCCCGGCGGGGGGGTGAGGCATCaaa                                                                                                                           >  1:479110/1‑143 (MQ=255)
                           ttCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGGGAGGGGGTGAGACATCAAAACCgg                                                                                                                      >  2:89577/1‑148 (MQ=255)
                                 gAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCt                                                                                                              <  1:398914/150‑1 (MQ=255)
                                  aGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTaa                                                                                                            <  2:461889/151‑1 (MQ=255)
                                         tCCGCCAGCCTCGACATCCTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAgcc                                                                                                      <  1:509887/150‑1 (MQ=255)
                                            gCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAgccgcc                                                                                                   <  1:34735/150‑1 (MQ=255)
                                             ccAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATTTTGGCAAAGGGCCGGGGGGGGGGTGAGACATCAAACCCCGCGAGATCTAAGTAgccgcc                                                                                                   >  2:131132/1‑149 (MQ=255)
                                             ccAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGATTGAATGTTGGGAAATGGCGGGGGAGAGGTTGTGACAACCAACCCCGCGCGATCTAATTAaccgcc                                                                                                   >  2:5815/1‑149 (MQ=255)
                                              cAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGGAATGGCCGGGCGGGGGGGTGGGACATCAAACCCGGCGAGATCTAAGTAGCCGCCCgg                                                                                                >  1:187168/1‑151 (MQ=255)
                                                       cATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGGGGGTGAGACATCAAACCCGGCGGGATCTAAGTAGCCGCCTGGCGGGGCa                                                                                         >  2:584521/1‑149 (MQ=255)
                                                         tACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGTATGTTGGCAATGGCCGGGCGGGGGGGTGAGACATCAAAACCGGCGGGGTCTAAGTAGCCGCCCGGCGTGGCATaa                                                                                      >  1:618816/1‑150 (MQ=255)
                                                          aCTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGc                                                                                    <  1:751999/151‑1 (MQ=255)
                                                                cATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTTATATTGGCAATGGCCGGGCGAGGGGTTGAGACATCAAACCCGGCGGGATCTAAATAGCCGCCCGGCGTGGCATAAACGtttt                                                                               >  1:569621/1‑150 (MQ=255)
                                                                cATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGtttt                                                                               <  1:654177/150‑1 (MQ=255)
                                                                        ccATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATcca                                                                      <  2:711994/151‑1 (MQ=255)
                                                                          aTTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATccacc                                                                    <  2:578010/151‑1 (MQ=255)
                                                                                 aaCCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCa                                                             <  1:105291/151‑1 (MQ=255)
                                                                                 aaCCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAAAGGCCGGGCGAGGGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCCGGCGTGGCATAAACGTTTTGCATATTcaccac                                                                  >  2:155672/1‑146 (MQ=255)
                                                                                   ccAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGGATGGTTGGCATTGCCGGGCGGGGGGGTGAGACATCAAACCCCGCGGGATCTAAGTAGCCGCCCTGCGTGGCCTAAGCGGTTTTGATATCCACCCCggaca                                                             >  1:562582/1‑145 (MQ=255)
                                                                                   ccAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAgc                                                           <  1:41178/151‑1 (MQ=255)
                                                                                       atgatCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTGGCAATGGCCGGGCGGGTGGTTGGGACATCCAACCCGGCGGGATCTAAGTAGCCGCCTGGCGTGGGCTAAGCGTTTTGCATATCCACCACGATCAGCGCAc                                                       >  1:137396/1‑151 (MQ=255)
                                                                                       atgatCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTTCATATCCACCACGATCAgcg                                                          >  2:611712/1‑148 (MQ=255)
                                                                                              gCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTtgct                                                 <  2:101921/150‑1 (MQ=255)
                                                                                                aTCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCg                                              <  1:122605/151‑1 (MQ=255)
                                                                                                 tCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTTAAGTTTGCACAGGCCGGGCGTGGGGGTGGGACAACAAACCCCGCGAGAATATAGTAGCCGGCTGGCGGGGAATAAAGGGGTTGCATATTCCACACGCTCAGCGCACctt                                                    >  1:610376/1‑144 (MQ=255)
                                                                                                  cccTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCgg                                             <  1:681969/150‑1 (MQ=255)
                                                                                                   ccTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGGTGGGAATGGGCGGGGGAGGGGTTGGGACATAAAAACCCGCGAGATTTAAGAAGGCGCCCGGGGTGGCATAAGAGTTTTGCATATCCCCCACGATTAGCGCACTTTGCTGCgg                                             >  2:214569/1‑149 (MQ=255)
                                                                                                     tGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATTGCCGGGGGGGGGGGTGGGACCTCAAACCCCGCGGGATCTTAGTAGCCGCCCGGCGGGGCCAAAACGGTTTTCCTATCCCCCCCCAACAACGCCCTTTTCTGCGGGTCg                                         >  1:277731/1‑151 (MQ=255)
                                                                                                     tGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAg                                                            <  1:710309/132‑1 (MQ=255)
                                                                                                     tGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGGTGAGACATCAAACCCCGCGAGATCTAAGTAGCCGCCCGGCGTGGCATAAACGGTTTTCATATCCACCACGATCa                                                             >  2:710383/1‑131 (MQ=255)
                                                                                                        cGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAg                                      <  2:594246/151‑1 (MQ=255)
                                                                                                            cgcgCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGt                                    <  2:438251/149‑1 (MQ=255)
                                                                                                                cTGCGGTCACGGCGGTTTGAATGGTGGCAAAGGCCGGGCGAGGGGGTGAGACATCAAACCCCGGGAGATCTAAGGAGCCGCCCGGCGTGGGCTAAGCGTTTTGCATATCCACCACGACAAGCGCACCTTGCTGCGGGGCGGAGGGAATgg                               >  1:78037/1‑150 (MQ=255)
                                                                                                                                  gAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTa             >  2:723294/1‑150 (MQ=255)
                                                                                                                                   aaTGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTa             >  1:520093/1‑149 (MQ=255)
                                                                                                                                       gTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGgtc        <  1:453641/151‑1 (MQ=255)
                                                                                                                                            cAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGtcat   <  2:499379/151‑1 (MQ=255)
                                                                                                                                             aaTGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGtcatc  >  1:166026/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGTCATC  >  NC_000913/1072574‑1072864

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