Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	447,954	A→G	I31I (ATT→ATC)	cyoD ←	cytochrome o ubiquinol oxidase subunit IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	447,954	0	A	G	100.0%	250.8 / NA	66	I31I (ATT→ATC)	cyoD	cytochrome o ubiquinol oxidase subunit IV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/33);  total (33/33)

GCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAA  >  NC_000913/447808‑448092                                                                                                                                                   |                                                                                                                                           gCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATc                                                                                                                                         >  2:277080/1‑150 (MQ=255)  ccGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCAc                                                                                                                                       >  1:79705/1‑151 (MQ=255)   cGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCAcc                                                                                                                                      <  1:740574/151‑1 (MQ=255)   cGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCAcc                                                                                                                                      <  1:183273/151‑1 (MQ=255)    gTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCg                                                                                                                                     >  1:639761/1‑151 (MQ=255)         gTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGt                                                                                                                                    <  2:680036/147‑1 (MQ=255)         gTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTc                                                                                                                                   >  2:195611/1‑148 (MQ=255)         gTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGa                                                                                                                                <  2:235125/151‑1 (MQ=255)            ccAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAAt                                                                                                                              <  1:38666/150‑1 (MQ=255)              aGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGAt                                                                                                                           <  1:431432/151‑1 (MQ=255)               gCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATc                                                                                                                          <  1:442805/151‑1 (MQ=255)                    cATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAg                                                                                                                     <  1:227261/151‑1 (MQ=255)                      tCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAgg                                                                                                                    >  2:35461/1‑150 (MQ=255)                        tGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGAt                                                                                                                  >  1:226763/1‑150 (MQ=255)                        tGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGAt                                                                                                                  >  2:648107/1‑150 (MQ=255)                        tGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGAt                                                                                                                  >  1:45383/1‑150 (MQ=255)                         gATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATa                                                                                                                 <  2:280838/150‑1 (MQ=255)                           tttGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATaa                                                                                                                >  1:726726/1‑149 (MQ=255)                           tttGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATaa                                                                                                                >  1:722538/1‑149 (MQ=255)                                   tCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCAt                                                                                                      >  1:168541/1‑151 (MQ=255)                                    cATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCAt                                                                                                      <  2:669922/150‑1 (MQ=255)                                       gtgCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTa                                                                                                   <  1:199445/150‑1 (MQ=255)                                          cAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGtt                                                                                               <  1:226257/151‑1 (MQ=255)                                             gAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGttttt                                                                                            >  1:551701/1‑151 (MQ=255)                                               aGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTAc                                                                                          >  1:673923/1‑151 (MQ=255)                                               aGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTAc                                                                                          >  1:537643/1‑151 (MQ=255)                                                gCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACg                                                                                         <  2:643359/151‑1 (MQ=255)                                                 cacaCCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGc                                                                                        >  1:540097/1‑151 (MQ=255)                                                  acacCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCt                                                                                       >  1:245691/1‑151 (MQ=255)                                                  acacCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCt                                                                                       <  1:34242/151‑1 (MQ=255)                                                    acCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGc                                                                                     <  1:748175/151‑1 (MQ=255)                                                        gATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATg                                                                                 <  2:176181/151‑1 (MQ=255)                                                          tGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATggg                                                                               >  1:85121/1‑151 (MQ=255)                                                            cACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGAc                                                                             <  2:264046/151‑1 (MQ=255)                                                            cACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGAc                                                                             >  1:605810/1‑151 (MQ=255)                                                              ccAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGAcgc                                                                           >  1:94002/1‑151 (MQ=255)                                                              ccAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGAcgc                                                                           >  2:260282/1‑151 (MQ=255)                                                                    ccTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCt                                                                     >  1:17067/1‑151 (MQ=255)                                                                      tGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCtgt                                                                   >  2:196021/1‑151 (MQ=255)                                                                          ccACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATc                                                               <  2:663825/151‑1 (MQ=255)                                                                                ccATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCgg                                                             >  1:356703/1‑147 (MQ=255)                                                                                    tGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATg                                                      >  1:249506/1‑150 (MQ=255)                                                                                    tGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATg                                                      <  1:650968/150‑1 (MQ=255)                                                                                         ggATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCAt                                                <  2:611396/151‑1 (MQ=255)                                                                                          gATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCAtt                                               >  2:646904/1‑151 (MQ=255)                                                                                          gATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCAtt                                               >  2:567896/1‑151 (MQ=255)                                                                                                 cccAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCg                                         <  2:551044/150‑1 (MQ=255)                                                                                                 cccAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGc                                        >  1:434460/1‑151 (MQ=255)                                                                                                  ccAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGcc                                       >  2:708269/1‑151 (MQ=255)                                                                                                    aGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCAt                                                <  1:39337/140‑1 (MQ=255)                                                                                                    aGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCAt                                                >  2:39336/1‑140 (MQ=255)                                                                                                       aTTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCAt                                  <  1:513030/151‑1 (MQ=255)                                                                                                           cGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGa                              <  2:444478/151‑1 (MQ=255)                                                                                                             gCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATa                            <  1:567927/151‑1 (MQ=255)                                                                                                              ccGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATaa                           <  2:226754/151‑1 (MQ=255)                                                                                                                           tCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAac               >  1:130898/1‑150 (MQ=255)                                                                                                                           tCCTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAac               >  1:112505/1‑150 (MQ=255)                                                                                                                            ccTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAacac             <  1:238896/151‑1 (MQ=255)                                                                                                                            ccTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAaca              <  2:487419/150‑1 (MQ=255)                                                                                                                            ccTGTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAac               <  2:315962/149‑1 (MQ=255)                                                                                                                               gTCATCACCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAAcacaca          <  1:259598/151‑1 (MQ=255)                                                                                                                                 catcaCCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAaacaac                      <  2:649130/137‑1 (MQ=255)                                                                                                                                 catcaCCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAaacaac                      >  1:649169/1‑137 (MQ=255)                                                                                                                                 catcaCCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAacacac           >  2:753334/1‑148 (MQ=255)                                                                                                                                  atcaCCATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAg         <  2:661517/149‑1 (MQ=255)                                                                                                                                       cATCCAGAACGGGATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCaaa  <  2:207308/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                           GCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAA  >  NC_000913/447808‑448092

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