Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,608,413	T→C	T204A (ACG→GCG)	yneE ←	bestrophin family putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,608,413	0	T	C	100.0%	199.8 / NA	54	T204A (ACG→GCG)	yneE	bestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/35);  total (19/35)

CCAGCGAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGC  >  NC_000913/1608269‑1608557                                                                                                                                                 |                                                                                                                                                 ccAGCGAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGc                                                                                                                                            <  2:44753/151‑1 (MQ=255)  cAGCGAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGc                                                                                                                                            >  1:654828/1‑150 (MQ=255)      gAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATgc                                                                                                                                       <  1:686498/151‑1 (MQ=255)            aaaGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGc                                                                                                                                 <  1:71436/151‑1 (MQ=255)                  tAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTg                                                                                                                           >  2:117384/1‑151 (MQ=255)                     gAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGcc                                                                                                                         >  1:198212/1‑150 (MQ=255)                      aaaTCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGcc                                                                                                                         >  2:808559/1‑149 (MQ=255)                          cAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGAcc                                                                                                                   <  2:49748/151‑1 (MQ=255)                          cAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGAcc                                                                                                                   <  2:378704/151‑1 (MQ=255)                          cAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGAcc                                                                                                                   <  2:270021/151‑1 (MQ=255)                          cAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGAc                                                                                                                    <  2:317684/150‑1 (MQ=255)                           aGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCg                                                                                                                  <  1:431414/151‑1 (MQ=255)                           aGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCg                                                                                                                  <  1:258943/151‑1 (MQ=255)                                  agaTAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAat                                                                                                           <  1:37183/151‑1 (MQ=255)                                      aaaaGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATc                                                                                                        >  1:407439/1‑150 (MQ=255)                                      aaaaGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATc                                                                                                        <  1:142870/150‑1 (MQ=255)                                      aaaaGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCa                                                                                                       <  2:836703/151‑1 (MQ=255)                                       aaaGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCAt                                                                                                      <  1:93118/151‑1 (MQ=255)                                        aaGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCAtt                                                                                                     <  1:319613/151‑1 (MQ=255)                                          ggCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTaa                                                                                                   <  1:720138/151‑1 (MQ=255)                                          ggCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTaa                                                                                                   <  2:451330/151‑1 (MQ=255)                                             gTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGcc                                                                                                >  1:814098/1‑151 (MQ=255)                                             gTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGcc                                                                                                <  2:684481/151‑1 (MQ=255)                                                  gTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATc                                                                                            <  2:291087/150‑1 (MQ=255)                                                       gCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAg                                                                                      <  1:217402/151‑1 (MQ=255)                                                          ggTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCt                                                                                   <  2:266478/151‑1 (MQ=255)                                                             ccACGAAAAGCGCGAACGACATCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGAAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCtaat                                                                                <  1:762532/151‑1 (MQ=255)                                                              cACGACCAGCGCGAACGGCAGCATGATCCAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGAAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCtaat                                                                                >  1:859037/1‑150 (MQ=255)                                                              cACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCtaata                                                                               <  1:808049/151‑1 (MQ=255)                                                                    cAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAgg                                                                         >  2:724000/1‑151 (MQ=255)                                                                      gcgcGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGat                                                                       <  2:350582/151‑1 (MQ=255)                                                                       cgcgAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGat                                                                       >  2:436470/1‑150 (MQ=255)                                                                          gAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCt                                                                   <  1:198095/151‑1 (MQ=255)                                                                           aaCGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTg                                                                  <  2:155215/151‑1 (MQ=255)                                                                           aaCGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTg                                                                  >  1:666255/1‑151 (MQ=255)                                                                                  gcaTGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTg                                                           <  2:62850/151‑1 (MQ=255)                                                                                      gATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTcc                                                        >  2:733974/1‑150 (MQ=255)                                                                                        tACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCAtt                                                     >  1:400681/1‑151 (MQ=255)                                                                                        tACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCAtt                                                     >  1:141084/1‑151 (MQ=255)                                                                                             aaCAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTgcggc                                                <  1:35639/151‑1 (MQ=255)                                                                                             aaCAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTgcggc                                                >  1:603233/1‑151 (MQ=255)                                                                                               cAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGgcgc                                              >  1:292534/1‑151 (MQ=255)                                                                                                 gATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTg                                            <  1:499122/151‑1 (MQ=255)                                                                                                  aTAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTg                                            >  2:645166/1‑150 (MQ=255)                                                                                                    aaaCGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAAc                                         >  1:659430/1‑151 (MQ=255)                                                                                                       cGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGc                                      <  1:203104/151‑1 (MQ=255)                                                                                                         gTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCa                                    <  2:20350/151‑1 (MQ=255)                                                                                                            cGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAAc                                  <  1:674155/150‑1 (MQ=255)                                                                                                                     aTCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCa                        <  2:664379/151‑1 (MQ=255)                                                                                                                         gggTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCattat                    <  2:832233/151‑1 (MQ=255)                                                                                                                          ggTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCattatt                   <  1:680691/151‑1 (MQ=255)                                                                                                                                       ggAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGt      >  2:236083/1‑151 (MQ=255)                                                                                                                                         aaTTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGtt     >  1:358419/1‑150 (MQ=255)                                                                                                                                            tGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGc  <  2:801446/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 CCAGCGAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGC  >  NC_000913/1608269‑1608557

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