Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,714,680	T→C	N572S (AAC→AGC)	bisC ←	biotin sulfoxide reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,714,680	0	T	C	100.0%	163.9 / NA	47	N572S (AAC→AGC)	bisC	biotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/26);  total (21/26)

CCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGTTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATAAAACGTTTCCAGCCATTGCAGCTCACTTTTTCC  >  NC_000913/3714545‑3714810                                                                                                                                        |                                                                                                                                   ccAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTc                                                                                                                                >  2:435907/1‑140 (MQ=255) ccAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTc                                                                                                                                >  1:637516/1‑140 (MQ=255)   aGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCg                                                                                                                              <  1:435907/140‑1 (MQ=255)    gggCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCgg                                                                                                                             <  1:887493/140‑1 (MQ=255)    gggCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCgg                                                                                                                             <  1:752687/140‑1 (MQ=255)      gCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCa                                                                                                                           >  1:21812/1‑140 (MQ=255)         aTCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCAtct                                                                                                                        <  2:21812/140‑1 (MQ=255)          tCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCAtctc                                                                                                                       >  2:909368/1‑140 (MQ=255)          tCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCAtctc                                                                                                                       >  1:25699/1‑140 (MQ=255)            cGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCaa                                                                                                                     <  1:909368/140‑1 (MQ=255)               gTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTa                                                                                                                  >  1:733590/1‑140 (MQ=255)                   ccGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTg                                                                                                              <  2:637516/140‑1 (MQ=255)                   ccGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTg                                                                                                              <  1:60107/140‑1 (MQ=255)                             cAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGc                                                                                                    <  1:897380/140‑1 (MQ=255)                             cAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGc                                                                                                    >  2:249551/1‑140 (MQ=255)                             cAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGc                                                                                                    >  1:645723/1‑140 (MQ=255)                               aTACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCa                                                                                                  <  2:25699/140‑1 (MQ=255)                                 aCGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGa                                                                                                <  1:598194/140‑1 (MQ=255)                                               tttCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATggc                                                                                  >  1:668739/1‑140 (MQ=255)                                                ttCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATggcg                                                                                 <  2:668739/140‑1 (MQ=255)                                                 tCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATggcgg                                                                                <  2:475063/140‑1 (MQ=255)                                                           cGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAAcc                                                                      <  2:421522/140‑1 (MQ=255)                                                             cTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACctg                                                                    <  1:738218/140‑1 (MQ=255)                                                             cTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACctg                                                                    <  1:942421/140‑1 (MQ=255)                                                                            ggATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCcgcg                                                     <  1:783700/140‑1 (MQ=255)                                                                             gATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGTAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCgcgc                                                    <  1:345166/140‑1 (MQ=255)                                                                              aTGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCcgcg                                                     >  1:117156/1‑138 (MQ=255)                                                                              aTGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCcgcg                                                     <  2:117156/138‑1 (MQ=255)                                                                              aTGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCcgcg                                                     <  2:260751/138‑1 (MQ=255)                                                                              aTGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCcgcg                                                     >  1:260751/1‑138 (MQ=255)                                                                              aTGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCt                                                   >  1:481422/1‑140 (MQ=255)                                                                              aTGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCt                                                   >  2:308320/1‑140 (MQ=255)                                                                                    cAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCggg                                             <  2:733590/140‑1 (MQ=255)                                                                                    cAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCggg                                             <  2:715979/140‑1 (MQ=255)                                                                                      gCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCa                                           >  2:76630/1‑140 (MQ=255)                                                                                          aTCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGt                                       >  1:659707/1‑140 (MQ=255)                                                                                            cgcgGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTta                                     <  2:481422/140‑1 (MQ=255)                                                                                                 cAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATaaaa                                >  2:526154/1‑140 (MQ=255)                                                                                                           cgaagcgaaTAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATaaaa                                <  2:71799/130‑1 (MQ=255)                                                                                                           cgaagcgaaTAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATaaaa                                >  1:71799/1‑130 (MQ=255)                                                                                                                cgaaTAAACCGCTCGCTGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATAAAACGTTTCCAGCCATTg                 <  2:659707/140‑1 (MQ=255)                                                                                                                           ctcgctGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATAAAACGTTTCCAGCCATTGCAGCTCACttt      >  1:469891/1‑140 (MQ=255)                                                                                                                           ctcgctGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATAAAACGTTTCCAGCCATTGCAGCTCACttt      >  1:831882/1‑140 (MQ=255)                                                                                                                            tcgctGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATAAAACGTTTCCAGCCATTGCAGCTCACtttt     >  1:206333/1‑140 (MQ=255)                                                                                                                              gctGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATAAAACGTTTCCAGCCATTGCAGCTCACTTTttc   <  2:206333/140‑1 (MQ=255)                                                                                                                              gctGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATAAAACGTTTCCAGCCATTGCAGCTCACTTTttc   >  1:186432/1‑140 (MQ=255)                                                                                                                               ctGTCCGGGCTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATAAAACGTTTCCAGCCATTGCAGCTCACTTTttcc  <  2:469891/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   CCAGGGCAATCCGGGTAACCGTAATCGGCAATACGCTGTGAGAAGATTTCAATCTTGCCGCTGGCGGTTTTTAACGGATGCGCCAGCGGATCGCGGCAAAAATCAGCGAAGCGAATAAACCGCTCGCTGTCCGGGTTTTCCGGCATCTCAATTAACTGGTTGGCTTGCCAGAACTCAGCAAATGGCGGCAATTCAACCTGCTGGCTTGCCCCGCGCTGTCGGGCAACGTTATAAAACGTTTCCAGCCATTGCAGCTCACTTTTTCC  >  NC_000913/3714545‑3714810

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