Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,810,304	G→A	*160* (TAG→TAA)	coaD →	pantetheine‑phosphate adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,810,304	0	G	A	100.0%	124.5 / NA	39	*160* (TAG→TAA)	coaD	pantetheine‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/19);  total (20/19)

GGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACA  >  NC_000913/3810195‑3810423                                                                                                              |                                                                                                                        ggTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGgc                                                                                            >  2:1252376/1‑139 (MQ=255)          tctcTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGaa                                                                                   >  1:378904/1‑139 (MQ=255)            tctTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAtt                                                                                 <  2:378904/139‑1 (MQ=255)            tctTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAtt                                                                                 >  1:523180/1‑139 (MQ=255)                cATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGcc                                                                                                               >  2:461408/1‑105 (MQ=255)                cATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGcc                                                                                                               <  1:1321918/105‑1 (MQ=255)                cATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGcc                                                                                                               <  1:461408/105‑1 (MQ=255)                cATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGcc                                                                                                               >  2:1321918/1‑105 (MQ=255)                cATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACtt                                                                             >  1:364256/1‑139 (MQ=255)                      tGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCa                                                                       <  2:600691/139‑1 (MQ=255)                      tGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCa                                                                       >  1:375223/1‑139 (MQ=255)                        gTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACt                                                                     >  2:18449/1‑139 (MQ=255)                           aaaGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGcc                                                                  <  2:327635/139‑1 (MQ=255)                             agagGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGa                                                                <  1:1252376/139‑1 (MQ=255)                                    gcgcgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGgcgcg                                                                                         <  1:1208841/107‑1 (MQ=255)                                    gcgcgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGgcgcg                                                                                         >  2:1208841/1‑107 (MQ=255)                                         ccATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGc                                                                   >  1:1308593/1‑124 (MQ=255)                                         ccATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGc                                                                   <  2:1308593/124‑1 (MQ=255)                                         ccATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGtt                                                    <  2:364256/139‑1 (MQ=255)                                                    aTGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTgcgc                                         >  2:800674/1‑139 (MQ=255)                                                        cACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGt                                     <  2:375223/139‑1 (MQ=255)                                                             aTTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGt                                >  2:994505/1‑139 (MQ=255)                                                               ttCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTg                                            >  2:156901/1‑125 (MQ=255)                                                               ttCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTg                                            <  1:156901/125‑1 (MQ=255)                                                                   tGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGcc                                                                  <  2:12808/99‑1 (MQ=255)                                                                   tGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGcc                                                                  >  1:12808/1‑99 (MQ=255)                                                                             tGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATg                                                                                                          >  1:622633/1‑49 (MQ=255)                                                                             tGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATg                                                                                                          <  2:622633/49‑1 (MQ=255)                                                                               tCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCgg                                                                                                             >  2:597156/1‑44 (MQ=255)                                                                               tCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCgg                                                                                                             <  1:597156/44‑1 (MQ=255)                                                                                ccATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCg             <  2:523180/139‑1 (MQ=255)                                                                                    cAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTTCCGcaca         <  1:994505/139‑1 (MQ=255)                                                                                     aGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCgg                                                                                                             >  1:476385/1‑38 (MQ=255)                                                                                     aGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCgg                                                                                                             <  2:476385/38‑1 (MQ=255)                                                                                        cgcTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCg     >  1:1226971/1‑139 (MQ=255)                                                                                         gcTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAt                                                                                  >  2:598247/1‑61 (MQ=255)                                                                                         gcTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAt                                                                                  <  1:598247/61‑1 (MQ=255)                                                                                         gcTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGa    >  2:910942/1‑139 (MQ=255)                                                                                           tGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGaca  <  1:910942/139‑1 (MQ=255)                                                                                                        agcgtaacgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCcac                          <  2:105679/102‑1 (MQ=255)                                                                                                        agcgtaacgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCcac                          >  1:105679/1‑102 (MQ=255)                                                                                                          cgtaacgtTTATGCCGGATGGTATGCCATCCGGCGCGCa                                                                                       >  2:780570/1‑39 (MQ=255)                                                                                                          cgtaacgtTTATGCCGGATGGTATGCCATCCGGCGCGCa                                                                                       <  1:780570/39‑1 (MQ=255)                                                                                                              |                                                                                                                        GGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACA  >  NC_000913/3810195‑3810423

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