Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	447,954	A→G	I31I (ATT→ATC)	cyoD ←	cytochrome o ubiquinol oxidase subunit IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	447,954	0	A	G	100.0%	205.7 / NA	55	I31I (ATT→ATC)	cyoD	cytochrome o ubiquinol oxidase subunit IV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/33);  total (22/33)

CGCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATC  >  NC_000913/447807‑448099                                                                                                                                                    |                                                                                                                                                  cGCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATc                                                                                                                                                >  2:239238/1‑151 (MQ=255)  gCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATc                                                                                                                                                >  1:205510/1‑150 (MQ=255)  gCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCa                                                                                                                                               <  2:585176/151‑1 (MQ=255)    cGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCAcc                                                                                                                                             <  1:950137/151‑1 (MQ=255)    cGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCAcc                                                                                                                                             >  1:291067/1‑151 (MQ=255)       cATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGt                                                                                                                                           <  2:729520/150‑1 (MQ=255)       cATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGt                                                                                                                                           <  1:825365/150‑1 (MQ=255)       cATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTc                                                                                                                                          <  1:157819/151‑1 (MQ=255)       cATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTc                                                                                                                                          >  1:304298/1‑151 (MQ=255)          gTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAg                                                                                                                                        <  2:22728/150‑1 (MQ=255)          gTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAg                                                                                                                                        <  1:75717/150‑1 (MQ=255)                gCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATc                                                                                                                                 <  1:532519/151‑1 (MQ=255)                    tCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGAc                                                                                                                              <  2:95236/150‑1 (MQ=255)                        cTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGAt                                                                                                                         >  2:462170/1‑151 (MQ=255)                          gATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAgg                                                                                                                           <  2:96104/147‑1 (MQ=255)                          gATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGa                                                                                                                          <  1:1063394/148‑1 (MQ=255)                          gATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGAt                                                                                                                         <  1:790237/149‑1 (MQ=255)                           aTTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATaaa                                                                                                                      <  2:679288/151‑1 (MQ=255)                               ggTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCt                                                                                                                  <  2:687778/151‑1 (MQ=255)                                gTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTg                                                                                                                 >  2:190157/1‑151 (MQ=255)                                                aGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTAc                                                                                                 >  1:764631/1‑151 (MQ=255)                                                aGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTAc                                                                                                 >  1:356138/1‑151 (MQ=255)                                                 gCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACg                                                                                                >  2:358688/1‑151 (MQ=255)                                                      ctaGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGcc                                                                                           >  2:904610/3‑151 (MQ=255)                                                        aGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCAt                                                                                         >  2:1056937/1‑151 (MQ=255)                                                          aTGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATgg                                                                                       <  2:823238/151‑1 (MQ=255)                                                          aTGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATgg                                                                                       >  1:466854/1‑151 (MQ=255)                                                            gCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATggg                                                                                      <  2:142827/150‑1 (MQ=255)                                                              aCCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGAcg                                                                                   <  1:577252/151‑1 (MQ=255)                                                                cAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGAcgcg                                                                                 <  1:950236/151‑1 (MQ=255)                                                                cAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGAcgcg                                                                                 >  1:357530/1‑151 (MQ=255)                                                                     ccTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCt                                                                            >  1:171947/1‑151 (MQ=255)                                                                       tGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCtgt                                                                          >  2:988653/1‑151 (MQ=255)                                                                              cTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTTCCATGGGACGCGCCGCTGTGATCgg                                                                    >  2:1014389/1‑150 (MQ=255)                                                                              cTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGt                                                                   >  1:945633/1‑151 (MQ=255)                                                                                     tGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATg                                                             <  2:152594/150‑1 (MQ=255)                                                                                      gCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGAc                                                           <  1:579860/151‑1 (MQ=255)                                                                                        cAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTc                                                         <  1:638848/151‑1 (MQ=255)                                                                                           gATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCAtt                                                      <  1:419436/151‑1 (MQ=255)                                                                                               gTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACAt                                                  >  1:451432/1‑151 (MQ=255)                                                                                                               ccGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATaa                                  <  1:601959/151‑1 (MQ=255)                                                                                                                cGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATaaa                                 <  2:146744/151‑1 (MQ=255)                                                                                                                cGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATaa                                  >  1:175370/1‑150 (MQ=255)                                                                                                                 ggAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAaac                                <  1:832533/151‑1 (MQ=255)                                                                                                                    gagGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAaacaac                             <  2:390495/151‑1 (MQ=255)                                                                                                                       gCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGt                          >  1:704180/1‑151 (MQ=255)                                                                                                                         aGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGa                        <  1:198831/151‑1 (MQ=255)                                                                                                                             ccTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAacac                    <  2:566579/151‑1 (MQ=255)                                                                                                                             ccTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAac                      <  2:230570/149‑1 (MQ=255)                                                                                                                               tGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAacacac                  >  1:431861/1‑151 (MQ=255)                                                                                                                                   atcaCCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGAt              <  2:131477/151‑1 (MQ=255)                                                                                                                                        cATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCaaa         <  2:277061/151‑1 (MQ=255)                                                                                                                                               aaCGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATc  <  2:534290/151‑1 (MQ=255)                                                                                                                                               aaCGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATc  <  1:431341/151‑1 (MQ=255)                                                                                                                                                 cGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATc  >  2:493412/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CGCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATC  >  NC_000913/447807‑448099

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