Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 812,947 G→A *226* (TAG→TAA bioD → dethiobiotin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913812,9470GA100.0% 81.1 / NA 26*226* (TAG→TAAbioDdethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/16);  total (10/16)

GAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAT  >  NC_000913/812814‑813072
                                                                                                                                     |                                                                                                                             
gAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGTCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCAtt                                                                                                                          <  1:932541/139‑1 (MQ=255)
gAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATACCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCAtt                                                                                                                          <  1:224140/139‑1 (MQ=255)
 aaaaCGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTc                                                                                                                         <  2:724696/139‑1 (MQ=255)
       tCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTAtt                                                                                                                   <  2:904951/139‑1 (MQ=255)
                  atatGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTg                                                                                                        <  1:1065240/139‑1 (MQ=255)
                    atGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcg                                                                                                      <  2:120519/139‑1 (MQ=255)
                           cGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCgaga                                                                                               >  1:1327127/1‑139 (MQ=255)
                             cTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATc                                                                                             <  2:173588/139‑1 (MQ=255)
                             cTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATc                                                                                             <  1:1192310/139‑1 (MQ=255)
                               cACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATcgc                                                                                           <  1:677343/139‑1 (MQ=255)
                                                    gctgctGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAAttt                                                                      >  1:188712/1‑139 (MQ=255)
                                                                  tCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              >  1:227556/1‑117 (MQ=255)
                                                                  tCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              <  2:227556/117‑1 (MQ=255)
                                                                  tCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTc                                                        >  2:1209471/1‑139 (MQ=255)
                                                                   ccccTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCt                                                       >  2:334430/1‑139 (MQ=255)
                                                                              cAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCAc                                            >  2:285342/1‑139 (MQ=255)
                                                                               aGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgag                                                                                                    <  1:848744/82‑1 (MQ=255)
                                                                               aGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgag                                                                                                    >  2:848744/1‑82 (MQ=255)
                                                                                   aaTCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                                 >  1:238774/1‑97 (MQ=255)
                                                                                   aaTCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                                 <  2:238774/97‑1 (MQ=255)
                                                                                         gAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAAttt                                 <  1:334430/139‑1 (MQ=255)
                                                                                              tGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACtt                            <  1:1274829/139‑1 (MQ=255)
                                                                                                                cATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcgg               >  2:847633/1‑134 (MQ=255)
                                                                                                                cATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcgg               <  1:847633/134‑1 (MQ=255)
                                                                                                                      ccttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGa    <  2:137041/139‑1 (MQ=255)
                                                                                                                        ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt  >  2:504636/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                             
GAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAT  >  NC_000913/812814‑813072

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