Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,478,527	Δ1 bp	intergenic (+36/+23)	yjgN → / ← valS	hypothetical protein/valine tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,478,527	0	G	.	100.0%	404.3 / NA	88	intergenic (+36/+23)	yjgN/valS	hypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (53/35);  total (53/35)

TTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGC  >  NC_000913/4478454‑4478597                                                                          |                                                                       tttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGttt                                                                   >  1:707463/1‑76 (MQ=255)  ttCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGtttt                                                                  >  1:2343422/1‑76 (MQ=255)   tCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTtc                                                                 >  2:2033553/1‑76 (MQ=255)    cACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTtca                                                                <  2:3636799/76‑1 (MQ=255)     aCGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAACGCCGG‑CGCATGCTCCGGCCTT‑CGTTTtcac                                                               >  2:1519463/1‑75 (MQ=255)      cGGGGAGTGATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTtcatc                                                              <  1:4388812/76‑1 (MQ=255)       ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTtcatcc                                                             >  2:56451/1‑75 (MQ=255)       ggggAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTtcatca                                                             <  2:985894/76‑1 (MQ=255)        gggAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCAc                                                            <  1:1948184/76‑1 (MQ=255)         ggAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACt                                                           <  1:548008/76‑1 (MQ=255)          gAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACtg                                                          <  1:1256798/76‑1 (MQ=255)          gAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACtg                                                          >  2:4741804/1‑76 (MQ=255)          gAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACtg                                                          <  1:4321047/76‑1 (MQ=255)           aGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACtgt                                                         >  1:842519/1‑76 (MQ=255)            gTTAGGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAGACGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACtgtg                                                        <  1:3970872/76‑1 (MQ=255)             ttATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTgtgt                                                       >  2:2329213/1‑76 (MQ=255)              tATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGtt                                                      >  2:422308/1‑76 (MQ=255)              tATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGtt                                                      <  2:2912015/76‑1 (MQ=255)               aTGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGttt                                                     <  1:2631366/76‑1 (MQ=255)                 gCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTg                                                   <  2:4405138/76‑1 (MQ=255)                 gCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTg                                                   >  2:3120087/1‑76 (MQ=255)                  ccTTCTTTCCCATTTCTGTAAGCGCT‑AAAAACGGCCGGA‑GCATGCTCCGGCATT‑CGTTTTCATCACTGTGtttttt                                                  >  2:1315322/1‑74 (MQ=255)                  ccTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGa                                                  >  2:4342939/1‑76 (MQ=255)                   cttcttTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGAt                                                 <  2:4163031/76‑1 (MQ=255)                   cttcttTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGAt                                                 >  2:2234281/1‑76 (MQ=255)                    ttcttTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGAtt                                                <  2:3576600/76‑1 (MQ=255)                     tcttTACCATTTCTGTAAGAGAT‑AAAATAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGAttc                                               >  2:1330722/1‑75 (MQ=255)                     tcttTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTa                                               >  2:340851/1‑76 (MQ=255)                     tcttTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCAGGCTCCGGCCTT‑CGTTTTCAGCACTGCGTTTTGCttc                                               >  2:303361/1‑75 (MQ=255)                       tttACCATTTCTGTAAGAGAT‑NAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACa                                             >  1:3276041/1‑76 (MQ=255)                       tttACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACa                                             >  2:4433068/1‑76 (MQ=255)                        ttACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAg                                            <  2:2229050/76‑1 (MQ=255)                         tACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAgc                                           >  2:4703852/1‑76 (MQ=255)                         tACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAgc                                           <  1:733602/76‑1 (MQ=255)                         tACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAgc                                           <  2:367421/76‑1 (MQ=255)                         tACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAgc                                           >  2:2246381/1‑76 (MQ=255)                           ccATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAgcgc                                         <  2:234122/76‑1 (MQ=255)                           ccATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAgcgc                                         <  1:2380398/76‑1 (MQ=255)                            cATTTCTGTAAGAGAT‑AACAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGctcg                                        >  1:1218092/1‑76 (MQ=255)                            cATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGcgcg                                        >  1:1205146/1‑76 (MQ=255)                             aTTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCgcgg                                       >  1:4253217/1‑76 (MQ=255)                              tttCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCgcggc                                      >  2:2894803/1‑76 (MQ=255)                              tttCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCgcggc                                      >  2:2302779/1‑76 (MQ=255)                                tCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGa                                    <  1:419864/76‑1 (MQ=255)                                   gTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                 >  1:831207/1‑76 (MQ=255)                                   gTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                 <  2:603731/76‑1 (MQ=255)                                   gTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                 >  1:1225971/1‑76 (MQ=255)                                   gTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa                                 <  2:4432840/76‑1 (MQ=255)                                    tAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAAc                                >  2:2358754/1‑76 (MQ=255)                                     aaGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACa                               <  1:1866192/76‑1 (MQ=255)                                      agagAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACac                              >  2:2236830/1‑75 (MQ=255)                                          gataaaaaaGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct                           <  1:594893/73‑1 (MQ=255)                                          gataaaaaaGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCct                           <  2:3902864/73‑1 (MQ=255)                                           ataaaaaaGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctg                          >  1:1802882/3‑76 (MQ=255)                                            taaaaaaGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc                         >  1:1912756/2‑76 (MQ=255)                                            taaaaaaGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCctgc                         >  1:2013989/2‑76 (MQ=255)                                               aaaaGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGt                      >  1:1284515/1‑76 (MQ=255)                                                aaaGGCGGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGtt                     <  2:222011/76‑1 (MQ=255)                                                aaaGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATCACAGCctgctgct                     >  2:383760/1‑74 (MQ=255)                                                 aaGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc                    <  1:2059/76‑1 (MQ=255)                                                 aaGGCCGGA‑GCATGCGCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTc                    <  2:1843363/76‑1 (MQ=255)                                                  aGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCa                   >  2:1447437/1‑76 (MQ=255)                                                   ggCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCaa                  >  2:4337619/1‑76 (MQ=255)                                                    gCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCCAt                 >  2:2707458/1‑76 (MQ=255)                                                    gCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAAt                 >  2:2675674/1‑76 (MQ=255)                                                         cagcatgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGACTACAGCGCGGCGATAACAGCCTGCCGTTCAATCAGt             >  1:3225065/3‑76 (MQ=255)                                                          ggcatgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt            <  1:1604405/75‑1 (MQ=255)                                                          agcatgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt            <  1:340851/75‑1 (MQ=255)                                                          agcatgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt            >  1:4246223/2‑76 (MQ=255)                                                          agcatgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtt            <  1:3711614/75‑1 (MQ=255)                                                          agcatgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGCTCAATCAGtt            >  1:790124/2‑76 (MQ=255)                                                           gcatgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGttt           >  1:2418602/1‑76 (MQ=255)                                                            catgcTCCGGCGTT‑CGTTTTCATCACTGTGTTTTGATTACCGGGCGTCTATAGCAGCCTGCTGTTCAATCAGtttg          <  2:4299282/76‑2 (MQ=255)                                                             atgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcg         >  2:3401943/1‑76 (MQ=255)                                                             atgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcg         >  1:3613319/1‑76 (MQ=255)                                                              tgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAg              >  1:1681309/1‑70 (MQ=255)                                                              tgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAg              <  2:1681308/70‑1 (MQ=255)                                                              tgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgc        >  1:4460000/1‑76 (MQ=255)                                                              tgcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTAAATCAGtttcgc        >  1:2641174/1‑76 (MQ=255)                                                               gcTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgct       >  2:2819442/1‑76 (MQ=255)                                                                cTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctt      <  2:1193930/76‑1 (MQ=255)                                                                 tCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgct       >  2:2136424/1‑74 (MQ=255)                                                                 tCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgct       <  1:2136425/74‑1 (MQ=255)                                                                  ccGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttc    <  1:649488/76‑1 (MQ=255)                                                                  ccGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttc    >  1:1309502/1‑76 (MQ=255)                                                                   cGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttcg   >  2:1995917/1‑76 (MQ=255)                                                                    ggCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttcgc  >  2:4710/1‑76 (MQ=255)                                                                    ggCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGtttcgctttcgc  >  1:4695843/1‑76 (MQ=255)                                                                          |                                                                       TTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAGCCTGCTGTTCAATCAGTTTCGCTTTCGC  >  NC_000913/4478454‑4478597

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