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_000913	540,515	0	C	T	100.0%	241.5 / NA	74	intergenic (+7/‑50)	allB/ybbY	allantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/38);  total (36/38)

ATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGT  >  NC_000913/540384‑540641                                                                                                                                    |                                                                                                                               aTGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCaa                                                                                                                         <  1:143437/139‑1 (MQ=255)      ggCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGccc                                                                                                                    >  1:1314951/1‑139 (MQ=255)           tACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTCATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCt                                                                                                               >  2:400838/1‑139 (MQ=255)            aCCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCtt                                                                                                              >  1:446462/1‑139 (MQ=255)                      cgcgTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCa                                                                                                    <  1:1568468/139‑1 (MQ=255)                               cGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGt                                                                                           <  2:1693969/139‑1 (MQ=255)                                gAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                          <  2:446462/139‑1 (MQ=255)                                 aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                         <  2:437385/139‑1 (MQ=255)                                 aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTa                                                                                         <  2:1314951/139‑1 (MQ=255)                                   aaCCATCTTACGTGGTGATGTGATTTACGACATTGAGCAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAg                                                                                       >  2:1466289/1‑139 (MQ=255)                                        tCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgtgtg                                                                                  >  1:1511046/1‑139 (MQ=255)                                             cGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATACTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTAt                                                                             >  2:132172/1‑139 (MQ=255)                                                 gtgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgt                                                                                     >  2:1676236/1‑127 (MQ=255)                                                 gtgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgt                                                                                     <  1:1676236/127‑1 (MQ=255)                                                  tgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCa                                                                        >  1:582569/1‑139 (MQ=255)                                                    atgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTg                                                                                                                            <  2:830476/85‑1 (MQ=255)                                                    atgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTg                                                                                                                            >  1:830476/1‑85 (MQ=255)                                                    atgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAAt                                                                      >  2:1854329/1‑139 (MQ=255)                                                           ttACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAg                                                               >  1:1533706/1‑139 (MQ=255)                                                            tACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTAATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGt                                                              >  1:1332560/1‑139 (MQ=255)                                                              cGACATTGAGCAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCa                                                            <  1:1466289/139‑1 (MQ=255)                                                               gACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAg                                                           <  1:1141700/139‑1 (MQ=255)                                                                    tGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGc                                                                                                                           >  2:398854/1‑70 (MQ=255)                                                                    tGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGc                                                                                                                           <  1:398854/70‑1 (MQ=255)                                                                      aaCAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGt                                                                                           >  1:1571233/1‑100 (MQ=255)                                                                      aaCAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGt                                                                                           <  2:1571233/100‑1 (MQ=255)                                                                         aGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                 <  2:1332560/139‑1 (MQ=255)                                                                           ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCggg                                                                                                      >  1:187238/1‑84 (MQ=255)                                                                           ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCggg                                                                                                      <  2:187238/84‑1 (MQ=255)                                                                             cTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCa                                                                                                                          >  1:808568/1‑62 (MQ=255)                                                                             cTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCa                                                                                                                          <  2:808568/62‑1 (MQ=255)                                                                             cTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAg                                                                                       <  2:1407013/97‑1 (MQ=255)                                                                             cTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAg                                                                                       >  1:1407013/1‑97 (MQ=255)                                                                              ttCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                     >  2:473406/1‑66 (MQ=255)                                                                              ttCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                     <  1:473406/66‑1 (MQ=255)                                                                              ttCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGtt                                            <  1:504620/139‑1 (MQ=255)                                                                                cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                          <  2:1511046/139‑1 (MQ=255)                                                                                 ccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATc                                         <  1:1854329/139‑1 (MQ=255)                                                                                     ttGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGc                                                                                                                           <  1:550142/53‑1 (MQ=255)                                                                                     ttGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGc                                                                                                                           >  2:550142/1‑53 (MQ=255)                                                                                       gcgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGcggcgg                                                                                                       >  1:1246528/1‑71 (MQ=255)                                                                                       gcgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGcggcgg                                                                                                       <  2:1246528/71‑1 (MQ=255)                                                                                        cgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                 <  2:1463380/124‑1 (MQ=255)                                                                                        cgcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                 >  1:1463380/1‑124 (MQ=255)                                                                                          ccGAAAGGTCATTTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTgc                                                                                                            <  2:1438993/63‑1 (MQ=255)                                                                                          ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                     >  2:1711342/1‑54 (MQ=255)                                                                                          ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGcc                                                                                                                     <  1:1711342/54‑1 (MQ=255)                                                                                          ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTgc                                                                                                            >  1:1438993/1‑63 (MQ=255)                                                                                          ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgtgtgt                                                                                 <  2:1103743/90‑1 (MQ=255)                                                                                          ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgtgtgt                                                                                 >  1:1103743/1‑90 (MQ=255)                                                                                             aaaGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTc                                                             <  1:1025012/107‑1 (MQ=255)                                                                                             aaaGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTc                                                             >  2:1025012/1‑107 (MQ=255)                                                                                               aGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGt                           >  2:469168/1‑139 (MQ=255)                                                                                                ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGtt                          <  1:469168/139‑1 (MQ=255)                                                                                                           ccTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAttt               >  2:1740056/1‑139 (MQ=255)                                                                                                            cTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCa                                                            >  2:664770/1‑93 (MQ=255)                                                                                                            cTTAAACATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCa                                                            <  1:664770/93‑1 (MQ=255)                                                                                                             ttAAACATCTGCAGTAATCTGGCTCCTGCGATGACGGTCCTTGCGGGGGGCATTTTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAttttt             <  1:132172/139‑1 (MQ=255)                                                                                                               aaaCATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCa                                                                                                    <  2:1547420/50‑1 (MQ=255)                                                                                                               aaaCATCAGCAGTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCa                                                                                                    >  1:1547420/1‑50 (MQ=255)                                                                                                                     cagcagTAATCTGGCTCCTGTAATGCCCGTCCTTGCGGCGGGc                                                                                                     <  1:1441250/43‑1 (MQ=255)                                                                                                                     cagcagTAATCTGGCTCCTGTAATGCCCGTCCTTGCGGCGGGc                                                                                                     >  2:1441250/1‑43 (MQ=255)                                                                                                                     cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTgcgg                                                                                                          <  2:1022513/38‑1 (MQ=255)                                                                                                                     cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTgcgg                                                                                                          >  1:1022513/1‑38 (MQ=255)                                                                                                                     cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                          <  2:394069/54‑1 (MQ=255)                                                                                                                     cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGtt                                                                                          >  1:394069/1‑54 (MQ=255)                                                                                                                     cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                 >  2:1796114/1‑95 (MQ=255)                                                                                                                     cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGc                                                 <  1:1796114/95‑1 (MQ=255)                                                                                                                     cagcagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGGTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac     >  1:1415521/1‑139 (MQ=255)                                                                                                                        cagTAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGt  <  2:582569/139‑1 (MQ=255)                                                                                                                           tAATCTGGCTCCTGCAATGGCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGtt                          >  2:1253344/1‑112 (MQ=255)                                                                                                                           tAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATtctc                                                                                               >  1:796330/1‑43 (MQ=255)                                                                                                                           tAATCTGGCTCCTGCAATGCCCGTCCTTGCGGCGGGCATtctc                                                                                               <  2:796330/43‑1 (MQ=255)                                                                                                                           tAATCTGGCTCCTGCAATGCCCGTACTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGtt                          <  1:1253344/112‑1 (MQ=255)                                                                                                                                    |                                                                                                                               ATGTTGGCCGTACCATTGGCGCGCGTATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGT  >  NC_000913/540384‑540641

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