Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,304,235	A→G	K111E (AAG→GAG)	eco →	ecotin, a serine protease inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,304,235	0	A	G	100.0%	222.8 / NA	61	K111E (AAG→GAG)	eco	ecotin, a serine protease inhibitor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/35);  total (26/35)

CAGACGCTGGAAGTCGATTGCAATTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCAAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCG  >  NC_000913/2304094‑2304381                                                                                                                                              |                                                                                                                                                   cAGACGCTGGAAGTCGATTGCAATTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAgaga                                                                                                                                           <  2:803500/151‑1 (MQ=255)   gACGCTGGAAGTCGATTGCAATTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGagaag                                                                                                                                         <  1:58773/151‑1 (MQ=255)         ggAAGTCGATTGCAATTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAAttt                                                                                                                                   <  1:3926/151‑1 (MQ=255)           aaGTCGATTGCAATTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGt                                                                                                                                 <  1:348561/151‑1 (MQ=255)               cGATTGCAATTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCAc                                                                                                                              <  2:33332/150‑1 (MQ=255)               cGATTGCAATTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCAc                                                                                                                              <  1:707056/150‑1 (MQ=255)                    gCAATTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGt                                                                                                                         <  1:176198/150‑1 (MQ=255)                     cAATTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTAt                                                                                                                       >  2:454602/1‑151 (MQ=255)                      aaTTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATc                                                                                                                      <  2:262568/151‑1 (MQ=255)                          tGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGGATCTggg                                                                                                                  >  1:508356/1‑151 (MQ=255)                            cATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCg                                                                                                                >  1:116742/1‑151 (MQ=255)                               cGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGAt                                                                                                              <  2:644316/150‑1 (MQ=255)                               cGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGGTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATg                                                                                                             >  1:188612/1‑151 (MQ=255)                                      gcggGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAAt                                                                                                      >  2:808798/1‑151 (MQ=255)                                       cggGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAAtg                                                                                                     >  1:75767/1‑151 (MQ=255)                                           aaGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAAtgctgc                                                                                                 <  2:346344/151‑1 (MQ=255)                                           aaGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAAtgctgc                                                                                                 >  1:473667/1‑151 (MQ=255)                                           aaGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAAtgctgc                                                                                                 <  1:536175/151‑1 (MQ=255)                                           aaGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAAtgctgc                                                                                                 >  1:12551/1‑151 (MQ=255)                                           aaGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAAtgctgc                                                                                                 <  1:120764/151‑1 (MQ=255)                                                    aacaaaacGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTacaaca                                                                                        <  2:559765/151‑1 (MQ=255)                                                    aacaaaacGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTacaac                                                                                         >  1:697154/1‑150 (MQ=255)                                                      caaaacGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGc                                                                                      <  1:210562/151‑1 (MQ=255)                                                      caaaacGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGc                                                                                      <  1:61609/151‑1 (MQ=255)                                                       aaaacGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCa                                                                                     <  2:762337/151‑1 (MQ=255)                                                         aacGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAg                                                                                   <  2:206377/151‑1 (MQ=255)                                                            gCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTg                                                                                <  2:76616/151‑1 (MQ=255)                                                            gCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTg                                                                                <  1:738021/151‑1 (MQ=255)                                                            gCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTg                                                                                >  1:94071/1‑151 (MQ=255)                                                               ggAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGcc                                                                              <  2:570765/150‑1 (MQ=255)                                                               ggAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCg                                                                             <  1:96062/151‑1 (MQ=255)                                                                           cTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGtat                                                                 >  1:588798/1‑151 (MQ=255)                                                                               gATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACg                                                              <  1:566171/150‑1 (MQ=255)                                                                               gATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGc                                                             <  1:735266/151‑1 (MQ=255)                                                                                  tattatGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGcc                                                            >  2:507844/1‑149 (MQ=255)                                                                                   attatGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGtt                                                                                              >  2:126004/1‑114 (MQ=255)                                                                                   attatGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGtt                                                                                              <  1:125981/114‑1 (MQ=255)                                                                                    ttatGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGAc                                                        >  2:135350/1‑151 (MQ=255)                                                                                    ttatGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGAc                                                        >  2:158020/1‑151 (MQ=255)                                                                                     tatGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGAc                                                        <  1:46066/150‑1 (MQ=255)                                                                                        gTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAAt                                                     <  1:524109/150‑1 (MQ=255)                                                                                            ttGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGt                                                   >  2:615810/1‑148 (MQ=255)                                                                                              gATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGAtg                                              <  2:579203/151‑1 (MQ=255)                                                                                              gATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGAt                                               <  2:419131/150‑1 (MQ=255)                                                                                                 aaaGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGAtgtg                                            >  1:597352/1‑150 (MQ=255)                                                                                                  aaGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGaa                                          <  2:154716/151‑1 (MQ=255)                                                                                                   aGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGa                                           >  1:388821/1‑149 (MQ=255)                                                                                                    gTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGt                                        >  2:221727/1‑151 (MQ=255)                                                                                                        gTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACcg                                    >  1:597657/1‑151 (MQ=255)                                                                                                          tCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACcgc                                   <  2:828640/150‑1 (MQ=255)                                                                                                           ccccGGTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACcgca                                  <  1:378108/150‑2 (MQ=255)                                                                                                                gTTTCAACGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTgg                             >  2:433101/1‑150 (MQ=255)                                                                                                                      aCGATGATGGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCgg                      >  1:431502/1‑151 (MQ=255)                                                                                                                          tgatgGCCTGCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGa                     >  2:389356/1‑148 (MQ=255)                                                                                                                                   gCCCGGATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAAtt           <  1:433068/149‑1 (MQ=255)                                                                                                                                       ggATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACa       <  1:780507/149‑1 (MQ=255)                                                                                                                                       ggATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACa       >  2:780563/1‑149 (MQ=255)                                                                                                                                        gATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAcg    >  1:87326/1‑151 (MQ=255)                                                                                                                                        gATGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAcg    <  1:223609/151‑1 (MQ=255)                                                                                                                                          tGGCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAcgcg  >  2:261673/1‑151 (MQ=255)                                                                                                                                           ggCGAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAAcgcg  <  2:785982/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   CAGACGCTGGAAGTCGATTGCAATTTGCATCGTCTCGGCGGGAAGCTGGAAAACAAAACGCTGGAAGGCTGGGGCTATGATTATTATGTCTTTGATAAAGTCAGTTCCCCGGTTTCAACGATGATGGCCTGCCCGGATGGCAAGAAAGAGAAGAAATTTGTCACCGCGTATCTGGGCGATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCG  >  NC_000913/2304094‑2304381

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