Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 540,515 C→T intergenic (+7/‑50) allB → / → ybbY allantoinase/putative uracil/xanthine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913540,5150CT100.0% 196.5 / NA 61intergenic (+7/‑50)allB/ybbYallantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/31);  total (30/31)

GCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGG  >  NC_000913/540390‑540640
                                                                                                                             |                                                                                                                             
gCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCg                                                                                                                  >  2:60762/1‑139 (MQ=255)
gCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCg                                                                                                                  >  2:1337448/1‑139 (MQ=255)
   gTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTTTCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                    >  2:1608425/1‑134 (MQ=255)
   gTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTTTCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                    <  1:1608425/134‑1 (MQ=255)
     aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCtt                                                                                                             <  1:391064/139‑1 (MQ=255)
                 cgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCAtt                                                                                                 >  1:519489/1‑139 (MQ=255)
                 cgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCAtt                                                                                                 >  2:146366/1‑139 (MQ=255)
                 cgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCAtt                                                                                                 >  1:1302964/1‑139 (MQ=255)
                   tATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATtct                                                                                               <  1:625615/139‑1 (MQ=255)
                       aCGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCgg                                                                                           <  1:922089/139‑1 (MQ=255)
                        cGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGt                                                                                          <  2:519489/139‑1 (MQ=255)
                         gAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                         <  1:1337448/139‑1 (MQ=255)
                               cATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgtg                                                                                   <  2:1031628/139‑1 (MQ=255)
                                 tCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                        >  1:1357120/1‑132 (MQ=255)
                                 tCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                        <  2:1357120/132‑1 (MQ=255)
                                  cTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgtgtgt                                                                                >  2:883240/1‑139 (MQ=255)
                                         ggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATGCTCCGGTTAAGGTGAGTTTATgct                                                                         >  1:949978/1‑137 (MQ=255)
                                               gtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGc                                                                                                    <  1:1175277/106‑1 (MQ=255)
                                               gtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGc                                                                                                    >  2:1175277/1‑106 (MQ=255)
                                                    ttACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAg                                                              >  1:648555/1‑139 (MQ=255)
                                                          cATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGcc                                                        <  2:648555/139‑1 (MQ=255)
                                                             tGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCcgcg                                                     >  2:657724/1‑139 (MQ=255)
                                                             tGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCcgcg                                                     <  2:400315/139‑1 (MQ=255)
                                                                  aGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                <  1:883240/139‑1 (MQ=255)
                                                                  aGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                <  2:1302964/139‑1 (MQ=255)
                                                                      cTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGt                                            >  1:275338/1‑139 (MQ=255)
                                                                         cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                <  2:938875/132‑1 (MQ=255)
                                                                         cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                         >  1:1251810/1‑139 (MQ=255)
                                                                         cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                         >  2:255659/1‑139 (MQ=255)
                                                                         cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCATGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCACTCAGCCGCGAAAGc                                                >  1:938875/1‑132 (MQ=255)
                                                                            tGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAg                                      <  1:60762/139‑1 (MQ=255)
                                                                              ttGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                        <  2:399592/87‑1 (MQ=255)
                                                                              ttGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                        >  1:399592/1‑87 (MQ=255)
                                                                                   ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                    >  2:1275733/1‑54 (MQ=255)
                                                                                   ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                    <  1:1275733/54‑1 (MQ=255)
                                                                                   ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                         >  2:407818/1‑81 (MQ=255)
                                                                                   ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                         <  1:407818/81‑1 (MQ=255)
                                                                                    cGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAg                              <  1:255659/139‑1 (MQ=255)
                                                                                         ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGtt                         <  2:275338/139‑1 (MQ=255)
                                                                                                    ccTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATtctc                                                                                              <  1:1415463/59‑1 (MQ=255)
                                                                                                    ccTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATtctc                                                                                              >  2:1415463/1‑59 (MQ=255)
                                                                                                    ccTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGa                                                    >  1:749518/1‑101 (MQ=255)
                                                                                                    ccTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGa                                                    <  2:749518/101‑1 (MQ=255)
                                                                                                       tAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCggg                                                                                                     <  1:724982/49‑1 (MQ=255)
                                                                                                       tAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCggg                                                                                                     >  2:724982/1‑49 (MQ=255)
                                                                                                        aaaCATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTg          <  2:1433016/139‑1 (MQ=255)
                                                                                                        aaaCATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTg          >  1:1433016/1‑139 (MQ=255)
                                                                                                        aaaCATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTg          <  2:824612/139‑1 (MQ=255)
                                                                                                         aaCATCAGCAGTAATCCGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                         >  1:1569601/1‑59 (MQ=255)
                                                                                                         aaCATCAGCAGTAATCCGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                         <  2:1569601/59‑1 (MQ=255)
                                                                                                              cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                               <  2:308294/112‑1 (MQ=255)
                                                                                                              cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa                               >  1:308294/1‑112 (MQ=255)
                                                                                                              cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac    >  2:560068/1‑138 (MQ=255)
                                                                                                              cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac    <  1:560068/138‑1 (MQ=255)
                                                                                                               agcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAg                                      >  1:1239961/1‑104 (MQ=255)
                                                                                                               agcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAg                                      <  2:1239961/104‑1 (MQ=255)
                                                                                                                gcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                         >  2:1335337/1‑100 (MQ=255)
                                                                                                                gcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                         <  1:1335337/100‑1 (MQ=255)
                                                                                                                gcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACgg  >  1:783308/1‑139 (MQ=255)
                                                                                                                      aTCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCa                                                                       <  2:948079/64‑1 (MQ=255)
                                                                                                                      aTCTGGCTCCTGCAATGCCCATCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCa                                                                       >  1:948079/1‑64 (MQ=255)
                                                                                                                             |                                                                                                                             
GCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGG  >  NC_000913/540390‑540640

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