Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,715,561	G→T	F278L (TTC→TTA)	bisC ←	biotin sulfoxide reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,715,561	0	G	T	100.0%	81.8 / NA	25	F278L (TTC→TTA)	bisC	biotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/12);  total (13/12)

CTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCATTTT  >  NC_000913/3715425‑3715692                                                                                                                                         |                                                                                                                                    cTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTTGCCCAGCAAATAAGAGGCTAAg                                                                                                                                  <  1:303799/140‑1 (MQ=255) cTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAg                                                                                                                                  <  1:44105/140‑1 (MQ=255)          gCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATaa                                                                                                                         >  2:703703/1‑140 (MQ=255)                   tATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTTGCTATTTCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTgtggtg                                                                                                                >  1:486506/1‑140 (MQ=255)                             gAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGGTTTCGCTATTCCGTCACTCTCGCCCAACAAATAAGAAGCCAAGACGGCACAAACTGTAGTGCAACGccccc                                                                                                      >  2:547081/1‑139 (MQ=255)                                     ccgccAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTATGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCATCGCGCCAGAAACGCt                                                                                              <  2:492596/140‑1 (MQ=255)                                         cAGCACGCGGTTTTTGGCTGCGCTAACACCACAAAGCTCTGCTGCTCATTCGGCAGTTTTCGCTATTCCGTCTCTCTCGCCCAGCAAATAAGTGGCTAAGACGGAATAACTTGTGGTGCAACGCGCCAGAAACGCTtcgt                                                                                          <  2:600654/140‑1 (MQ=255)                                                cgGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCt                                                                                              >  1:409801/1‑129 (MQ=255)                                                cgGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCt                                                                                              <  2:409801/129‑1 (MQ=255)                                                  gATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAAc                                                                                 <  2:384998/140‑1 (MQ=255)                                                                                tgctgcCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATcccc                                                    >  1:50181/1‑139 (MQ=255)                                                                                tgctgcCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATcccc                                                    <  2:50181/139‑1 (MQ=255)                                                                                tgctgcCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCa                                                   >  2:439483/1‑140 (MQ=255)                                                                                     cccATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCAt                                               <  1:517154/139‑1 (MQ=255)                                                                                     cccATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCAt                                               >  2:517154/1‑139 (MQ=255)                                                                                                             cGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCAtgtg                       >  2:304810/1‑139 (MQ=255)                                                                                                             cGTCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCAtgtg                       <  1:304810/139‑1 (MQ=255)                                                                                                               tCACTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCgg                    >  1:170047/1‑140 (MQ=255)                                                                                                                 tctcaCGCCCAGCAAAGAAGGGGCTAAGGTGGCATAACCTATGGTGCAACGCGCCAGAAGCGGTGCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACAGCGGTGCCCATGTGCGGTg                  <  2:133581/139‑1 (MQ=255)                                                                                                                 aCTCTCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTg                  <  2:170047/140‑1 (MQ=255)                                                                                                                  ctctCGCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCACAAACGCTTCGTCGTGCCAACCATTTTCCACCCGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGc                 >  2:24660/1‑140 (MQ=255)                                                                                                                       gCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGcc                >  2:624480/1‑136 (MQ=255)                                                                                                                       gCCCAGCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGcc                <  1:624480/136‑1 (MQ=255)                                                                                                                            gCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCAtgtg                       >  1:561198/1‑124 (MQ=255)                                                                                                                            gCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCAtgtg                       <  2:561198/124‑1 (MQ=255)                                                                                                                            gCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGccacccac          <  1:436745/137‑1 (MQ=255)                                                                                                                            gCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGccacccac          >  2:436745/1‑137 (MQ=255)                                                                                                                            gCAAATAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTcc       >  2:488028/1‑140 (MQ=255)                                                                                                                                 tAAGAGGCTAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCAtttt  >  2:394174/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                    CTGCCATCAGCATGGTGGTATTTTGGTGGAAAATAGCCGCCAGCTCGCGGATTTTCGCTGCGCCAACACCACAAATCTCTGCTGCCCATTCGGCGGTTTTCGCTATTCCGTCACTCTCGCCCAGCAAATAAGAGGCGAAGACGGCATAACCTGTGGTGCAACGCGCCAGAAACGCTTCGTCGTGCCAACCATTTTCCACCAGCGTATGGGCGATCCCCAGCATCAGCGCAACATCGGTGCCCATGTGCGGTGCCACCCACTCCATTTT  >  NC_000913/3715425‑3715692

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