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%	323.2 / NA	66	coding (1740/1848 nt)	mutL	methyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (41/25);  total (41/25)

AACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGATTTTTTT  >  NC_000913/4399006‑4399295                                                                                                                                                  |                                                                                                                                                 aacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCg                                                                                                                                            >  2:633020/1‑151 (MQ=255)   caaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGa                                                                                                                                          <  2:509533/151‑1 (MQ=255)         ttACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTgg                                                                                                                                     >  2:336286/1‑150 (MQ=255)          tACCAATCTTGGTTCCTGAACTGATAGGCTACCTGGCGATGCAGTCAGAATTCGAACCTCGCAATATTGCGCAGTGGATAGCACGAAATCTGATAAGCGAACATGCGCAGCGGTCAATGGCACAGGCCATAACTCT‑CTGGCGGAc                                                                                                                                         >  2:509679/1‑145 (MQ=255)          tACCAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTg                                                                                                                                      >  2:759633/1‑148 (MQ=255)           aCAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGTAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGAGGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCGT‑CTGGCGGACGTGGaa                                                                                                                                   <  2:727359/150‑1 (MQ=255)              aaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGt                                                                                                                               >  1:489031/1‑151 (MQ=255)                tCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGtt                                                                                                                              >  2:276083/1‑150 (MQ=255)                   tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGtt                                                                                                                              <  1:327530/147‑1 (MQ=255)                   tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATg                                                                                                                           <  2:216363/150‑1 (MQ=255)                     aTTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTcc                                                                                                                        <  1:161592/151‑1 (MQ=255)                      ttCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCg                                                                                                                       >  1:426569/1‑151 (MQ=255)                           gAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACt                                                                                                                  >  1:447030/1‑151 (MQ=255)                                gATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGa                                                                                                             <  2:673336/151‑1 (MQ=255)                                 aTAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGa                                                                                                             >  2:554972/1‑150 (MQ=255)                                   aGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaa                                                                                                           >  2:318525/1‑150 (MQ=255)                                    ggCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATGACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACGTGTGaaaa                                                                                                          <  1:296702/150‑1 (MQ=255)                                     gCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTACGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcg                                                                                                        <  1:218211/151‑1 (MQ=255)                                       tACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgcc                                                                                                      >  2:350396/1‑151 (MQ=255)                                        aCCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccg                                                                                                     <  1:24851/151‑1 (MQ=255)                                          cTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcc                                                                                                   >  1:751457/1‑151 (MQ=255)                                           tGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGccgccg                                                                                                  >  1:520784/1‑151 (MQ=255)                                            ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGccgccg                                                                                                  <  1:472216/150‑1 (MQ=255)                                            ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCgg                                                                                                 >  2:707653/1‑151 (MQ=255)                                            ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCgg                                                                                                 >  1:528383/1‑151 (MQ=255)                                               gAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtg                                                                                              >  2:687225/1‑151 (MQ=255)                                               gAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtg                                                                                              >  2:99702/1‑151 (MQ=255)                                                aaGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtgg                                                                                             <  2:109481/151‑1 (MQ=255)                                                        cGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                                                                                     >  1:740988/1‑151 (MQ=255)                                                         gTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                                                                                     >  1:98148/1‑150 (MQ=255)                                                           aTTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACaa                                                                                   <  1:501345/150‑1 (MQ=255)                                                             tCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCt                                                                                >  1:263307/1‑151 (MQ=255)                                                                aaCCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGtt                                                                             >  1:115524/1‑151 (MQ=255)                                                                      gCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAttt                                                                        >  2:761460/1‑150 (MQ=255)                                                                      gCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAttt                                                                        >  1:335168/1‑150 (MQ=255)                                                                         atatTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACAt                                                                    <  1:765772/151‑1 (MQ=255)                                                                           atTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATc                                                                   <  2:495749/150‑1 (MQ=255)                                                                              gcgcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCgg                                                                <  1:637629/150‑1 (MQ=255)                                                                               cgcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCg                                                              >  1:26191/1‑151 (MQ=255)                                                                                 cAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGAt                                                            <  1:382619/151‑1 (MQ=255)                                                                                   gTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATa                                                           >  1:522028/1‑150 (MQ=255)                                                                                      gATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAg                                                       <  1:358545/151‑1 (MQ=255)                                                                                       aTTGCACGAAAACTGATAAGCGAACATGCGAAGAGGTCAATGGCACAGGCCATAAACCT‑CTGGCGGACGTGGTAAGGGTATGTCCGAAACTTGTGAAATCGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCAGGCGATaaaa                                                        >  2:718209/1‑149 (MQ=255)                                                                                          gCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCt                                                   <  2:124598/151‑1 (MQ=255)                                                                                           cACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTg                                                  <  2:524723/151‑1 (MQ=255)                                                                                            aCGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTg                                                  >  2:684605/1‑150 (MQ=255)                                                                                             cGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGa                                                 >  2:586269/1‑150 (MQ=255)                                                                                                    tgatgaGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATg                                           >  1:2131/1‑149 (MQ=255)                                                                                                    tgatgaGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATg                                           >  2:129927/1‑149 (MQ=255)                                                                                                        gaGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGa                                      <  2:281730/150‑1 (MQ=255)                                                                                                           cGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATc                                  <  1:399306/151‑1 (MQ=255)                                                                                                              aCATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGt                               <  1:276094/151‑1 (MQ=255)                                                                                                                  gcgcAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAg                            <  1:644184/150‑1 (MQ=255)                                                                                                                    gcAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCg                         <  2:195059/151‑1 (MQ=255)                                                                                                                     cAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCg                         >  2:737339/1‑150 (MQ=255)                                                                                                                      aGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCg                         >  2:117795/1‑149 (MQ=255)                                                                                                                       gTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAGGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAg                       <  2:447009/150‑1 (MQ=255)                                                                                                                            cAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATTTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcctgcc                 >  1:225307/1‑151 (MQ=255)                                                                                                                            cAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcctgcc                 >  2:599671/1‑151 (MQ=255)                                                                                                                            cAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcctgcc                 >  2:523190/1‑151 (MQ=255)                                                                                                                               tGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTa               >  1:445351/1‑150 (MQ=255)                                                                                                                                  cacaGGCCATAACCCT‑CTGGCGGATGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAgg            >  2:686026/1‑150 (MQ=255)                                                                                                                                    caGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGTTATCAGTATGGCGAGCCTGCCTAAGGc           >  2:140368/1‑149 (MQ=255)                                                                                                                                     aGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGAt        >  1:509308/1‑151 (MQ=255)                                                                                                                                        ccATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGAtttt     >  2:380897/1‑151 (MQ=255)                                                                                                                                           tAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGAttttttt  >  2:493271/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                 AACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGGCGATTTTTTT  >  NC_000913/4399006‑4399295

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