Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,391	+C	coding (474/516 nt)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,391	1	.	C	100.0%	147.0 / NA	47	K158N (AAG→AAC)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (22/25);  total (22/25)

TCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCG‑TT‑‑G‑TGACGGGATGAACGACGCCGAAATGCT  >  NC_000913/4008301‑4008457                                                                                             |                                                                              tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                               <  1:1906222/93‑2 (MQ=255) tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgc                                                                               <  1:3133448/93‑2 (MQ=255) tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAg‑g                                                                              >  2:3840507/1‑93 (MQ=255) tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGCTTCAAAAGGCCATGCGCTGGGAg‑g                                                                              >  2:1041082/1‑93 (MQ=255)   gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAgcgg                                                                             >  1:1516109/1‑90 (MQ=255)    ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAgcggt                                                                            <  2:4049354/93‑5 (MQ=255)         gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGGGCTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                   >  2:3326159/1‑93 (MQ=255)          cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                  <  1:1782384/93‑1 (MQ=255)          cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                                  <  2:3966728/93‑1 (MQ=255)           gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                 <  1:1177000/93‑2 (MQ=255)            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                >  1:46973/1‑91 (MQ=255)            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                >  1:3511183/1‑91 (MQ=255)            aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                                >  1:3215869/1‑91 (MQ=255)             tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                               >  2:2281928/1‑90 (MQ=255)             tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                               <  1:3329958/93‑4 (MQ=255)             tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                               >  2:3445805/1‑90 (MQ=255)                                              tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tt                                 <  2:1756261/93‑1 (MQ=255)                                               gtGGAAGTGATGGCGGGCGGCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑ttt                                <  2:1612774/92‑2 (MQ=255)                                               cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑ttt                                <  2:2585349/93‑2 (MQ=255)                                                tGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tttg                               >  1:3760002/1‑91 (MQ=255)                                                 ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tttgg                              >  1:4328546/1‑90 (MQ=255)                                                 ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tttgg                              >  1:1504553/1‑90 (MQ=255)                                                  gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tttggt                             >  2:362797/1‑89 (MQ=255)                                                  gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑tttggt                             <  2:3225814/93‑5 (MQ=255)                                                        gtGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACggg                      <  1:3290333/92‑1 (MQ=255)                                                        aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACggg                      <  1:980410/93‑1 (MQ=255)                                                         tGGCGGGCGGCGTTTCAAAAGGCAATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGa                     <  2:2605564/93‑1 (MQ=255)                                                          ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGAt                    <  2:1250992/93‑1 (MQ=255)                                                           gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATg                   >  2:4734727/1‑93 (MQ=255)                                                            cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGa                  <  2:1009336/93‑1 (MQ=255)                                                            cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGa                  <  2:2641560/93‑1 (MQ=255)                                                             gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGaa                 <  1:1004691/93‑1 (MQ=255)                                                             gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGaa                 <  2:594135/93‑1 (MQ=255)                                                              ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAac                >  2:2025077/1‑93 (MQ=255)                                                              ggcggcGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTAGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAac                >  2:1242468/1‑93 (MQ=255)                                                               gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacg               <  1:3802206/93‑1 (MQ=255)                                                                cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacga              <  1:2414018/93‑1 (MQ=255)                                                                 ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacgac             <  1:321200/93‑1 (MQ=255)                                                                 ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacgac             <  2:907022/93‑1 (MQ=255)                                                                  gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacgacg            >  1:1848270/1‑93 (MQ=255)                                                                  gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAacgacg            >  1:978131/1‑93 (MQ=255)                                                                   cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGc           <  2:1666186/93‑1 (MQ=255)                                                                    gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGcc          <  2:4673438/93‑1 (MQ=255)                                                                    gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGcc          >  1:1095282/1‑93 (MQ=255)                                                                       tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGaa       >  2:3850580/4‑93 (MQ=255)                                                                       gcaaaaGGCCCTGCGCTGTAAGCGCTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGaa       <  2:3861309/90‑1 (MQ=255)                                                                        caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGaaa      >  1:2116878/3‑93 (MQ=255)                                                                        caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGaaa      <  2:643823/91‑1 (MQ=255)                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTT‑‑GTTGACGGGAGGAACTACGCCGAAAt     >  2:3805337/2‑93 (MQ=255)                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAAt     >  2:2436174/2‑93 (MQ=255)                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAAt     >  2:208377/2‑93 (MQ=255)                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAAt     <  2:2079373/92‑1 (MQ=255)                                                                         aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAAt     <  1:3137890/92‑1 (MQ=255)                                                                          aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAATg    >  1:3644293/1‑93 (MQ=255)                                                                          aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAATg    >  2:991364/1‑93 (MQ=255)                                                                            aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCG‑TTTGG‑TGACGGGATGAACGACGCCGAAATGCt  >  1:1869543/1‑93 (MQ=255)                                                                                             |                                                                              TCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCG‑TT‑‑G‑TGACGGGATGAACGACGCCGAAATGCT  >  NC_000913/4008301‑4008457

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