Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,609,015	A→G	V3A (GTT→GCT)	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,609,015	0	A	G	100.0%	148.6 / NA	41	V3A (GTT→GCT)	yneE	bestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/15);  total (26/15)

GAATTTGATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCAT  >  NC_000913/1608867‑1609158                                                                                                                                                     |                                                                                                                                                gAATTTGATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCa                                                                                                                                               <  2:420873/151‑1 (MQ=255)                 tGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTc                                                                                                                                >  1:636445/1‑150 (MQ=255)                  gCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTc                                                                                                                                >  2:372264/1‑149 (MQ=255)                     tgtACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCg                                                                                                                           >  2:222726/1‑151 (MQ=255)                        aCCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTgg                                                                                                                        >  2:351857/1‑151 (MQ=255)                          cAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGaa                                                                                                                      <  2:75700/151‑1 (MQ=255)                            gggCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGaaaa                                                                                                                    <  2:273983/151‑1 (MQ=255)                                agcaTGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCg                                                                                                                >  2:37550/1‑151 (MQ=255)                                   aTGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGccc                                                                                                             >  2:553747/1‑151 (MQ=255)                                       aaaTAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTAtt                                                                                                         >  1:492109/1‑151 (MQ=255)                                             cAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGt                                                                                                   <  1:253287/151‑1 (MQ=255)                                                cGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGAt                                                                                                >  1:654149/1‑151 (MQ=255)                                                cGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGAt                                                                                                >  2:416221/1‑151 (MQ=255)                                                   tAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTAc                                                                                             >  2:310959/1‑151 (MQ=255)                                                    aGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACg                                                                                            <  1:568538/151‑1 (MQ=255)                                                          gagaAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATcac                                                                                       >  2:665434/1‑150 (MQ=255)                                                                    agtagtAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGAc                                                                            <  2:558787/151‑1 (MQ=255)                                                                     gtagtaAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACg                                                                           <  2:112671/151‑1 (MQ=255)                                                                           aGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACtt                                                                     >  1:474153/1‑151 (MQ=255)                                                                           aGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACtt                                                                     >  1:350039/1‑151 (MQ=255)                                                                                  ggATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAAtt                                                              <  1:35417/151‑1 (MQ=255)                                                                                   gATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTa                                                             <  2:654160/151‑1 (MQ=255)                                                                                      atCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCt                                                          >  1:221516/1‑151 (MQ=255)                                                                                          tGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACa                                                       <  1:15161/150‑1 (MQ=255)                                                                                             ataataCTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCt                                                   <  2:35064/151‑1 (MQ=255)                                                                                             ataataCTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCt                                                   >  2:506564/1‑151 (MQ=255)                                                                                                   cTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTacaac                                             >  1:6789/1‑151 (MQ=255)                                                                                                      aGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACAtt                                          >  1:295162/1‑151 (MQ=255)                                                                                                      aGCCGTGCCAGACAAAAATACGGCGCAGACAGTGTTGCTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGAAAAAGTCACCCTATTAAACGTGATTACGATCACATTCTCAAGACCCTAATTACAAATTACCTACAAGCTTACAACAtt                                          >  1:519032/1‑151 (MQ=255)                                                                                                            gCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAgg                                    >  2:28983/1‑151 (MQ=255)                                                                                                                 aCAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTg                               >  2:179450/1‑151 (MQ=255)                                                                                                                   aaaaaTACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGc                              >  1:429230/1‑150 (MQ=255)                                                                                                                      aaTACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCAtgt                          <  2:492111/151‑1 (MQ=255)                                                                                                                      aaTACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCAtgt                          >  2:200575/1‑151 (MQ=255)                                                                                                                      aaTACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCAtgt                          <  2:577728/151‑1 (MQ=255)                                                                                                                        tACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCtt                                                  <  1:251210/125‑1 (MQ=255)                                                                                                                        tACGGCGCAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCtt                                                  >  2:251211/1‑125 (MQ=255)                                                                                                                              gcAGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTg                  <  2:295164/151‑1 (MQ=255)                                                                                                                                atccAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGtt                >  2:120314/3‑151 (MQ=255)                                                                                                                                aGCCAGTGTTGTTGTGGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGtt                >  2:183174/1‑151 (MQ=255)                                                                                                                                              tgGACGAGCAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCat  >  1:630472/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                GAATTTGATGCCCAGATGCGTGTACCAGGGCAGCATGAAAATAACAGCGATAGAAAAGAGAAAATTGAGTAGTAAGCGCGAGGATATCTTGGATAATACTGAGCCGTGCCAGACAAAAATACGGCGCAGCCAGTGTTGTTGTGGACGAACAATCATGGTTATCTTCAGGCGTGGAAAAGTCGCCCTATTAAACGTGATTACGATCACATTCTCAAGACGCTACTTACAAATTACCTACAAGCTTACAACATTAGCAGGCGCTGCATGTGGCACCTGTTAATGATAAAGGCAT  >  NC_000913/1608867‑1609158

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