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%	286.9 / NA	86	pseudogene (1038/2253 nt)	ygaQ	uncharacterized protein, putative enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (48/38);  total (48/38)

GATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAT  >  NC_000913/2787315‑2787564                                                                                                                        |                                                                                                                                   gATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAGATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTc                                                                                                                 >  2:1227599/1‑139 (MQ=255)  aTTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTcc                                                                                                                <  2:1075372/139‑1 (MQ=255)       aCCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattca                                                                                                           >  2:691677/1‑139 (MQ=255)         cTTTAACTTTTTGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATGCGAGCGAAGATAAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                       <  1:1729211/125‑1 (MQ=255)         cTTTAACTTTATGGGAATTGGGGGTAGTGGCAAAAATGACCTGGATTATCCCTCTGCTGATATGCGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                       >  2:1729211/1‑125 (MQ=255)                   aTGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAg                                                                                               >  1:138299/1‑139 (MQ=255)                    tGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGt                                                                                              >  1:1163309/1‑139 (MQ=255)                     gggAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTa                                                                                             <  1:1227599/139‑1 (MQ=255)                          ttGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCa                                                                                                                            >  1:757347/1‑103 (MQ=255)                          ttGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCa                                                                                                                            <  2:757347/103‑1 (MQ=255)                              gggTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATCTACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGaa                                                                                    >  1:1155261/1‑139 (MQ=255)                               ggTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                   <  2:1756517/139‑1 (MQ=255)                                  aGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                 >  2:1316900/1‑122 (MQ=255)                                  aGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                 <  1:1316900/122‑1 (MQ=255)                                   gTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                   <  1:442208/135‑1 (MQ=255)                                   gTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                   >  2:442208/1‑135 (MQ=255)                                    tGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTaaa                                                                              >  1:749165/1‑139 (MQ=255)                                         aaaaTGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                         <  1:691677/139‑1 (MQ=255)                                              gACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATc                                                                                                  >  2:1802219/1‑109 (MQ=255)                                              gACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATc                                                                                                  <  1:1802219/109‑1 (MQ=255)                                              gACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                    >  2:537114/1‑139 (MQ=255)                                              gACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                    >  1:116144/1‑139 (MQ=255)                                               aCCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                   <  2:1155261/139‑1 (MQ=255)                                                    gATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                       >  2:1303658/1‑82 (MQ=255)                                                    gATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGttt                                                                                                                       <  1:1303658/82‑1 (MQ=255)                                                         tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                >  2:1650059/1‑132 (MQ=255)                                                         tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                <  1:1650059/132‑1 (MQ=255)                                                         tCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGAtata                                                         >  2:380209/1‑139 (MQ=255)                                                              ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                    >  2:693250/1‑139 (MQ=255)                                                              ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCACTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                    >  1:951347/1‑139 (MQ=255)                                                               tgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                              <  1:434307/64‑1 (MQ=255)                                                               tgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                              >  2:434307/1‑64 (MQ=255)                                                                 ctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                   <  2:1248879/105‑1 (MQ=255)                                                                 ctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                   >  1:1248879/1‑105 (MQ=255)                                                                    atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTaaa                                                                                                      >  2:1405299/1‑83 (MQ=255)                                                                    atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTaaa                                                                                                      <  1:1405299/83‑1 (MQ=255)                                                                    atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                   <  1:475181/118‑1 (MQ=255)                                                                    atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                   >  2:475181/1‑118 (MQ=255)                                                                         cGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGCTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTCCAAAAg                                         >  2:797953/1‑139 (MQ=255)                                                                           aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                              >  1:133834/1‑52 (MQ=255)                                                                           aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                              <  2:133834/52‑1 (MQ=255)                                                                           aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGt                                                                                                                         >  1:1442383/1‑57 (MQ=255)                                                                           aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGt                                                                                                                         <  2:1442383/57‑1 (MQ=255)                                                                           aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                >  1:1427948/1‑114 (MQ=255)                                                                           aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                <  1:1233951/114‑1 (MQ=255)                                                                           aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                >  2:1233951/1‑114 (MQ=255)                                                                           aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                                <  2:1427948/114‑1 (MQ=255)                                                                                  aTCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                                <  2:749165/139‑1 (MQ=255)                                                                                     aGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCaa                             >  1:52303/1‑139 (MQ=255)                                                                                      gTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAg                            >  2:36488/1‑139 (MQ=255)                                                                                            atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattcatt                                                                                                         >  1:1520507/1‑56 (MQ=255)                                                                                            atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattcatt                                                                                                         <  2:1520507/56‑1 (MQ=255)                                                                                            atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                 <  1:72878/64‑1 (MQ=255)                                                                                            atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                 >  2:72878/1‑64 (MQ=255)                                                                                                aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  <  2:52303/139‑1 (MQ=255)                                                                                                aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  <  1:380209/139‑1 (MQ=255)                                                                                                aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  <  2:116144/139‑1 (MQ=255)                                                                                                      tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                   >  2:1595499/1‑68 (MQ=255)                                                                                                      tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAAc                                                                                   <  1:1595499/68‑1 (MQ=255)                                                                                                      tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                              <  1:1638981/89‑1 (MQ=255)                                                                                                      tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaag                                                              >  2:1638981/1‑89 (MQ=255)                                                                                                      tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTa                          <  2:1515109/125‑1 (MQ=255)                                                                                                      tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTa                          >  1:1515109/1‑125 (MQ=255)                                                                                                      tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATctct            >  2:1789760/1‑139 (MQ=255)                                                                                                       ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTa           >  1:1200376/1‑139 (MQ=255)                                                                                                        gtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAAGAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTaa          >  1:1078595/1‑139 (MQ=255)                                                                                                        gtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTaa          >  1:612619/1‑139 (MQ=255)                                                                                                        gtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTaa          >  2:940688/1‑139 (MQ=255)                                                                                                        gtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAG‑AGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAg         >  2:98951/1‑139 (MQ=255)                                                                                                         tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  >  2:152718/1‑130 (MQ=255)                                                                                                         tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                  <  1:152718/130‑1 (MQ=255)                                                                                                           tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACAt                <  2:235907/130‑1 (MQ=255)                                                                                                           tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACAt                >  1:235907/1‑130 (MQ=255)                                                                                                                ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                 >  2:1623163/1‑44 (MQ=255)                                                                                                                ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                 <  1:1623163/44‑1 (MQ=255)                                                                                                                ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                        >  2:227461/1‑53 (MQ=255)                                                                                                                ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACaaa                                                                                        <  1:227461/53‑1 (MQ=255)                                                                                                                ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                   >  1:675606/1‑74 (MQ=255)                                                                                                                ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                   <  2:675606/74‑1 (MQ=255)                                                                                                                ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                    <  2:908727/89‑1 (MQ=255)                                                                                                                ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACAtct                                                    >  1:908727/1‑89 (MQ=255)                                                                                                                ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAt  <  1:537114/139‑1 (MQ=255)                                                                                                                ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAt  >  1:357191/1‑139 (MQ=255)                                                                                                                   tATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGt                           <  2:1013197/111‑1 (MQ=255)                                                                                                                   tATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGt                           >  1:1013197/1‑111 (MQ=255)                                                                                                                     ttAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTaaa                        <  1:1490028/112‑1 (MQ=255)                                                                                                                        |                                                                                                                                   GATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAAT  >  NC_000913/2787315‑2787564

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