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

ACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGG  >  NC_000913/4399007‑4399285                                                                                                                                                 |                                                                                                                                       acaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCg                                                                                                                                  <  1:623946/150‑1 (MQ=255)       tttACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTgg                                                                                                                           >  1:916150/1‑151 (MQ=255)         tACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGa                                                                                                                          >  1:167164/1‑150 (MQ=255)           cAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACg                                                                                                                       <  2:132580/151‑1 (MQ=255)            aaaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACgg                                                                                                                      >  1:488930/1‑151 (MQ=255)            aaaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACgg                                                                                                                      <  1:375282/151‑1 (MQ=255)             aaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGt                                                                                                                     >  2:178281/1‑151 (MQ=255)             aaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGt                                                                                                                     <  2:995235/151‑1 (MQ=255)              aTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGtt                                                                                                                    >  1:24156/1‑151 (MQ=255)                  tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATg                                                                                                                 <  2:432399/150‑1 (MQ=255)                  tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGt                                                                                                                >  2:513680/1‑151 (MQ=255)                   gATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTc                                                                                                               >  1:559373/1‑151 (MQ=255)                    aTTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTcc                                                                                                              >  1:1044284/1‑151 (MQ=255)                     ttCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCg                                                                                                             >  2:715296/1‑151 (MQ=255)                             cTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTtg                                                                                                      <  2:733856/150‑1 (MQ=255)                               gATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGa                                                                                                   <  2:716145/151‑1 (MQ=255)                                 tAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaa                                                                                                  >  2:806427/1‑150 (MQ=255)                                   ggCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaaa                                                                                                <  2:987224/150‑1 (MQ=255)                                   ggCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAc                                                                                               <  2:1056830/151‑1 (MQ=255)                                     cTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgc                                                                                             >  2:515401/1‑151 (MQ=255)                                       aCCTGGCGAAGCAGTCCGTATTCGAACCTGGTAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgc                                                                                             <  1:774423/149‑1 (MQ=255)                                       aCCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccg                                                                                           >  1:96250/1‑151 (MQ=255)                                        ccTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccgc                                                                                          >  2:678627/1‑151 (MQ=255)                                         cTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcc                                                                                         >  1:535788/1‑151 (MQ=255)                                         cTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcc                                                                                         >  2:243408/1‑151 (MQ=255)                                           ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcc                                                                                         <  1:521333/149‑1 (MQ=255)                                             cGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggt                                                                                     >  1:704085/1‑151 (MQ=255)                                                      ccGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGtt                                                                             >  2:609462/1‑150 (MQ=255)                                                       cGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                                                                           <  2:997895/151‑1 (MQ=255)                                                       cGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                                                                           >  1:752527/1‑151 (MQ=255)                                                          aTTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAAt                                                                        <  1:765340/151‑1 (MQ=255)                                                            tCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATc                                                                       <  1:513694/150‑1 (MQ=255)                                                            tCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATc                                                                       <  1:16756/150‑1 (MQ=255)                                                            tCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCt                                                                      >  1:814756/1‑151 (MQ=255)                                                             cGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTg                                                                     >  2:195100/1‑151 (MQ=255)                                                              gAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGt                                                                    >  1:705951/1‑151 (MQ=255)                                                              gAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTATGAAAACGCCGCCGGGTGGTCTGTTACAATCTGt                                                                    >  2:685896/1‑151 (MQ=255)                                                               aaCCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGtt                                                                   >  2:484700/1‑151 (MQ=255)                                                               aaCCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGtt                                                                   >  1:326990/1‑151 (MQ=255)                                                                aCCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg                                                                  >  2:539949/1‑151 (MQ=255)                                                                    ggCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAtt                                                               <  1:1003332/150‑1 (MQ=255)                                                                     gCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAttt                                                              >  2:1050578/1‑150 (MQ=255)                                                                       aaTATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTAc                                                            >  1:89627/1‑150 (MQ=255)                                                                        atatTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACAt                                                          >  1:835527/1‑151 (MQ=255)                                                                        atatTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACAt                                                          <  1:789312/151‑1 (MQ=255)                                                                          atTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATcc                                                        >  2:1018212/1‑151 (MQ=255)                                                                            tGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCgg                                                      >  1:666457/1‑151 (MQ=255)                                                                            tGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCg                                                       <  1:615439/150‑1 (MQ=255)                                                                              cgcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCg                                                    >  1:935532/1‑151 (MQ=255)                                                                              cgcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCg                                                    >  1:167392/1‑151 (MQ=255)                                                                              cgcAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCg                                                    >  1:184763/1‑151 (MQ=255)                                                                                  gTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATaa                                                <  1:328600/151‑1 (MQ=255)                                                                                   tGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATaa                                                >  2:834267/1‑150 (MQ=255)                                                                                   tGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATCTCCGCAACTTGTGAAAACGCCGCCGGGTTGTCTGTTACAATCTGTTGATTTACATCCGGCGATaa                                                >  2:191325/1‑150 (MQ=255)                                                                                                aTCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGaa                                      >  1:314135/1‑147 (MQ=255)                                                                                                    gatgaGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGt                              >  1:601993/1‑151 (MQ=255)                                                                                                     atgaGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTg                             >  1:261058/1‑151 (MQ=255)                                                                                                       gaGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGat                           >  2:274890/1‑151 (MQ=255)                                                                                                        aGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGat                           >  2:655614/1‑150 (MQ=255)                                                                                                            aaCATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAg                      <  1:452378/151‑1 (MQ=255)                                                                                                             aCATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGt                     <  1:305150/151‑1 (MQ=255)                                                                                                             aCATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGt                     <  1:1005440/151‑1 (MQ=255)                                                                                                              cATGCGCAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTa                    <  2:387197/151‑1 (MQ=255)                                                                                                                 gcgcAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAgg                 >  1:223800/1‑151 (MQ=255)                                                                                                                   gcAGTGGTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCg               >  1:917903/1‑151 (MQ=255)                                                                                                                        ggTCAATGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAgcct          >  2:661580/1‑151 (MQ=255)                                                                                                                             aTGGCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTa     <  1:412250/151‑1 (MQ=255)                                                                                                                                gCACAGGCCATAACCCT‑CTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAgg  <  1:137025/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                       ACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCCGGCGATAAAAGCCCTGAAAGATGAGTGATATCAGTAAGGCGAGCCTGCCTAAGG  >  NC_000913/4399007‑4399285

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