Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,454,292	C→A	G122G (GGC→GGA)	paaA →	ring 1,2‑phenylacetyl‑CoA epoxidase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,454,292	0	C	A	100.0%	239.5 / NA	74	G122G (GGC→GGA)	paaA	ring 1,2‑phenylacetyl‑CoA epoxidase subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (39/35);  total (39/35)

TATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGCTGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTTTGAAG  >  NC_000913/1454158‑1454419                                                                                                                                       |                                                                                                                                tATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCggatggc                                                                                                                             >  2:710135/1‑139 (MQ=255) tATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCggatggc                                                                                                                             >  1:1295622/1‑139 (MQ=255)        ataGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGa                                                                                                                      <  1:710135/139‑1 (MQ=255)         taGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGCGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGAt                                                                                                                     >  1:1178446/1‑139 (MQ=255)           gcgcCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGa                                                                                                                      <  1:1035641/136‑1 (MQ=255)           gcgcCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGa                                                                                                                      >  2:1035641/1‑136 (MQ=255)           gcgcCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATgg                                                                                                                   >  1:529658/1‑139 (MQ=255)           gcgcCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATgg                                                                                                                   >  2:367728/1‑139 (MQ=255)           gcgcCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACAGTGAGTTGTGCTGATAGCGGTGTTATCGGATGTCTGGTGGTTgg                                                                                                                   >  1:1195266/1‑139 (MQ=255)           gcgcCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGAAGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGATGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATgg                                                                                                                   >  2:302015/1‑139 (MQ=255)              ccgccgAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGgcgc                                                                                                                >  1:1180286/1‑139 (MQ=255)                 ccgAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGa                                                                                                                      <  1:860247/130‑1 (MQ=255)                 ccgAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGa                                                                                                                      >  2:860247/1‑130 (MQ=255)                  cgAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCg                                                                                                            <  2:1295622/139‑1 (MQ=255)                    aaaCACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGAt                                                                                                          <  1:367728/139‑1 (MQ=255)                     aaCACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATc                                                                                                         <  2:1271841/139‑1 (MQ=255)                           gggCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAAt                                                                                                   <  1:302015/139‑1 (MQ=255)                              cTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCaa                                                                                                    >  2:159506/1‑135 (MQ=255)                              cTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCaa                                                                                                    <  1:159506/135‑1 (MQ=255)                                 cgcCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTg                                                                                             >  2:735635/1‑139 (MQ=255)                                           gaCATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGcc                                                                                   >  1:770797/1‑139 (MQ=255)                                             cATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGc                                                                                 <  2:529658/139‑1 (MQ=255)                                                cTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCg                                                                                  <  2:544071/135‑1 (MQ=255)                                                cTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCg                                                                                  >  1:544071/1‑135 (MQ=255)                                                      aaaaaTGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTa                                                                                       >  1:524650/1‑124 (MQ=255)                                                      aaaaaTGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTa                                                                                       <  2:524650/124‑1 (MQ=255)                                                      aaaaaTGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTAt                                                                        >  1:1228177/1‑139 (MQ=255)                                                       aaaaTGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATg                                                                       <  1:413228/139‑1 (MQ=255)                                                            gCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCg                                                                  <  2:1228177/139‑1 (MQ=255)                                                              tCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTa                                                                <  2:1180286/139‑1 (MQ=255)                                                                                 ttCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGaaa                                             <  1:735635/139‑1 (MQ=255)                                                                                 ttCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGaaa                                             >  2:1020389/1‑139 (MQ=255)                                                                                  tcctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                               <  2:878837/104‑1 (MQ=255)                                                                                  tcctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCAc                                                                               >  1:878837/1‑104 (MQ=255)                                                                                   cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATc                                                                                                  >  1:962157/1‑84 (MQ=255)                                                                                   cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATc                                                                                                  <  2:962157/84‑1 (MQ=255)                                                                                   cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTa                                                                                       <  1:778437/95‑1 (MQ=255)                                                                                   cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTa                                                                                       >  2:778437/1‑95 (MQ=255)                                                                                   cctccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAAt                                           >  1:88498/1‑139 (MQ=255)                                                                                     tccATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCt                                         <  1:804743/139‑1 (MQ=255)                                                                                       cATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCg                                                                                          <  1:1070973/88‑1 (MQ=255)                                                                                       cATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCg                                                                                          >  2:1070973/1‑88 (MQ=255)                                                                                       cATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGc                                                             >  2:1287013/1‑117 (MQ=255)                                                                                       cATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGc                                                             <  1:1287013/117‑1 (MQ=255)                                                                                       cATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCa                                                      >  2:527613/1‑124 (MQ=255)                                                                                       cATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCa                                                      <  1:527613/124‑1 (MQ=255)                                                                                       cATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGc                                       >  1:1040142/1‑139 (MQ=255)                                                                                         tCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGcc                                                                                   <  1:197586/93‑1 (MQ=255)                                                                                         tCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGcc                                                                                   >  2:197586/1‑93 (MQ=255)                                                                                             tAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGc                                                                     <  2:483528/103‑1 (MQ=255)                                                                                             tAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGc                                                                     >  1:483528/1‑103 (MQ=255)                                                                                               aTTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAagaaga                               <  2:770797/139‑1 (MQ=255)                                                                                                   tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGc                                                                     <  2:1183061/97‑1 (MQ=255)                                                                                                   tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGc                                                                     >  1:1183061/1‑97 (MQ=255)                                                                                                   tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGt                           >  2:404005/1‑139 (MQ=255)                                                                                                   tCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGt                           >  1:1011291/1‑139 (MQ=255)                                                                                                             aGTTGGGCTGATATCGGTGTTATCGGATGGCTGGAGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGc                 >  1:82311/1‑139 (MQ=255)                                                                                                              gTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATg                                                                                     >  2:583975/1‑70 (MQ=255)                                                                                                              gTTGGGCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATg                                                                                     <  1:583975/70‑1 (MQ=255)                                                                                                                  ggCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGc                                                       >  2:512803/1‑96 (MQ=255)                                                                                                                  ggCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGc                                                       <  1:512803/96‑1 (MQ=255)                                                                                                                  ggCTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAg            >  2:671260/1‑139 (MQ=255)                                                                                                                    cTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCagag                                                        >  1:1331114/1‑93 (MQ=255)                                                                                                                    cTGATATCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCagag                                                        <  2:1331114/93‑1 (MQ=255)                                                                                                                       atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAAt                                                                                                   <  2:503097/47‑1 (MQ=255)                                                                                                                       atatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAAt                                                                                                   >  1:503097/1‑47 (MQ=255)                                                                                                                        tatCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGATGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGttt      <  2:1011291/139‑1 (MQ=255)                                                                                                                         atCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAgg                                                                                               <  2:615954/49‑1 (MQ=255)                                                                                                                         atCGGTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAgg                                                                                               >  1:615954/1‑49 (MQ=255)                                                                                                                            ggTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAaagaag                                >  1:1171025/1‑109 (MQ=255)                                                                                                                            ggTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAaagaag                                <  2:1171025/109‑1 (MQ=255)                                                                                                                            ggTGTTATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTTTCAAg  >  1:921660/1‑139 (MQ=255)                                                                                                                                 tATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGc                                       >  1:275800/1‑97 (MQ=255)                                                                                                                                 tATCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGc                                       <  2:275800/97‑1 (MQ=255)                                                                                                                                   tCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGc                                                             >  1:1227250/1‑73 (MQ=255)                                                                                                                                   tCGGATGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGc                                                             <  2:1227250/73‑1 (MQ=255)                                                                                                                                       |                                                                                                                                TATCTCTATAGCGCCGCCGAAACACTGGGCTGCGCCAGAGAAGACATCTACCAAAAAATGCTCGACGGACGGATGAAATATTCCTCCATCTTTAATTATCCGACATTGAGTTGGGCTGATATCGGTGTTATCGGCTGGCTGGTGGATGGCGCAGCGATCGTCAATCAGGTGGCGTTATGCCGCACTTCTTATGGCCCGTATGCCAGAGCAATGGTGAAAATCTGCAAAGAAGAGAGTTTTCACCAGCGTCAGGGGTTTGAAG  >  NC_000913/1454158‑1454419

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