Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,500,043	C→T	pseudogene (492/942 nt)	yjgX ←	pseudogene fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,500,043	0	C	T	100.0%	253.6 / NA	79	pseudogene (492/942 nt)	yjgX	pseudogene fragment
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/43);  total (36/43)

AAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAG  >  NC_000913/4499907‑4500176                                                                                                                                         |                                                                                                                                      aaaaGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGTTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTa                                                                                                                                     >  2:917904/1‑139 (MQ=255)      tAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGt                                                                                                                                >  1:1271748/1‑139 (MQ=255)               gCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGt                                                                                                                       >  1:1374603/1‑139 (MQ=255)                   aaTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                   <  1:830187/139‑1 (MQ=255)                   aaTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                   <  2:494010/139‑1 (MQ=255)                   aaTGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                   <  2:1449341/139‑1 (MQ=255)                        cGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTa                                                                                                                                     >  2:800941/1‑116 (MQ=255)                        cGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTa                                                                                                                                     <  1:800941/116‑1 (MQ=255)                        cGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGAt                                                                                                              >  2:426379/1‑139 (MQ=255)                            aGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTgg                                                                                                          >  1:1013069/1‑139 (MQ=255)                              tACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGAt                                                                                                        <  2:453283/139‑1 (MQ=255)                               aCCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATg                                                                                                       <  2:1013069/139‑1 (MQ=255)                                 cGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                                     <  1:577922/139‑1 (MQ=255)                                 cGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                                     <  1:197873/139‑1 (MQ=255)                                   gCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGaa                                                                                                   >  1:483358/1‑139 (MQ=255)                                              cATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTaa                                                                                                                             >  2:660589/1‑102 (MQ=255)                                              cATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTaa                                                                                                                             <  1:660589/102‑1 (MQ=255)                                                tCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTacac                                                                                                                  >  2:1472825/1‑111 (MQ=255)                                                tCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTacac                                                                                                                  <  1:1472825/111‑1 (MQ=255)                                                        gCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGc                                                                                                                          >  2:112971/1‑95 (MQ=255)                                                        gCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGc                                                                                                                          <  1:112971/95‑1 (MQ=255)                                                        gCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAAtt                                                                              >  2:1641328/1‑139 (MQ=255)                                                        gCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAAtt                                                                              >  2:340284/1‑139 (MQ=255)                                                             gAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCtt                                                                         >  1:1309537/1‑139 (MQ=255)                                                               aaGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCg                                                                       <  1:1086129/139‑1 (MQ=255)                                                                  ttATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTacac                                                                                                                  >  2:968698/1‑93 (MQ=255)                                                                  ttATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTacac                                                                                                                  <  1:968698/93‑1 (MQ=255)                                                                  ttATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGgag                                                                    <  2:129305/139‑1 (MQ=255)                                                                  ttATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGgag                                                                    <  2:80274/139‑1 (MQ=255)                                                                     tCCAGTTTACTGATGTTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                 <  1:917904/139‑1 (MQ=255)                                                                     tCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAgg                                                                 <  2:1374603/139‑1 (MQ=255)                                                                       cAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGAt                                                               <  2:1191165/139‑1 (MQ=255)                                                                          tttACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTa                                                            >  2:1245188/1‑139 (MQ=255)                                                                            tACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATc                                                          >  2:758222/1‑139 (MQ=255)                                                                                  tGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAg                                                    <  2:1533429/139‑1 (MQ=255)                                                                                   gCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgc                                                   >  1:766545/1‑139 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                   <  2:909266/74‑1 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                   >  1:909266/1‑74 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                   <  2:218200/74‑1 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTtaca                                                                                                                   >  1:218200/1‑74 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGAttt                                                                                                            >  2:1350708/1‑81 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGAttt                                                                                                            <  1:1350708/81‑1 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAAt                                                                               <  2:1574163/110‑1 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAAt                                                                               >  1:1574163/1‑110 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcg                                                  >  2:1614589/1‑139 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcg                                                  >  2:606849/1‑139 (MQ=255)                                                                                    cTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAACTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcg                                                  >  2:1561818/1‑139 (MQ=255)                                                                                     tGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                 <  1:1614589/139‑1 (MQ=255)                                                                                      ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAgg                                                                                                                           >  1:599328/1‑64 (MQ=255)                                                                                      ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAgg                                                                                                                           <  2:599328/64‑1 (MQ=255)                                                                                      ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCg                                                        <  2:138767/131‑1 (MQ=255)                                                                                      ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCg                                                        >  1:138767/1‑131 (MQ=255)                                                                                      ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGcc                                                      <  1:1546513/133‑1 (MQ=255)                                                                                      ggCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGcc                                                      >  2:1546513/1‑133 (MQ=255)                                                                                          aGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGtt                                            <  1:1245188/139‑1 (MQ=255)                                                                                          aGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGtt                                            <  2:77343/139‑1 (MQ=255)                                                                                          aGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGtt                                            <  1:1641328/139‑1 (MQ=255)                                                                                            gTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGAt                                                               >  2:1549224/1‑118 (MQ=255)                                                                                            gTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGAt                                                               <  1:1549224/118‑1 (MQ=255)                                                                                                  ggCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTcc                                    <  1:606849/139‑1 (MQ=255)                                                                                                     aaTTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                                     >  2:129058/1‑71 (MQ=255)                                                                                                     aaTTTTTCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATgag                                                                                                     <  1:129058/71‑1 (MQ=255)                                                                                                       tttttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGtt                                                                                                                               >  2:1428586/1‑43 (MQ=255)                                                                                                       tttttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGtt                                                                                                                               <  1:1428586/43‑1 (MQ=255)                                                                                                       tttttCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTaa                               <  1:340284/139‑1 (MQ=255)                                                                                                           tCGGGGATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTaacaac                           <  2:483358/139‑1 (MQ=255)                                                                                                                gATACTGCTGCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATTAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGttt                      <  2:979736/139‑1 (MQ=255)                                                                                                                gATACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGttt                      <  1:758222/139‑1 (MQ=255)                                                                                                                  tACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGaa                                                                                                   <  1:847727/60‑1 (MQ=255)                                                                                                                  tACTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGaa                                                                                                   >  2:847727/1‑60 (MQ=255)                                                                                                                   aCTGCTCCAGGTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTc                   <  2:1309537/139‑1 (MQ=255)                                                                                                                            ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGt                                        >  2:334002/1‑109 (MQ=255)                                                                                                                            ggTAATTATTCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGt                                        <  1:334002/109‑1 (MQ=255)                                                                                                                                 ttattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                 >  1:145555/1‑95 (MQ=255)                                                                                                                                 ttattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAgcgc                                                 <  2:145555/95‑1 (MQ=255)                                                                                                                                  tattCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGga    >  1:479478/1‑139 (MQ=255)                                                                                                                                   attCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaa   >  1:601256/1‑139 (MQ=255)                                                                                                                                    ttCGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGgaag  <  1:220573/139‑1 (MQ=255)                                                                                                                                      cGGTTAGGAGTTAAGGCTGTTACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCtt                                                                         >  1:238065/1‑66 (MQ=255)                                                                                                                                         |                                                                                                                                      AAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTAATTATTCGGCTAGGAGTTAAGGCTGTCACACGGATTTGGATGAGAACCCATCATGTGCAGGAAAATTATCTTCGGAGAGGATGTATCCGCCAGCGCACGTTCTGTTTCCTGTAACAACAGTTTGTCATCTGTTTTACGGGAAG  >  NC_000913/4499907‑4500176

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