Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,393	C→G	A159G (GCT→GGT)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,393	0	C	G	100.0%	199.6 / NA	56	A159G (GCT→GGT)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (18/38);  total (18/38)

GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑C‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008309‑4008458                                                                                        |                                                                           gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                                  <  2:1055363/93‑1 (MQ=255)   gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                >  2:5666223/1‑92 (MQ=255)   gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                >  1:2245207/1‑92 (MQ=255)   gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                                >  1:383212/1‑92 (MQ=255)    aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTTCC‑AAAGGCCATGCGCTGGGAGCGGTGGCG‑‑‑‑aaga                                                               >  2:5218446/1‑91 (MQ=255)    aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               <  1:4210164/93‑3 (MQ=255)    aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               >  2:3437448/1‑91 (MQ=255)    aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGAGATT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                               >  2:5286486/1‑91 (MQ=255)     tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                              <  1:2952901/93‑4 (MQ=255)     tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                              <  2:1335952/93‑4 (MQ=255)     tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                              >  1:6051407/1‑90 (MQ=255)                                      tCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAATGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  1:1901182/93‑1 (MQ=255)                                      tCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 >  2:3190494/1‑93 (MQ=255)                                      tCTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                                 <  2:2842303/93‑1 (MQ=255)                                       cTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  1:2334120/93‑2 (MQ=255)                                       cTGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                                <  1:4465851/93‑2 (MQ=255)                                        tGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:1737083/93‑3 (MQ=255)                                        tGGAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                               <  1:4989388/93‑3 (MQ=255)                                          gAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  2:226086/93‑5 (MQ=255)                                          gAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  1:4817457/93‑5 (MQ=255)                                          gAAGTGATGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                             <  2:6069960/93‑5 (MQ=255)                                                 tGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  2:4924566/93‑1 (MQ=255)                                                 tGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  1:2302993/93‑1 (MQ=255)                                                 tGGCGGGCGGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa                      <  2:2950850/93‑1 (MQ=255)                                                  ggcgggcgGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:2143856/93‑1 (MQ=255)                                                  ggcgggcgGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  2:743236/93‑1 (MQ=255)                                                  ggcgggcgGCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt                     <  1:1022573/93‑1 (MQ=255)                                                    ggggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  2:3840682/92‑1 (MQ=255)                                                    cgggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  2:3318504/93‑1 (MQ=255)                                                    cgggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa                   <  1:1615914/93‑1 (MQ=255)                                                     gggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  <  2:1588234/93‑1 (MQ=255)                                                     gggcggCGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa                  >  1:3361384/1‑93 (MQ=255)                                                      ggcggcGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  1:2267000/93‑1 (MQ=255)                                                      ggcggcGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  2:4815600/93‑1 (MQ=255)                                                      ggcggcGTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac                 <  2:4580266/93‑1 (MQ=255)                                                       gcggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGCTTGTATTGCGTTTGGTGACGGGATGAacg                >  2:4396292/1‑93 (MQ=255)                                                       gcggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                <  2:4018872/93‑1 (MQ=255)                                                       gcggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg                <  2:2466906/93‑1 (MQ=255)                                                        cggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  2:670750/93‑1 (MQ=255)                                                        cggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  1:4477690/93‑1 (MQ=255)                                                        cggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga               <  2:4254565/93‑1 (MQ=255)                                                        cggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAGGAacga               <  1:5107890/93‑1 (MQ=255)                                                         ggcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGTGTTTGGTGACGGGAGGAacgac              >  2:2158791/1‑93 (MQ=255)                                                          gcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  2:2336717/1‑93 (MQ=255)                                                          gcgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg             >  1:3018123/1‑93 (MQ=255)                                                           cgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            >  2:3055681/1‑93 (MQ=255)                                                           cgTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc            <  1:4718515/93‑1 (MQ=255)                                                            gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc           <  2:6070951/93‑1 (MQ=255)                                                                tcaaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaa        >  2:1177485/4‑93 (MQ=255)                                                                 caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       >  2:4662578/3‑93 (MQ=255)                                                                 caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:3624707/91‑1 (MQ=255)                                                                 caaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa       <  1:824842/91‑1 (MQ=255)                                                                  aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:4803706/2‑93 (MQ=255)                                                                  aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:2985712/2‑93 (MQ=255)                                                                  aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      <  1:2804031/92‑1 (MQ=255)                                                                  aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  2:4072808/2‑93 (MQ=255)                                                                  aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:5024781/2‑93 (MQ=255)                                                                  aaaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt      >  1:1681157/2‑93 (MQ=255)                                                                   aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAtc     >  2:4531756/1‑92 (MQ=255)                                                                   aaaGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg     >  2:662695/1‑93 (MQ=255)                                                                     aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:4274699/93‑1 (MQ=255)                                                                     aGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt   <  2:910103/93‑1 (MQ=255)                                                                      ggCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTg  >  2:5616454/1‑93 (MQ=255)                                                                                        |                                                                           GCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGTTT‑C‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTG  >  NC_000913/4008309‑4008458

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