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

GTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGTCATC  >  NC_000913/1072574‑1072864
                                                                                                                                                 |                                                                                                                                                  
gTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGc                                                                                                                                              <  2:356593/151‑1 (MQ=255)
 ttttCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAAT‑GTTGGCAATTGCCggg                                                                                                                                               >  1:502095/1‑149 (MQ=255)
   ttCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAAT‑GTTGGGAATGGGCggggg                                                                                                                                             >  2:719264/1‑147 (MQ=255)
            ttCGATTTATGAAAAGTCGGTGAGCCGGGCCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGTTCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTgag                                                                                                                                   <  1:606387/150‑1 (MQ=255)
             tCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTgaga                                                                                                                                  <  2:403458/150‑1 (MQ=255)
                    aTGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCaaa                                                                                                                           <  1:412233/150‑1 (MQ=255)
                     tGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGGAT‑GTTGGGAATGGCCGGGCGGGGGGGTGAGACATCaaa                                                                                                                           >  1:130147/1‑149 (MQ=255)
                           ttCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTTAT‑TTTTGCAATGGCCCGGCGGGGGGGTGAGGCATCaaa                                                                                                                           >  1:479110/1‑143 (MQ=255)
                           ttCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGGGAGGGGGTGAGACATCAAAACCgg                                                                                                                      >  2:89577/1‑148 (MQ=255)
                                 gAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCt                                                                                                              <  1:398914/150‑1 (MQ=255)
                                  aGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTaa                                                                                                            <  2:461889/151‑1 (MQ=255)
                                         tCCGCCAGCCTCGACATCCTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAgcc                                                                                                      <  1:509887/150‑1 (MQ=255)
                                            gCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAgccgcc                                                                                                   <  1:34735/150‑1 (MQ=255)
                                             ccAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAAT‑TTTGGCAAAGGGCCGGGGGGGGGGTGAGACATCAAACCCCGCGAGATCTAAGTAgccgcc                                                                                                   >  2:131132/1‑149 (MQ=255)
                                             ccAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGATTGAAT‑GTTGGGAAATGGCGGGGGAGAGGTTGTGACAACCAACCCCGCGCGATCTAATTAaccgcc                                                                                                   >  2:5815/1‑149 (MQ=255)
                                              cAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGGAATGGCCGGGCGGGGGGGTGGGACATCAAACCCGGCGAGATCTAAGTAGCCGCCCgg                                                                                                >  1:187168/1‑151 (MQ=255)
                                                       cATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGGGGGTGAGACATCAAACCCGGCGGGATCTAAGTAGCCGCCTGGCGGGGCa                                                                                         >  2:584521/1‑149 (MQ=255)
                                                         tACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGTAT‑GTTGGCAATGGCCGGGCGGGGGGGTGAGACATCAAAACCGGCGGGGTCTAAGTAGCCGCCCGGCGTGGCATaa                                                                                      >  1:618816/1‑150 (MQ=255)
                                                          aCTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGc                                                                                    <  1:751999/151‑1 (MQ=255)
                                                                cATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTTAT‑ATTGGCAATGGCCGGGCGAGGGGTTGAGACATCAAACCCGGCGGGATCTAAATAGCCGCCCGGCGTGGCATAAACGtttt                                                                               >  1:569621/1‑150 (MQ=255)
                                                                cATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGtttt                                                                               <  1:654177/150‑1 (MQ=255)
                                                                        ccATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATcca                                                                      <  2:711994/151‑1 (MQ=255)
                                                                          aTTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATccacc                                                                    <  2:578010/151‑1 (MQ=255)
                                                                                 aaCCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCa                                                             <  1:105291/151‑1 (MQ=255)
                                                                                 aaCCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAAAGGCCGGGCGAGGGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCCGGCGTGGCATAAACGTTTTGCATATTcaccac                                                                  >  2:155672/1‑146 (MQ=255)
                                                                                   ccAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGGATGGTTGGCATTG‑CCGGGCGGGGGGGTGAGACATCAAACCCCGCGGGATCTAAGTAGCCGCCCTGCGTGGCCTAAGCGGTTTTGATATCCACCCCggaca                                                             >  1:562582/1‑145 (MQ=255)
                                                                                   ccAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAgc                                                           <  1:41178/151‑1 (MQ=255)
                                                                                       atgatCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAAT‑GTTGGCAATGGCCGGGCGGGTGGTTGGGACATCCAACCCGGCGGGATCTAAGTAGCCGCCTGGCGTGGGCTAAGCGTTTTGCATATCCACCACGATCAGCGCAc                                                       >  1:137396/1‑151 (MQ=255)
                                                                                       atgatCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTTCATATCCACCACGATCAgcg                                                          >  2:611712/1‑148 (MQ=255)
                                                                                              gCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTtgct                                                 <  2:101921/150‑1 (MQ=255)
                                                                                                aTCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCg                                              <  1:122605/151‑1 (MQ=255)
                                                                                                 tCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTTAA‑GTTTGCACAGGCCGGGCGTGGGGGTGGGACAACAAACCCCGCGAGAATATAGTAGCCGGCTGGCGGGGAATAAAGGGGTTGCATATTCCACACGCTCAGCGCACctt                                                    >  1:610376/1‑144 (MQ=255)
                                                                                                  cccTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCgg                                             <  1:681969/150‑1 (MQ=255)
                                                                                                   ccTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GGTGGGAATGGGCGGGGGAGGGGTTGGGACATAAAAACCCGCGAGATTTAAGAAGGCGCCCGGGGTGGCATAAGAGTTTTGCATATCCCCCACGATTAGCGCACTTTGCTGCgg                                             >  2:214569/1‑149 (MQ=255)
                                                                                                     tGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATTGCCGGGGGGGGGGGTGGGACCTCAAACCCCGCGGGATCTTAGTAGCCGCCCGGCGGGGCCAAAACGGTTTTCCTATCCCCCCCCAACAACGCCCTTTTCTGCGGGTCg                                         >  1:277731/1‑151 (MQ=255)
                                                                                                     tGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAg                                                            <  1:710309/132‑1 (MQ=255)
                                                                                                     tGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGGTGAGACATCAAACCCCGCGAGATCTAAGTAGCCGCCCGGCGTGGCATAAACGGTTTTCATATCCACCACGATCa                                                             >  2:710383/1‑131 (MQ=255)
                                                                                                        cGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAg                                      <  2:594246/151‑1 (MQ=255)
                                                                                                            cgcgCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGt                                    <  2:438251/149‑1 (MQ=255)
                                                                                                                cTGCGGTCACGGCGGTTTGAAT‑GGTGGCAAAGGCCGGGCGAGGGGGTGAGACATCAAACCCCGGGAGATCTAAGGAGCCGCCCGGCGTGGGCTAAGCGTTTTGCATATCCACCACGACAAGCGCACCTTGCTGCGGGGCGGAGGGAATgg                               >  1:78037/1‑150 (MQ=255)
                                                                                                                                  gAAT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTa             >  2:723294/1‑150 (MQ=255)
                                                                                                                                   aaT‑GTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTa             >  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)
                                                                                                                                                 |                                                                                                                                                  
GTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAAT‑GTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGTCATC  >  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: ‑