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%	250.0 / NA	67	T204A (ACG→GCG)	yneE	bestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/43);  total (24/43)

TCCAGCGAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTG  >  NC_000913/1608268‑1608559                                                                                                                                                  |                                                                                                                                                   tCCAGCGAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAACTGGCGCATAgg                                                                                                                                               >  1:633297/1‑151 (MQ=255)      cGAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATg                                                                                                                                          >  2:564204/1‑151 (MQ=255)        aaaTAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATgcg                                                                                                                                        <  2:3144/151‑1 (MQ=255)        aaaTAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATgcg                                                                                                                                        <  2:272803/151‑1 (MQ=255)            aaaaGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCg                                                                                                                                    >  1:64349/1‑151 (MQ=255)              aaGTGTAGGAAATCAGCGCAGAGATACAATGCGTCATTTAATTCAGGTCCACGACCATAGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTTTAGGCAAAGGGAATTATCGCATAGTCAATGCGCTct                                                                                                                                    <  2:620496/149‑2 (MQ=255)                  gtAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATcc                                                                                                                               <  1:720913/150‑1 (MQ=255)                   tAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAAt                                                                                                                                           >  2:735024/1‑137 (MQ=255)                   tAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAAt                                                                                                                                           <  1:735073/137‑1 (MQ=255)                        aaTAAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAgg                                                                                                                        >  2:813892/1‑151 (MQ=255)                           cAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGAcc                                                                                                                     <  1:758296/151‑1 (MQ=255)                                  gagaTAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGaa                                                                                                               >  2:574591/1‑150 (MQ=255)                                  gagaTAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGaa                                                                                                               >  1:28379/1‑150 (MQ=255)                                       aaaaGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATc                                                                                                          <  1:702784/150‑1 (MQ=255)                                        aaaGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCAt                                                                                                        >  2:429186/1‑151 (MQ=255)                                           ggCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTaa                                                                                                     <  1:528387/151‑1 (MQ=255)                                                cATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCg                                                                                                 <  2:144275/150‑1 (MQ=255)                                                cATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGa                                                                                                <  2:870268/151‑1 (MQ=255)                                                   gTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATc                                                                                              <  2:269939/150‑1 (MQ=255)                                                     aaTGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGtt                                                                                           <  1:564454/151‑1 (MQ=255)                                                      aTGCAGGTTCACGACCAGCGCGAACGGCAGCATGATACAAAAAAGATAAGCGGTACGATGCAAAATCAGGGGGTAGGCAAAGGGAATTCGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTg                                                                                          >  2:131844/1‑151 (MQ=255)                                                       tGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTg                                                                                          >  1:585965/1‑150 (MQ=255)                                                           ggTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCt                                                                                     <  1:239007/151‑1 (MQ=255)                                                            gTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCta                                                                                    <  1:305133/151‑1 (MQ=255)                                                            gTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCta                                                                                    <  1:360580/151‑1 (MQ=255)                                                              ccACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCtaat                                                                                  >  2:398802/1‑151 (MQ=255)                                                               cccGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCtaata                                                                                 <  2:323234/149‑1 (MQ=255)                                                                aCGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCtaataa                                                                                <  1:679927/151‑1 (MQ=255)                                                                 cGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATaaa                                                                               <  1:629568/151‑1 (MQ=255)                                                                       gcgcGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGat                                                                         <  2:689/151‑1 (MQ=255)                                                                          cgAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATc                                                                      <  2:639152/151‑1 (MQ=255)                                                                          cgAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATc                                                                      >  2:693988/1‑151 (MQ=255)                                                                             aCGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTg                                                                    >  1:535337/1‑150 (MQ=255)                                                                              cGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGaa                                                                  <  1:580169/151‑1 (MQ=255)                                                                                  agcaTGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCt                                                              <  1:389044/151‑1 (MQ=255)                                                                                   gcaTGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTg                                                             <  2:667139/151‑1 (MQ=255)                                                                                     aTGATACAAAACAGATAAACGGTACGATGCAAATTCAGGGTGTAGGCAAAGGGAATTAGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACAGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCt                                                              >  2:851253/1‑148 (MQ=255)                                                                                      tGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTcc                                                          >  1:217555/1‑151 (MQ=255)                                                                                              aaCAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTgcggc                                                  <  2:410097/151‑1 (MQ=255)                                                                                              aaCAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTgcgg                                                   >  2:173881/1‑150 (MQ=255)                                                                                                 aGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCt                                               <  1:883084/151‑1 (MQ=255)                                                                                                 aGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTAAAAGCTGTCCATTGCGGCGCt                                               <  2:873377/151‑1 (MQ=255)                                                                                                   aTAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTg                                              <  2:557630/150‑1 (MQ=255)                                                                                                     aaaCGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAAc                                           >  2:701742/1‑151 (MQ=255)                                                                                                        cGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGc                                        >  2:410908/1‑151 (MQ=255)                                                                                                        cGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGc                                        <  2:762921/151‑1 (MQ=255)                                                                                                               aTGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAGCCGCCAACCAt                                 <  2:783407/151‑1 (MQ=255)                                                                                                               aTGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCAt                                 <  2:228900/151‑1 (MQ=255)                                                                                                                 gCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATtc                               >  2:19444/1‑151 (MQ=255)                                                                                                                  cAAAATCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATtct                              <  2:250316/151‑1 (MQ=255)                                                                                                                     aaTCAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTccc                           <  1:781679/151‑1 (MQ=255)                                                                                                                        cAGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCat                         <  2:528969/150‑1 (MQ=255)                                                                                                                         aGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCatta                       <  1:839777/151‑1 (MQ=255)                                                                                                                         aGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCatta                       <  1:492429/151‑1 (MQ=255)                                                                                                                         aGGGTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCatta                       <  1:564254/151‑1 (MQ=255)                                                                                                                           ggTGTAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCAGCCGGACAGGTCCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGTTATCTGAAAGCTGTCCATTGCGGGGCTGAACCGCCATCCATTCTCCCat                         <  2:561302/147‑1 (MQ=255)                                                                                                                            gtgtAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCAttatta                    <  1:549426/151‑1 (MQ=255)                                                                                                                              gtAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTaa                   <  2:121586/150‑1 (MQ=255)                                                                                                                              gtAGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAAc                  >  2:389080/1‑151 (MQ=255)                                                                                                                                aGGCAAAGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACaa                <  2:220402/151‑1 (MQ=255)                                                                                                                                    aaaGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACa                 >  2:478826/1‑146 (MQ=255)                                                                                                                                    aaaGGGAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGAt             >  2:53151/1‑150 (MQ=255)                                                                                                                                         gAATTGGCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGtt       >  1:326988/1‑151 (MQ=255)                                                                                                                                          aaTTGGCGCATAGGCAATGCGCTCTCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGtt       >  2:847050/1‑150 (MQ=255)                                                                                                                                            ttGGCGCATAGTCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAg     <  1:745411/150‑1 (MQ=255)                                                                                                                                              ggCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCt   <  1:748262/150‑1 (MQ=255)                                                                                                                                              ggCGCATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTg  <  1:905648/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TCCAGCGAAATAAAAGTGTAGGAAATCAGCACAGAGATAAAAGGCGTCATGTAATGCAGGTCCACGACCAGCGCGAACGGCAGCATGATACAAAACAGATAAACGGTACGATGCAAAATCAGGGTGTAGGCAAAGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTG  >  NC_000913/1608268‑1608559

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