Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,994,415	G→A	pseudogene (34/60 nt)	pbl →	pseudogene, peptidoglycan‑binding enzyme family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,994,415	0	G	A	100.0%	248.2 / NA	77	pseudogene (34/60 nt)	pbl	pseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (45/32);  total (45/32)

TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGA  >  NC_000913/2994287‑2994550                                                                                                                                 |                                                                                                                                         ttttCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTaa                                                                                                                                >  2:484970/1‑139 (MQ=255)        gTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCa                                                                                                                         >  1:563145/1‑139 (MQ=255)         ttGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGCGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAg                                                                                                                        <  2:486452/139‑1 (MQ=255)               aGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAAt                                                                                                                  <  2:563145/139‑1 (MQ=255)               aGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAAt                                                                                                                  <  1:484970/139‑1 (MQ=255)                    cGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGcc                                                                                                             >  2:384576/1‑139 (MQ=255)                              ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc‑t                                                                                                  <  1:384576/139‑1 (MQ=255)                              ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc‑t                                                                                                  <  2:528819/139‑1 (MQ=255)                              ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc‑t                                                                                                  >  1:528819/1‑139 (MQ=255)                              ccTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc‑t                                                                                                  <  1:722049/139‑1 (MQ=255)                                        tataATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTtata                                                                                        <  1:447455/139‑1 (MQ=255)                                                 ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaaa                                                                               <  1:672062/139‑1 (MQ=255)                                                 ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaaa                                                                               <  2:381655/139‑1 (MQ=255)                                                 ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaaa                                                                               <  2:212606/139‑1 (MQ=255)                                                 ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaaa                                                                               <  2:116229/139‑1 (MQ=255)                                                 ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaaa                                                                               >  2:124690/1‑139 (MQ=255)                                                 ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaaa                                                                               >  1:116229/1‑139 (MQ=255)                                                 ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaaa                                                                               >  1:381655/1‑139 (MQ=255)                                                 ggTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaaa                                                                               >  2:672062/1‑139 (MQ=255)                                                  gTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTTATATAGTGCaaa                                                                               <  1:132879/139‑1 (MQ=255)                                                   ttttGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAaca                                                                             >  2:41861/1‑139 (MQ=255)                                                                aaCCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAg                                                                >  2:259786/1‑139 (MQ=255)                                                                          agagACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGAt                                                      <  1:124690/139‑1 (MQ=255)                                                                           gagaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                                     <  1:614526/139‑1 (MQ=255)                                                                            agaCAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCa                                                    >  2:572859/1‑139 (MQ=255)                                                                               cAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                                 <  1:572859/139‑1 (MQ=255)                                                                                            cAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACAccc                                                                          <  1:250200/101‑1 (MQ=255)                                                                                            cAAAATATATATTGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACAccc                                                                          >  2:250200/1‑101 (MQ=255)                                                                                                        tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc‑t                                                                                                  >  1:215692/1‑65 (MQ=255)                                                                                                        tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc‑t                                                                                                  <  2:215692/65‑1 (MQ=255)                                                                                                        tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCct                                                                                                     <  1:168052/63‑1 (MQ=255)                                                                                                        tGTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCct                                                                                                     >  2:168052/1‑63 (MQ=255)                                                                                                         gTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                            <  1:661898/118‑1 (MQ=255)                                                                                                         gTATATACCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATa                                            >  2:661898/1‑118 (MQ=255)                                                                                                           atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTtat                                                                                         <  1:132919/71‑1 (MQ=255)                                                                                                           atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTtat                                                                                         >  2:132919/1‑71 (MQ=255)                                                                                                           atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                                   >  1:738584/1‑109 (MQ=255)                                                                                                           atataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCaa                                                   <  2:738584/109‑1 (MQ=255)                                                                                                             ataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc‑t                                                                                                  <  2:188676/60‑1 (MQ=255)                                                                                                             ataCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc‑t                                                                                                  >  1:188676/1‑60 (MQ=255)                                                                                                               acgaGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCctc‑t                                                                                                  <  1:245192/55‑1 (MQ=255)                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATg                 >  1:94774/1‑139 (MQ=255)                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATg                 >  1:399332/1‑139 (MQ=255)                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATg                 >  1:140539/1‑139 (MQ=255)                                                                                                               aCCAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGActc‑t                                                                                                  >  2:245192/1‑58 (MQ=255)                                                                                                                 cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTTTATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGat               >  2:502217/1‑139 (MQ=255)                                                                                                                 cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaa                                                                                >  1:760182/1‑74 (MQ=255)                                                                                                                 cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaa                                                                                <  2:760182/74‑1 (MQ=255)                                                                                                                 cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGat               >  1:250481/1‑139 (MQ=255)                                                                                                                 cAGGCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGat               >  2:668245/1‑139 (MQ=255)                                                                                                                    tcTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaaa                                                                               <  1:640231/71‑1 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATg                                                                                                                 <  2:150761/39‑1 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATg                                                                                                                 >  1:150761/1‑39 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGc                                                                                                       >  1:193234/1‑49 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGc                                                                                                       <  2:193234/49‑1 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGc                                                                                  >  1:77952/1‑69 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGc                                                                                  <  2:77952/69‑1 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCaaa                                                                               >  2:640231/1‑72 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                                 <  2:7088/102‑1 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                                 <  2:346058/102‑1 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                                 >  1:7088/1‑102 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACt                                                 >  1:346058/1‑102 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAAt                          <  2:86055/125‑1 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAAt                          >  1:86055/1‑125 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            >  2:607151/1‑139 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            >  2:450735/1‑139 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            >  1:661376/1‑139 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            >  2:674310/1‑139 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            >  1:70514/1‑139 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            >  1:114378/1‑139 (MQ=255)                                                                                                                    gCTTAATGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTGCATGATAtt            >  2:522539/1‑139 (MQ=255)                                                                                                                        aaTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGt        >  2:499716/1‑139 (MQ=255)                                                                                                                         aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGAt                                                      >  1:377190/1‑92 (MQ=255)                                                                                                                         aTGAACTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGAt                                                      <  2:377190/92‑1 (MQ=255)                                                                                                                              cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATc                                                     >  1:530925/1‑88 (MQ=255)                                                                                                                              cTAAACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGa  >  1:68738/1‑139 (MQ=255)                                                                                                                              cTAAACAATCGTAAGGCTACAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATAtt            >  2:428172/1‑129 (MQ=255)                                                                                                                                 |                                                                                                                                         TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATATATATTGTATATACCAGGCTTAATGAACTAGACAATCGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTC‑TTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGA  >  NC_000913/2994287‑2994550

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