Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,787,434	G→A	pseudogene (1038/2253 nt)	ygaQ →	uncharacterized protein, putative enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,787,434	0	G	A	100.0%	267.1 / NA	81	pseudogene (1038/2253 nt)	ygaQ	uncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (42/39);  total (42/39)

GATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGG  >  NC_000913/2787303‑2787566                                                                                                                                    |                                                                                                                                     gATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTc                                                                                                                               <  1:729473/139‑1 (MQ=255)         tCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTg                                                                                                                       >  1:114555/1‑139 (MQ=255)                   aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattca                                                                                                             >  2:102816/1‑139 (MQ=255)                   aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattca                                                                                                             >  2:289820/1‑139 (MQ=255)                     cTTTAACTTGATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattcatt                                                                                                           <  1:920266/139‑1 (MQ=255)                       ttAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTaa                                                                                                         <  1:289820/139‑1 (MQ=255)                        tAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTaaa                                                                                                        >  1:793751/1‑139 (MQ=255)                               aTGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAg                                                                                                 >  1:784833/1‑139 (MQ=255)                                                           aCCTGGATTATCCCTCTGCTGATATACGAGCGAATATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAActccatcc                                                                     <  2:697470/139‑1 (MQ=255)                                                           aCCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     <  2:567010/139‑1 (MQ=255)                                                           aCCTGGATTATCCCTCTGCTGATAAACGAGCGAAGATCAGAAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     >  1:567010/1‑139 (MQ=255)                                                            ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTTAGGTTTTGATTCcattc                                                                                                              <  2:533182/97‑1 (MQ=255)                                                            ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattc                                                                                                              >  1:533182/1‑97 (MQ=255)                                                              tGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                  <  1:321239/139‑1 (MQ=255)                                                                  ttATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGat                                                              <  1:102816/139‑1 (MQ=255)                                                                     tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattcatt                                                                                                           >  2:318337/1‑91 (MQ=255)                                                                     tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattcatt                                                                                                           <  1:318337/91‑1 (MQ=255)                                                                     tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTa                                                                                                          <  2:326510/92‑1 (MQ=255)                                                                     tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTa                                                                                                          >  1:326510/1‑92 (MQ=255)                                                                     tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     >  1:29500/1‑129 (MQ=255)                                                                     tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     <  2:29500/129‑1 (MQ=255)                                                                          ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                         >  1:112844/1‑72 (MQ=255)                                                                          ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                         <  2:112844/72‑1 (MQ=255)                                                                          ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                   >  2:513210/1‑94 (MQ=255)                                                                          ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                   <  1:513210/94‑1 (MQ=255)                                                                          ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                           >  2:105157/1‑118 (MQ=255)                                                                          ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                           <  1:105157/118‑1 (MQ=255)                                                                          ctactgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                      >  2:539232/4‑139 (MQ=255)                                                                           tgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtctc                                                     >  2:543855/1‑139 (MQ=255)                                                                               gATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                 <  2:367362/75‑1 (MQ=255)                                                                               gATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                 >  1:367362/1‑75 (MQ=255)                                                                               gATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                    >  1:849583/1‑120 (MQ=255)                                                                               gATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                    <  2:849583/120‑1 (MQ=255)                                                                                atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                  <  1:75104/121‑1 (MQ=255)                                                                                atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                  >  2:75104/1‑121 (MQ=255)                                                                                       aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                 <  1:100867/67‑1 (MQ=255)                                                                                       aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                                 >  2:100867/1‑67 (MQ=255)                                                                                       aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtc                                                       <  1:189935/125‑1 (MQ=255)                                                                                       aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtc                                                       >  2:189935/1‑125 (MQ=255)                                                                                       aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAt                                                        <  1:70141/124‑1 (MQ=255)                                                                                       aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAt                                                        >  2:70141/1‑124 (MQ=255)                                                                                       aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGAc                                         >  1:111329/1‑139 (MQ=255)                                                                                          gAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     <  1:153200/108‑1 (MQ=255)                                                                                          gAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     >  2:153200/1‑108 (MQ=255)                                                                                          gAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTc                                      >  1:707623/1‑139 (MQ=255)                                                                                              aTCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                                  <  2:793751/139‑1 (MQ=255)                                                                                                  gTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGTACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAg                              >  2:577979/1‑139 (MQ=255)                                                                                                   tAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGt                             <  1:854049/139‑1 (MQ=255)                                                                                                        atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGa                        >  2:264019/1‑139 (MQ=255)                                                                                                         tataGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     >  1:844630/1‑93 (MQ=255)                                                                                                         tataGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     <  2:844630/93‑1 (MQ=255)                                                                                                            aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                    <  2:784833/139‑1 (MQ=255)                                                                                                             gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                      >  1:474359/1‑104 (MQ=255)                                                                                                             gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                      <  2:474359/104‑1 (MQ=255)                                                                                                             gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                    <  2:797725/138‑1 (MQ=255)                                                                                                             gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                    >  1:797725/1‑138 (MQ=255)                                                                                                                ggtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGTACATAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACAtct                <  1:577979/139‑1 (MQ=255)                                                                                                                 gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     >  2:647161/1‑85 (MQ=255)                                                                                                                 gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     <  1:647161/85‑1 (MQ=255)                                                                                                                  tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                   <  2:43772/54‑1 (MQ=255)                                                                                                                  tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                   <  2:425733/54‑1 (MQ=255)                                                                                                                  tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                   >  1:425733/1‑54 (MQ=255)                                                                                                                  tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                   >  1:43772/1‑54 (MQ=255)                                                                                                                  tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAgtccatcc                                                                     <  2:609769/84‑1 (MQ=255)                                                                                                                  tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                     >  1:609769/1‑84 (MQ=255)                                                                                                                  tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                           >  1:84248/1‑78 (MQ=255)                                                                                                                  tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                           <  2:84248/78‑1 (MQ=255)                                                                                                                  tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                                                                   <  1:292188/86‑1 (MQ=255)                                                                                                                  tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                                                                   >  2:292188/1‑86 (MQ=255)                                                                                                                   ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGa                        >  2:422205/1‑128 (MQ=255)                                                                                                                   ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGa                        <  1:422205/128‑1 (MQ=255)                                                                                                                   ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTa             >  1:378385/1‑139 (MQ=255)                                                                                                                     tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACg                                                                                    >  1:769770/1‑66 (MQ=255)                                                                                                                     tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACg                                                                                    <  2:769770/66‑1 (MQ=255)                                                                                                                     tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTa                            >  2:203206/1‑122 (MQ=255)                                                                                                                     tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTa                            <  1:203206/122‑1 (MQ=255)                                                                                                                            ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCa                                                    >  1:494019/1‑91 (MQ=255)                                                                                                                            ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCa                                                    <  2:494019/91‑1 (MQ=255)                                                                                                                            ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                                  <  1:401418/109‑1 (MQ=255)                                                                                                                            ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                                  >  2:401418/1‑109 (MQ=255)                                                                                                                              cTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATgg  >  2:822671/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                     GATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAATGG  >  NC_000913/2787303‑2787566

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