Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,399,151	Δ1 bp	coding (1740/1848 nt)	mutL →	methyl‑directed mismatch repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,399,151	0	G	.	100.0%	348.3 / NA	72	coding (1740/1848 nt)	mutL	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (39/33);  total (39/33)

AACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGATTTTTTTG  >  NC_000913/4399006‑4399296                                                                                                                                                  |                                                                                                                                                  aacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCg                                                                                                                                             >  2:677862/1‑151 (MQ=255)  acaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCgg                                                                                                                                            >  1:559131/1‑151 (MQ=255)      aaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGt                                                                                                                                        <  1:697631/151‑1 (MQ=255)        tttACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTgg                                                                                                                                      >  1:30092/1‑151 (MQ=255)            cAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACg                                                                                                                                  >  1:342291/1‑151 (MQ=255)            cAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACg                                                                                                                                  >  2:62011/1‑151 (MQ=255)               aTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGtt                                                                                                                               >  1:358737/1‑151 (MQ=255)                 cTTGATTCCTGAACTGACAGGCTACCTGGCGAAGCAGTCCGCATTCTAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACTCT‑CTGGCGGACGTGGCACGGTTAt                                                                                                                             >  2:688173/1‑151 (MQ=255)                   tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTAt                                                                                                                             <  2:252819/149‑1 (MQ=255)                    gATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTc                                                                                                                          >  2:652567/1‑151 (MQ=255)                       tCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCg                                                                                                                        <  1:677896/150‑1 (MQ=255)                       tCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCg                                                                                                                        <  1:636191/150‑1 (MQ=255)                             aCTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTtg                                                                                                                 >  1:252842/1‑151 (MQ=255)                                gATAGGCTACCTGGCGAAGCAGTCCGTACTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCCCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTtgtg                                                                                                               <  2:103856/150‑1 (MQ=255)                                       tACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgcc                                                                                                       >  1:717996/1‑151 (MQ=255)                                          cTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcc                                                                                                    >  2:31730/1‑151 (MQ=255)                                           tGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGccgccg                                                                                                   >  1:461144/1‑151 (MQ=255)                                               gAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCGATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggt                                                                                                <  2:531856/150‑1 (MQ=255)                                                 aGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtggt                                                                                             >  1:910894/1‑151 (MQ=255)                                                    aGTCCGTATTCGAACCTGGCAACATTTCGCAGCGGATTGCACGCCATCTGATGAGCGAACATGCGCAGTGGCCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGATATGTCCGCAACTTGCGAAAACGCCGCCGGGTGGTCTg                                                                                          >  2:524844/1‑151 (MQ=255)                                                     gTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGt                                                                                         >  2:722316/1‑151 (MQ=255)                                                      tCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGtt                                                                                        >  1:756586/1‑151 (MQ=255)                                                       ccGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTa                                                                                       <  1:652602/151‑1 (MQ=255)                                                        cGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                                                                                      <  1:862058/151‑1 (MQ=255)                                                               gAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGt                                                                               <  1:741331/151‑1 (MQ=255)                                                                aaCCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGtt                                                                              >  1:558034/1‑151 (MQ=255)                                                                    tGGCAATATTGCGCAGTGGATTGCACGAAATCTGCTGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCAT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAt                                                                           >  2:478172/1‑150 (MQ=255)                                                                       cAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTAc                                                                       <  2:950549/151‑1 (MQ=255)                                                                        aaTATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTAc                                                                       >  1:729981/1‑150 (MQ=255)                                                                          tatTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATc                                                                    >  1:517561/1‑151 (MQ=255)                                                                              gcgcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGc                                                                <  2:138088/151‑1 (MQ=255)                                                                               cgcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCg                                                               >  1:585566/1‑151 (MQ=255)                                                                                gcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGa                                                              <  2:36408/151‑1 (MQ=255)                                                                                 cAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGAt                                                             >  1:64658/1‑151 (MQ=255)                                                                                 cAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGAt                                                             >  2:910041/1‑151 (MQ=255)                                                                                  aGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGAt                                                             >  2:74490/1‑150 (MQ=255)                                                                                  aGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATa                                                            <  2:220272/151‑1 (MQ=255)                                                                                   gTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATaa                                                           <  1:851697/151‑1 (MQ=255)                                                                                     ggATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATaaaa                                                         <  1:934873/151‑1 (MQ=255)                                                                                      gATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATaaa                                                          >  1:957500/1‑149 (MQ=255)                                                                                      gATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAg                                                        <  1:822089/151‑1 (MQ=255)                                                                                         tGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGcc                                                      <  1:644020/150‑1 (MQ=255)                                                                                           cACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTg                                                   <  2:47168/151‑1 (MQ=255)                                                                                            aCGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTg                                                   >  2:450198/1‑150 (MQ=255)                                                                                            aCGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGa                                                  <  2:541906/151‑1 (MQ=255)                                                                                              gAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGaaa                                                <  1:711509/151‑1 (MQ=255)                                                                                               aaaTCTGATGAGCGATCATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAg                                               <  2:675960/151‑1 (MQ=255)                                                                                               aaaTCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGcc                                                      >  2:764724/1‑144 (MQ=255)                                                                                               aaaTCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGcc                                                      <  1:764766/144‑1 (MQ=255)                                                                                               aaaTCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGaa                                                 >  2:399815/1‑149 (MQ=255)                                                                                               aaaTCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAg                                               <  1:460554/151‑1 (MQ=255)                                                                                               aaaTCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGGTATGTCCGCAAATTGTGAAAACGCCGCCGGGGGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGaa                                                 >  2:404664/1‑149 (MQ=255)                                                                                                    tgatgaGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATg                                            >  2:319515/1‑149 (MQ=255)                                                                                                     gatgaGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAg                                          <  2:276777/150‑1 (MQ=255)                                                                                                       tgaGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTg                                        >  1:478548/1‑150 (MQ=255)                                                                                                          gCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGatat                                    >  2:957695/1‑151 (MQ=255)                                                                                                          gCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGatat                                    <  1:506278/151‑1 (MQ=255)                                                                                                             aaCATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATc                                   >  1:728919/1‑149 (MQ=255)                                                                                                              aCATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGt                                >  1:227885/1‑151 (MQ=255)                                                                                                               cATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGt                                >  1:425765/1‑150 (MQ=255)                                                                                                               cATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTa                               <  1:327647/151‑1 (MQ=255)                                                                                                                   cgcAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGc                           >  2:36111/1‑151 (MQ=255)                                                                                                                      aGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCg                          >  1:877528/1‑149 (MQ=255)                                                                                                                         ggTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcc                      <  1:860358/150‑1 (MQ=255)                                                                                                                            cAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcctgc                   <  1:777560/150‑1 (MQ=255)                                                                                                                                 gCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAgg             <  1:948940/151‑1 (MQ=255)                                                                                                                                   acaGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCg           <  2:660083/151‑1 (MQ=255)                                                                                                                                   acaGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCg           <  1:722354/151‑1 (MQ=255)                                                                                                                                      ggCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGAtt        <  1:819300/151‑1 (MQ=255)                                                                                                                                       gCCATAACCCT‑CTGGCGGACGGGGAACGGTTAAGTCCGCAACTTGTGAAAACGCCGCCAGGGGGTCTGTTACAATCGGTTGATCTGCATACGGCGATGAAAGCCCTGAACGATGGGTGATACCAGTAAGGCAGGCCTGCCTAAGGGGAtat       >  2:1020865/1‑149 (MQ=255)                                                                                                                                            aaCCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGAttttttt   <  2:652394/150‑1 (MQ=255)                                                                                                                                            aaCCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGATTTTTTTg  >  2:493169/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                  AACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGATTTTTTTG  >  NC_000913/4399006‑4399296

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: ‑