Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,715,195	T→C	E400E (GAA→GAG)	bisC ←	biotin sulfoxide reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,715,195	0	T	C	100.0%	127.0 / NA	33	E400E (GAA→GAG)	bisC	biotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/17);  total (16/17)

CGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGG  >  NC_000913/3715051‑3715337                                                                                                                                                 |                                                                                                                                               cGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTcaac                                                                                                                                           >  2:129476/1‑150 (MQ=255)                       aGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTAt                                                                                                                    >  2:103218/1‑151 (MQ=255)                        ggCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATc                                                                                                                   <  1:120885/151‑1 (MQ=255)                          cGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCa                                                                                                                 <  2:264792/151‑1 (MQ=255)                              tGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCg                                                                                                              <  1:402031/150‑1 (MQ=255)                                                tAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCa                                                                                            >  1:647981/1‑150 (MQ=255)                                                   aGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCt                                                                                        <  2:103681/151‑1 (MQ=255)                                                                 gcccACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGa                                                                            >  1:18524/1‑149 (MQ=255)                                                                 gcccACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGa                                                                            >  2:605779/1‑149 (MQ=255)                                                                 gcccACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAag                                                                          <  2:521610/151‑1 (MQ=255)                                                                   ccaccaGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAagag                                                                        <  2:643736/151‑1 (MQ=255)                                                                     accaGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAagagag                                                                      <  1:571710/151‑1 (MQ=255)                                                                        aGATACAACGAATATCCGCGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCAc                                                                   <  2:477942/151‑1 (MQ=255)                                                                        aGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCGCCCGGCAAGCTGCCCTGCATGGAAGAGAGCAc                                                                   <  1:442011/151‑1 (MQ=255)                                                                          aTAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCg                                                                 <  2:94977/151‑1 (MQ=255)                                                                          aTAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCg                                                                 >  1:516647/1‑151 (MQ=255)                                                                            aaaaCGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGc                                                                >  2:18483/1‑150 (MQ=255)                                                                                 gAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAAcg                                                          >  1:252644/1‑151 (MQ=255)                                                                                     atCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAAcgccg                                                       <  1:181521/150‑1 (MQ=255)                                                                                        cGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTg                                                   >  2:407178/1‑151 (MQ=255)                                                                                            aaaTGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGt                                               <  2:237335/151‑1 (MQ=255)                                                                                               tGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTac                                            >  2:107529/1‑151 (MQ=255)                                                                                                 tCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTacc                                           >  1:189163/1‑150 (MQ=255)                                                                                                    gTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTAccacca                                       <  1:276546/151‑1 (MQ=255)                                                                                                                  gtTGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATGGt                          <  1:593131/150‑1 (MQ=255)                                                                                                                   ttGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCTCCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATgtt                          <  2:554546/149‑3 (MQ=255)                                                                                                                    tGATATGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATGGTaa                        <  2:637096/150‑1 (MQ=255)                                                                                                                       tatGCGCCACCAGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATGGTaagaa                     >  1:46142/1‑150 (MQ=255)                                                                                                                              caccaGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATGGTAAGAAAGACCaa              >  1:253587/1‑150 (MQ=255)                                                                                                                                 caGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAAcc          >  1:494069/1‑151 (MQ=255)                                                                                                                                  aGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAAccg         >  2:389335/1‑151 (MQ=255)                                                                                                                                  aGGGTTTTCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAAccg         <  2:650944/151‑1 (MQ=255)                                                                                                                                         tCCAGTGCCTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCgg  >  2:577574/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                               CGGTTTTTGCCAGGCACGGATCAGGCGATTGGTATCCTGATGATGAGTAAAGTTGGCACCGCCCGCCCACCAGATAAAACGAATATCCGGGAAATGTCGGTTCATACCGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGATTTTATCCACCGCATCGCAGCCACCCGGCAAGCTGCCCTGCATGGAAGAGAGCACCGCAGAACGCCGCGTGGGGTTACCACCATTGGCAAAATGGTAAGAAAGACCAAAACCGCCGCCGG  >  NC_000913/3715051‑3715337

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