Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,478,503	(A)7→6	intergenic (+12/+47)	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,497	0	A	.	100.0%	330.0 / NA	72	intergenic (+6/+53)	yjgN/valS	hypothetical protein/valine tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (39/33);  total (39/33)

AGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAG  >  NC_000913/4478422‑4478569                                                                            |                                                                         agaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATa                                                                          <  1:4133326/76‑1 (MQ=255) agaaCAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATa                                                                          <  1:4114567/76‑1 (MQ=255)   aacaacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaa                                                                        >  2:84907/1‑76 (MQ=255)   aacaacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaa                                                                        >  1:444812/1‑76 (MQ=255)     caacCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaa                                                                      <  1:2441216/76‑1 (MQ=255)       acCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaag                                                                    >  2:4651324/1‑75 (MQ=255)       acCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaag                                                                    >  1:3310521/1‑75 (MQ=255)       acCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaag                                                                    >  1:2331284/1‑75 (MQ=255)         cAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATaaaaaaggc                                                                  <  1:1222308/76‑4 (MQ=255)            aTAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCgg                                                              >  1:3070893/1‑76 (MQ=255)             tAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGa                                                             <  1:1093040/76‑1 (MQ=255)             tAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGa                                                             >  1:2251271/1‑76 (MQ=255)              aaaGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGag                                                            >  1:2722329/1‑75 (MQ=255)               aaGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGagc                                                           <  2:1207543/76‑3 (MQ=255)                aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGagca                                                          <  2:1221974/76‑4 (MQ=255)                aGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGagca                                                          >  2:2529601/1‑73 (MQ=255)                  gCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGga‑gcatg                                                       >  2:1501711/1‑76 (MQ=255)                   cttcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGga‑gcatgc                                                      >  1:1686183/1‑76 (MQ=255)                    gtcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCt                                                     <  2:1490776/75‑1 (MQ=255)                     tcttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTc                                                    <  1:3603145/76‑1 (MQ=255)                      cttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTcc                                                   <  1:4085748/76‑1 (MQ=255)                      cttAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTcc                                                   <  2:444812/76‑1 (MQ=255)                       ttAGGCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCg                                                  <  2:4094229/76‑1 (MQ=255)                          gCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGcc                                               >  1:125426/1‑76 (MQ=255)                           ccAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCt                                              <  1:3161219/76‑1 (MQ=255)                            cAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCtt                                             >  1:1402123/1‑76 (MQ=255)                             aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCttc                                            <  1:4726168/76‑2 (MQ=255)                             aGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCttc                                            >  1:4238140/1‑75 (MQ=255)                              gTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCttcg                                           >  1:367421/1‑74 (MQ=255)                               tATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCttcgt                                          <  2:4238140/76‑4 (MQ=255)                                 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                              >  2:4741804/1‑76 (MQ=255)                                          gAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACtg                              <  1:4321047/76‑1 (MQ=255)                                          gAGTTATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACtg                              <  1:1256798/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:2912015/76‑1 (MQ=255)                                              tATGCCTTCTTTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGtt                          >  2:422308/1‑76 (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:2234281/1‑76 (MQ=255)                                                   cttcttTACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGAt                     <  2:4163031/76‑1 (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:367421/76‑1 (MQ=255)                                                         tACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAgc               <  1:733602/76‑1 (MQ=255)                                                         tACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAgc               >  2:2246381/1‑76 (MQ=255)                                                         tACCATTTCTGTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAgc               >  2:4703852/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     <  2:4432840/76‑1 (MQ=255)                                                                   gTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa     >  1:1225971/1‑76 (MQ=255)                                                                   gTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa     <  2:603731/76‑1 (MQ=255)                                                                   gTAAGAGAT‑AAAAAAGGCCGGA‑GCATGCTCCGGCCTT‑CGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATaa     >  1:831207/1‑76 (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)                                                                            |                                                                         AGAACAACCAGATAAAGGCTTCTTAGCCAGTATTTCACGGGGAGTTATGCCTTCTTTACCATTTCTGTAAGAGATAAAAAAAGGCCGGATGCATGCTCCGGCCTTGCGTTTTCATCACTGTGTTTTGATTACAGCGCGGCGATAACAG  >  NC_000913/4478422‑4478569

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