Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	658,577	G→A	5.5%	pseudogene (353/594 nt)	ybeM →	pseudogene, C‑N hydrolase superfamily homology;putative enzyme; Not classified; putative amidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	658,577	0	G	A	5.5%	244.6 / 3.0	73	pseudogene (353/594 nt)	ybeM	pseudogene, C‑N hydrolase superfamily homology;putative enzyme; Not classified; putative amidase
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base A (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

GAGGTGGAAGGGATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTC  >  NC_000913/658435‑658722                                                                                                                                               |                                                                                                                                                  gAGGTGGAAGGGATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGCGCATCACTGGTCAAc                                                                                                                                             <  4:543498/149‑1 (MQ=255) gAGGTGGAAGGGATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAAc                                                                                                                                             >  2:344319/1‑149 (MQ=255)      ggAAGGGATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGttgc                                                                                                                                        <  3:163526/149‑1 (MQ=255)      ggAAGGGATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGttgc                                                                                                                                        <  3:398377/149‑1 (MQ=255)       gAAGGGATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGttgct                                                                                                                                       >  3:210318/1‑149 (MQ=255)       gAAGGGATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGttgct                                                                                                                                       >  4:408362/1‑149 (MQ=255)        aaGGGATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGttgctt                                                                                                                                      <  4:2144/149‑1 (MQ=255)         aGGGATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGttgcttg                                                                                                                                     >  1:143571/1‑149 (MQ=255)            gATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCg                                                                                                                                  >  4:205103/1‑149 (MQ=255)             aTGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGc                                                                                                                                 >  3:219215/1‑149 (MQ=255)             aTGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTAGTCAACGTTGCTTGCCGc                                                                                                                                 >  3:169547/1‑149 (MQ=255)             aTGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTAGTCAACGTTGCTTGCCGc                                                                                                                                 >  1:513214/1‑149 (MQ=255)              tGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCt                                                                                                                                >  1:33530/1‑149 (MQ=255)               gAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTc                                                                                                                               >  4:302711/1‑149 (MQ=255)                aaGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCg                                                                                                                              <  3:27825/149‑1 (MQ=255)                  ggtcggtcTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTg                                                                                                                            <  2:106472/149‑1 (MQ=255)                  ggtcggtcTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTg                                                                                                                            >  1:299549/1‑149 (MQ=255)                  ggtcggtcTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTg                                                                                                                            <  2:89435/149‑1 (MQ=255)                      ggtcTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCt                                                                                                                        >  4:342029/1‑149 (MQ=255)                       gtcTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTg                                                                                                                       <  2:156014/149‑1 (MQ=255)                          tgatgaCCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGAt                                                                                                                    >  3:73068/1‑149 (MQ=255)                           gatgaCCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATa                                                                                                                   >  3:14227/1‑149 (MQ=255)                              gaCCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATacca                                                                                                                >  2:368358/1‑149 (MQ=255)                                ccTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATaccacc                                                                                                              <  1:362470/149‑1 (MQ=255)                                  tGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTg                                                                                                            >  1:571254/1‑149 (MQ=255)                                   gTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGt                                                                                                           <  4:526588/149‑1 (MQ=255)                                   gTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCGGt                                                                                                           <  4:352263/149‑1 (MQ=255)                                    ttATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGtt                                                                                                          <  1:77477/149‑1 (MQ=255)                                       tGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTtata                                                                                                       <  3:408362/149‑1 (MQ=255)                                        gACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTatat                                                                                                      <  3:302711/149‑1 (MQ=255)                                         aCTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATAtg                                                                                                     >  4:549202/1‑149 (MQ=255)                                         aCTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATAtg                                                                                                     >  4:394053/1‑149 (MQ=255)                                         aCTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATAtg                                                                                                     >  4:396516/1‑149 (MQ=255)                                           ttACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATAtggt                                                                                                   <  1:423882/149‑1 (MQ=255)                                           ttACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATAtggt                                                                                                   >  4:224488/1‑149 (MQ=255)                                             aCGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATAtggtgg                                                                                                 <  2:143571/149‑1 (MQ=255)                                                     cAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCgggg                                                                                         >  1:520897/1‑149 (MQ=255)                                                     cAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCgggg                                                                                         >  4:431100/1‑149 (MQ=255)                                                      agagCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCggggg                                                                                        >  4:397008/1‑149 (MQ=255)                                                      agagCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCggggg                                                                                        >  1:439348/1‑149 (MQ=255)                                                                ctggcACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGaa                                                                              <  2:173744/149‑1 (MQ=255)                                                                ctggcACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGaa                                                                              <  4:14227/149‑1 (MQ=255)                                                                   gcACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACaa                                                                           >  4:116918/1‑149 (MQ=255)                                                                      caGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACaaaaa                                                                        <  2:43716/149‑1 (MQ=255)                                                                         gCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAatat                                                                     <  4:395269/149‑1 (MQ=255)                                                                         gCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAatat                                                                     <  2:221182/149‑1 (MQ=255)                                                                             tACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGt                                                                 <  3:493156/149‑1 (MQ=255)                                                                                 gggAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCaaa                                                             <  4:219215/149‑1 (MQ=255)                                                                                    aGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGcc                                                          <  4:73068/149‑1 (MQ=255)                                                                                      cTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCgg                                                        >  3:368658/1‑149 (MQ=255)                                                                                        gAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTAGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGAt                                                      <  4:169547/149‑1 (MQ=255)                                                                                        gAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTAGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGAt                                                      <  2:513214/149‑1 (MQ=255)                                                                                         aaaTCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGAtt                                                     <  3:205103/149‑1 (MQ=255)                                                                                            tCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATtata                                                  <  3:549202/149‑1 (MQ=255)                                                                                               tGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGAt                                               >  4:136683/1‑149 (MQ=255)                                                                                               tGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGAt                                               >  3:442557/1‑149 (MQ=255)                                                                                               tGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGAt                                               >  1:5995/1‑149 (MQ=255)                                                                                                     ttCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCttt                                         <  4:210318/149‑1 (MQ=255)                                                                                                       ccTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTgg                                       <  2:520897/149‑1 (MQ=255)                                                                                                        cTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGc                                      >  3:334625/1‑149 (MQ=255)                                                                                                                 gggTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCAtt                             <  4:362611/149‑1 (MQ=255)                                                                                                                    ttCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCg                          >  3:211426/1‑149 (MQ=255)                                                                                                                                tCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGaa              >  2:477228/1‑149 (MQ=255)                                                                                                                                tCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGaa              >  3:146411/1‑149 (MQ=255)                                                                                                                                  aaaGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAAt            >  3:547063/1‑149 (MQ=255)                                                                                                                                  aaaGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAAt            >  3:83583/1‑149 (MQ=255)                                                                                                                                   aaGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGcggcgg                                                                                           >  3:552881/1‑69 (MQ=255)                                                                                                                                   aaGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGcggcgg                                                                                           <  4:552881/69‑1 (MQ=255)                                                                                                                                     gagCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGcc         <  2:170714/149‑1 (MQ=255)                                                                                                                                         aTCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCa     >  4:127981/1‑149 (MQ=255)                                                                                                                                           cACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACt   >  2:521966/1‑149 (MQ=255)                                                                                                                                            aCTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTc  >  1:300079/1‑149 (MQ=255)                                                                                                                                            aCTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTc  >  3:173028/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                  GAGGTGGAAGGGATGAAGGTCGGTCTGATGACCTGTTATGACTTACGCTTTCCAGAGCTGGCGCTGGCACAGGCATTACAGGGAGCTGAAATCCTGGTACTTCCTGCCGCCTGGGTTCGCGGGCCGCTCAAAGAGCATCACTGGTCAACGTTGCTTGCCGCTCGTGCGCTGGATACCACCTGTTATATGGTGGCGGCGGGGGAGTGCGGGAACAAAAATATCGGTCAAAGCCGGATTATAGATCCCTTTGGCGTCACCATTGCGGCAGCGTCAGAAATGCCTGCACTC  >  NC_000913/658435‑658722

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