Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,375	+A	coding (458/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,372	1	.	A	100.0%	174.7 / NA	51	S152* (TCA→TAA)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (18/33);  total (18/33)

GCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑C‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCC  >  NC_000913/4008280‑4008449                                                                                                |                                                                                        gCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑c                                                                                          <  2:7242/93‑1 (MQ=255)   tGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑C‑aa                                                                                       <  2:2490160/93‑1 (MQ=255)   tGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑C‑aa                                                                                       >  1:525662/1‑93 (MQ=255)    ggAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑C‑aaa                                                                                      >  2:3084901/1‑93 (MQ=255)    ggAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑C‑aaa                                                                                      >  1:2066963/1‑93 (MQ=255)    ggAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑C‑aaa                                                                                      <  1:1945596/93‑1 (MQ=255)    ggAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GCTT‑C‑caa                                                                                      >  2:776144/1‑93 (MQ=255)     gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑C‑aaaa                                                                                     <  1:3977048/93‑2 (MQ=255)     gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑C‑aaaa                                                                                     <  1:2502323/93‑2 (MQ=255)     gAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑C‑aaaa                                                                                     <  2:1325041/93‑2 (MQ=255)          ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTTCC‑AAAGGCCa                                                                                 >  2:3850791/1‑93 (MQ=255)          ggCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCa                                                                                 <  2:1748311/93‑1 (MQ=255)           gCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCAt                                                                                <  2:2370098/93‑1 (MQ=255)           gCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCAt                                                                                <  1:4177035/93‑1 (MQ=255)             gATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGCGGTTT‑A‑AAAAGCCATgg                                                                              >  1:4079721/1‑92 (MQ=255)              aTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCCTgcc                                                                             >  2:3837211/1‑92 (MQ=255)              aTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATgcg                                                                             >  1:2791068/1‑93 (MQ=255)              aTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCAATgcg                                                                             <  1:3578933/93‑1 (MQ=255)              aTTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAAGCCATgcg                                                                             >  2:4277203/1‑93 (MQ=255)               ttAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATgcgc                                                                            <  1:4375908/93‑1 (MQ=255)                tAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCt                                                                           <  2:1283216/93‑1 (MQ=255)                tAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCt                                                                           <  1:3084383/93‑1 (MQ=255)                 aaCGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAAGACATGCGCTg                                                                          >  2:1584583/1‑93 (MQ=255)                   cGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGa                                                                        <  1:1207242/93‑1 (MQ=255)                   agCTCGTTGGGGCGATCGCGTCAACGACAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCAGC‑GTTT‑CAAAAGGCCATGCGCTGGa                                                                        <  2:1465182/92‑1 (MQ=255)                      tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAgc                                                                     <  1:3133448/93‑2 (MQ=255)                      tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAgc                                                                     <  1:1906222/93‑2 (MQ=255)                      tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAg‑g                                                                    >  2:3840507/1‑93 (MQ=255)                      tCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GCTT‑CAAAAGGCCATGCGCTGGGAg‑g                                                                    >  2:1041082/1‑93 (MQ=255)                        gTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAAGCCATGCGCTGGAAgcgg                                                                   >  1:1516109/1‑90 (MQ=255)                         ttGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAgcggt                                                                  <  2:4049354/93‑5 (MQ=255)                              gCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGG‑GCTT‑CAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑a                                                         >  2:3326159/1‑93 (MQ=255)                               cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                        <  2:3966728/93‑1 (MQ=255)                               cGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aa                                                        <  1:1782384/93‑1 (MQ=255)                                gATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aag                                                       <  1:1177000/93‑2 (MQ=255)                                 aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                      >  1:46973/1‑91 (MQ=255)                                 aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                      >  1:3511183/1‑91 (MQ=255)                                 aTCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAAGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aaga                                                      >  1:3215869/1‑91 (MQ=255)                                  tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                     <  1:3329958/93‑4 (MQ=255)                                  tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                     >  2:2281928/1‑90 (MQ=255)                                  tCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCG‑‑‑‑aagaa                                                     >  2:3445805/1‑90 (MQ=255)                                                                   tCTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtt                        <  2:1756261/93‑1 (MQ=255)                                                                    gtGGAAGTGATGGCGGGCGGC‑GTGT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                       <  2:1612774/92‑2 (MQ=255)                                                                    cTGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGttt                       <  2:2585349/93‑2 (MQ=255)                                                                     tGGAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttg                      >  1:3760002/1‑91 (MQ=255)                                                                      ggAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                     >  1:1504553/1‑90 (MQ=255)                                                                      ggAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttgg                     >  1:4328546/1‑90 (MQ=255)                                                                       gAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                    <  2:3225814/93‑5 (MQ=255)                                                                       gAAGTGATGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGtttggt                    >  2:362797/1‑89 (MQ=255)                                                                             gtGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg              <  1:3290333/92‑1 (MQ=255)                                                                             aTGGCGGGCGGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACggg              <  1:980410/93‑1 (MQ=255)                                                                              tGGCGGGCGGC‑GTTT‑CAAAAGGCAATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGa             <  2:2605564/93‑1 (MQ=255)                                                                               ggcgggcgGC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGAt            <  2:1250992/93‑1 (MQ=255)                                                                                gcgggcggC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATg           >  2:4734727/1‑93 (MQ=255)                                                                                 cgggcggC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa          <  2:1009336/93‑1 (MQ=255)                                                                                 cgggcggC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGa          <  2:2641560/93‑1 (MQ=255)                                                                                  gggcggC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa         <  2:594135/93‑1 (MQ=255)                                                                                  gggcggC‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGaa         <  1:1004691/93‑1 (MQ=255)                                                                                   ggcggc‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac        >  2:2025077/1‑93 (MQ=255)                                                                                   ggcggc‑GTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTAGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAac        >  2:1242468/1‑93 (MQ=255)                                                                                    gcggc‑gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacg       <  1:3802206/93‑1 (MQ=255)                                                                                     cggc‑gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacga      <  1:2414018/93‑1 (MQ=255)                                                                                      ggc‑gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac     <  2:907022/93‑1 (MQ=255)                                                                                      ggc‑gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgac     <  1:321200/93‑1 (MQ=255)                                                                                       gc‑gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg    >  1:978131/1‑93 (MQ=255)                                                                                       gc‑gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAacgacg    >  1:1848270/1‑93 (MQ=255)                                                                                        c‑gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGc   <  2:1666186/93‑1 (MQ=255)                                                                                          gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc  <  2:4673438/93‑1 (MQ=255)                                                                                          gTTT‑CAAAAGGCCATGCGCTGGAAGCGGTGGCGAAGAAACTGGGCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGcc  >  1:1095282/1‑93 (MQ=255)                                                                                                |                                                                                        GCTGGAGCAGGCGATTAACGCTCGTTGGGGCGATCGCGTCAACGTCAGTTTCTCTACCTTAACCTGTCTGGAAGTGATGGCGGGCGGC‑GTTT‑C‑AAAGGCCATGCGCTGGAAG‑GCTGGCG‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATTG‑ATTGCGTT‑‑GTGACGGGATGAACGACGCC  >  NC_000913/4008280‑4008449

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