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

ATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAA  >  NC_000913/2787304‑2787563                                                                                                                                   |                                                                                                                                  atatCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCa                                                                                                                           <  2:238408/139‑1 (MQ=255)    tCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGt                                                                                                                        <  1:388798/139‑1 (MQ=255)           tGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAGAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGAtt                                                                                                                 <  1:1262648/139‑1 (MQ=255)           tGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGGCCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGt                                                                                                                        <  2:852170/132‑1 (MQ=255)           tGATTATACCTTTAACTTAATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGt                                                                                                                        >  1:852170/1‑132 (MQ=255)            gATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTc                                                                                                                >  1:141004/1‑139 (MQ=255)              ttATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                             <  1:540550/124‑1 (MQ=255)              ttATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGt                                                                                                                             >  2:540550/1‑124 (MQ=255)                    cTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattcatt                                                                                                        >  1:886824/1‑139 (MQ=255)                                              gTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                <  2:85631/121‑1 (MQ=255)                                              gTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                >  1:85631/1‑121 (MQ=255)                                               tGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattc                                                                                                           <  2:1192420/109‑1 (MQ=255)                                               tGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcattc                                                                                                           >  1:1192420/1‑109 (MQ=255)                                                  tAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACaa                                                                          <  1:1016719/139‑1 (MQ=255)                                                          aCCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                  <  1:312926/139‑1 (MQ=255)                                                           ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                >  1:811495/1‑108 (MQ=255)                                                           ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                <  2:811495/108‑1 (MQ=255)                                                           ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                 >  2:400950/1‑139 (MQ=255)                                                           ccTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAAtccatcca                                                                 >  2:919964/1‑139 (MQ=255)                                                            cTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                                                                <  1:400950/139‑1 (MQ=255)                                                             tGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                               <  1:919964/139‑1 (MQ=255)                                                                     cccTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATAc                                                       >  1:279471/1‑139 (MQ=255)                                                                         ctgctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATAAAtct                                                   >  1:1039031/1‑139 (MQ=255)                                                                            ctgATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCAc                                                >  1:1074140/1‑139 (MQ=255)                                                                               atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaa                                                              <  1:295967/122‑1 (MQ=255)                                                                               atatACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCagaa                                                              >  2:295967/1‑122 (MQ=255)                                                                                      aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCa                                                                                                                           <  2:104503/54‑1 (MQ=255)                                                                                      aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCa                                                                                                                           >  1:104503/1‑54 (MQ=255)                                                                                      aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                                                                <  1:548052/113‑1 (MQ=255)                                                                                      aGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCag                                                                >  2:548052/1‑113 (MQ=255)                                                                                        cGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACtt                                    <  1:110639/139‑1 (MQ=255)                                                                                         gAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGATATTAAGAAGGTCAGGTTTTGAGTCCATTCATTAAAGATCCAGTAACAAAAGAACGGAAACAAATCCATCCAGAAGATATACATATAACTGCAAAACACtta                                   >  2:581729/1‑138 (MQ=255)                                                                                           aGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGa                                 <  1:128858/139‑1 (MQ=255)                                                                                             aTCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAg                               <  2:886824/139‑1 (MQ=255)                                                                                                 gTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccat                                                                    >  2:1231155/1‑98 (MQ=255)                                                                                                 gTAAAGATATAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccat                                                                    <  1:1231155/98‑1 (MQ=255)                                                                                                       atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                              >  1:174989/1‑50 (MQ=255)                                                                                                       atatAGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCca                                                                                                              <  2:174989/50‑1 (MQ=255)                                                                                                         ataGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATa                   >  1:641436/1‑139 (MQ=255)                                                                                                         ataGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATa                   >  2:631077/1‑139 (MQ=255)                                                                                                          taGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                <  1:1078066/61‑1 (MQ=255)                                                                                                          taGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                >  2:1078066/1‑61 (MQ=255)                                                                                                           aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                   <  1:544215/89‑1 (MQ=255)                                                                                                           aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                   >  2:544215/1‑89 (MQ=255)                                                                                                           aGAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAAc                 <  2:141004/139‑1 (MQ=255)                                                                                                            gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaa                                          >  2:136889/1‑113 (MQ=255)                                                                                                            gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCaaa                                          <  1:136889/113‑1 (MQ=255)                                                                                                            gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCaa                            >  2:301592/1‑127 (MQ=255)                                                                                                            gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCaa                            <  1:301592/127‑1 (MQ=255)                                                                                                            gAAGGTGGTATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACa                >  1:86759/1‑139 (MQ=255)                                                                                                                gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                  <  1:160104/85‑1 (MQ=255)                                                                                                                gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                  >  2:160104/1‑85 (MQ=255)                                                                                                                gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAat                                                                        <  1:955066/79‑1 (MQ=255)                                                                                                                gtggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACAtctc            >  2:669487/1‑139 (MQ=255)                                                                                                                gtggtATCCCTGGCTATTAACAAGGCCAGGTTTTGATTCCATTCATTAAACATCCAGAAACAAAAGAACGTAAACAAat                                                                        >  2:955066/1‑79 (MQ=255)                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTcc                                                                                                               >  1:812089/1‑39 (MQ=255)                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTcc                                                                                                               <  2:812089/39‑1 (MQ=255)                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcat                                                                                                             <  1:905763/41‑1 (MQ=255)                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCcat                                                                                                             >  2:905763/1‑41 (MQ=255)                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatcc                                                                      >  1:485233/1‑80 (MQ=255)                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatcc                                                                      <  2:485233/80‑1 (MQ=255)                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                               <  1:1105671/87‑1 (MQ=255)                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                               >  2:1105671/1‑87 (MQ=255)                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATctct           >  2:1209100/1‑139 (MQ=255)                                                                                                                 tggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATctct           >  1:1168665/1‑139 (MQ=255)                                                                                                                  ggtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTa          >  1:72571/1‑139 (MQ=255)                                                                                                                   gtATCGCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                               >  1:706656/1‑85 (MQ=255)                                                                                                                   gtATCGCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCaga                                                               <  2:706656/85‑1 (MQ=255)                                                                                                                   gtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGAt                                                                                                  >  1:658008/1‑50 (MQ=255)                                                                                                                   gtATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGAt                                                                                                  <  2:658008/50‑1 (MQ=255)                                                                                                                    tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                   >  1:1095407/1‑80 (MQ=255)                                                                                                                    tATCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatc                                                                   <  2:1095407/80‑1 (MQ=255)                                                                                                                     aTCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAgg       >  2:981190/1‑139 (MQ=255)                                                                                                                     aTCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAgg       <  1:1209100/139‑1 (MQ=255)                                                                                                                      tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                  >  2:599164/1‑79 (MQ=255)                                                                                                                      tCCCTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                  <  1:599164/79‑1 (MQ=255)                                                                                                                       cccTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                <  2:446359/48‑1 (MQ=255)                                                                                                                       cccTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATcc                                                                                                >  1:446359/1‑48 (MQ=255)                                                                                                                        ccTGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATg    <  2:86759/139‑1 (MQ=255)                                                                                                                          tGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGaa  >  1:62279/1‑139 (MQ=255)                                                                                                                          tGGCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGaa  >  1:1075298/1‑139 (MQ=255)                                                                                                                           ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                  >  2:928420/1‑74 (MQ=255)                                                                                                                           ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAatccatcc                                                                  <  1:928420/74‑1 (MQ=255)                                                                                                                           ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGt                          >  2:844343/1‑114 (MQ=255)                                                                                                                           ggCTATTAACAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGt                          <  1:844343/114‑1 (MQ=255)                                                                                                                                   |                                                                                                                                  ATATCATTCCTGATTATACCTTTAACTTTATGGGAATTGGGGGTAGTGGTAAAAATGACCTGGATTATCCCTCTGCTGATATACGAGCGAAGATCAGTAAAGATATAGAAGGTGGTATCCCTGGCTATTAGCAAGGTCAGGTTTTGATTCCATTCATTAAAGATCCAGTAACAAAAGAACGTAAACAAATCCATCCAGAAGATATACATCTCACTGCAAAAGACTTCGAAGCAAGTAAAGATAACATCTCTAAGGATGAA  >  NC_000913/2787304‑2787563

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: ‑