Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	729,260	C→T	5.4%	L43L (CTG→TTG)	ybfA →	DUF2517 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	729,260	0	C	T	5.4%	240.8 / 2.2	74	L43L (CTG→TTG)	ybfA	DUF2517 family protein
Reads supporting (aligned to +/- strand):  ref base C (36/34);  new base T (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

ATACAAAAGGGGGAGAAAAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATG  >  NC_000913/729113‑729408                                                                                                                                                    |                                                                                                                                                     aTACAAAAGGGGGAGAAAAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCt                                                                                                                                                     <  1:22144/149‑1 (MQ=255)    cAAAAGGGGGAGAAAAGTATGGAACTCTACAGAGAATATCCTGCAAGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCa                                                                                                                                                  <  3:101335/149‑1 (MQ=255)           gggAGAAAAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTc                                                                                                                                           >  4:313432/1‑149 (MQ=255)               gaAAAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGt                                                                                                                                       >  4:122443/1‑149 (MQ=255)                aaaaGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTc                                                                                                                                      >  4:223489/1‑149 (MQ=255)                aaaaGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTc                                                                                                                                      >  1:125895/1‑149 (MQ=255)                 aaaGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCg                                                                                                                                     <  8:21567/149‑1 (MQ=255)                    gTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACTTGCATCGCGTCTGGTCGaa                                                                                                                                   >  3:311586/1‑148 (MQ=255)                    gTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACTTGCATCGCGTCTGGTCGaa                                                                                                                                   >  5:58858/1‑148 (MQ=255)                      aTGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAAc                                                                                                                                >  4:123500/1‑149 (MQ=255)                          aaCTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGc                                                                                                                            <  4:255409/149‑1 (MQ=255)                               tACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATa                                                                                                                        >  5:205975/1‑148 (MQ=255)                                 cAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAAc                                                                                                                     >  4:280086/1‑149 (MQ=255)                                  agagAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAAcc                                                                                                                    >  6:15754/1‑149 (MQ=255)                                    agaATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCTTGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCgg                                                                                                                  >  8:318435/1‑149 (MQ=255)                                     gaAAATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGgt                                                                                                                 <  1:23714/149‑1 (MQ=255)                                        tatCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGtgtg                                                                                                              >  8:151376/1‑149 (MQ=255)                                        tatCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCATTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGtgtg                                                                                                              >  5:147922/1‑149 (MQ=255)                                              gCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGtggatgga                                                                                                        >  5:314093/1‑149 (MQ=255)                                               cctGGCTTATCTTTTTACTCCGTACTTATGCGGGTGCAGCGGGCGTTCTGTCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATAGCGTCTGGGCGAAAACCAGCGATAAACCGGTGTggatggat                                                                                                       <  1:225199/147‑1 (MQ=255)                                               cATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTggatggat                                                                                                       >  3:101318/1‑149 (MQ=255)                                                    cTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGc                                                                                                  <  2:108780/149‑1 (MQ=255)                                                     ttATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGcc                                                                                                 >  4:219043/1‑149 (MQ=255)                                                      tATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGcc                                                                                                 >  3:175396/1‑148 (MQ=255)                                                      tATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGcc                                                                                                 >  2:307566/1‑148 (MQ=255)                                                                  ccGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTg                                                                                    <  3:123500/149‑1 (MQ=255)                                                                    gTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGt                                                                                  >  6:19550/1‑149 (MQ=255)                                                                    gTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGAGGAAAACCAGCATAAAACCGGTGGGGATGGAGCAGGCCGAAAAGCAAAATGGt                                                                                  >  6:8063/1‑149 (MQ=255)                                                                          aTGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACTTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATTGATCAGGCCGAAAAGGCAACTGGTGAttt                                                                             <  6:58858/148‑1 (MQ=255)                                                                          aTGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACTTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGAttt                                                                             <  4:311586/148‑1 (MQ=255)                                                                            gCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTa                                                                          >  4:204794/1‑149 (MQ=255)                                                                             cGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTa                                                                          <  6:314093/148‑1 (MQ=255)                                                                                ttGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTattga                                                                      <  3:189555/149‑1 (MQ=255)                                                                                ttGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTattga                                                                      <  3:185202/149‑1 (MQ=255)                                                                                   cAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTCCCTGCATCGCGTCTGGTCGAAAACCAGCGTTAAACCGGTGTGGATGGAGCAGGCCGAAAAGTCAACTGGTGATTTTTattgatt                                                                    <  1:61963/148‑1 (MQ=255)                                                                                   cAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTattgattg                                                                   >  7:266610/1‑149 (MQ=255)                                                                                    aGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGc                                                                  >  4:159718/1‑149 (MQ=255)                                                                                    aGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGc                                                                  >  2:257839/1‑149 (MQ=255)                                                                                     gCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTCCCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAACGGCAACTGGTGATTTTTATTGATTGCt                                                                 >  5:74888/1‑149 (MQ=255)                                                                                       gggCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCt                                                                 >  8:165372/1‑147 (MQ=255)                                                                                       gggCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGa                                                               <  6:101918/149‑1 (MQ=255)                                                                                            ttCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTTGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGtttt                                                          <  7:318435/149‑1 (MQ=255)                                                                                            ttCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGtttt                                                          <  3:106287/149‑1 (MQ=255)                                                                                                   gcTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACa                                                   >  2:130133/1‑149 (MQ=255)                                                                                                       ccTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACAtt                                                 <  8:270158/147‑1 (MQ=255)                                                                                                       ccTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACAtt                                                 <  5:280424/147‑1 (MQ=255)                                                                                                       ccTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTg                                                >  5:24950/1‑148 (MQ=255)                                                                                                       ccTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGa                                               <  5:15754/149‑1 (MQ=255)                                                                                                           tCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGaaaa                                            <  3:122443/148‑1 (MQ=255)                                                                                                            cATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAAcc                                          >  8:249559/1‑149 (MQ=255)                                                                                                             aTGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAAcc                                          <  1:307566/148‑1 (MQ=255)                                                                                                             aTGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCg                                         >  7:217602/1‑149 (MQ=255)                                                                                                               gCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGc                                        >  4:54655/1‑148 (MQ=255)                                                                                                                cTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGc                                        >  4:60064/1‑147 (MQ=255)                                                                                                                   ttCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAg                                     <  3:204794/147‑1 (MQ=255)                                                                                                                      tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTAGAAAACCAGCGAGAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCAAGCGGa                                 >  1:266490/1‑148 (MQ=255)                                                                                                                         aaaGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACg                             <  8:217602/149‑1 (MQ=255)                                                                                                                             aCCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTgg                         <  3:280086/149‑1 (MQ=255)                                                                                                                              ccGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGc                        <  3:303243/149‑1 (MQ=255)                                                                                                                               cgcgCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAAGCGGTGTGGATGGATCAGGCCGAAAAGGCAAGTGGTGATTTTTCTTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCg                       <  3:313432/149‑1 (MQ=255)                                                                                                                                    ccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGa                                                               <  3:321934/104‑1 (MQ=255)                                                                                                                                    ccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGa                                                               >  4:321934/1‑104 (MQ=255)                                                                                                                                     cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGttttt                 <  5:146134/149‑1 (MQ=255)                                                                                                                                     cGCTTATACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAACGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATAGATTGAAAACCGCCAGCGGAAACGCGGGCGGttttt                 <  1:333324/149‑1 (MQ=255)                                                                                                                                      gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTg                                                         >  3:196520/1‑108 (MQ=255)                                                                                                                                      gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTg                                                         <  4:196520/108‑1 (MQ=255)                                                                                                                                      gCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                <  3:219043/149‑1 (MQ=255)                                                                                                                                       cTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTg               >  7:326353/1‑149 (MQ=255)                                                                                                                                        ttttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTg               >  8:200518/1‑148 (MQ=255)                                                                                                                                          ttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGt              <  3:223489/147‑1 (MQ=255)                                                                                                                                           tACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTAt            >  8:212304/1‑148 (MQ=255)                                                                                                                                               gCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGa       <  4:175396/149‑1 (MQ=255)                                                                                                                                                cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGaa      <  7:249559/149‑1 (MQ=255)                                                                                                                                                cTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGaa      <  1:257839/149‑1 (MQ=255)                                                                                                                                                   ccTGCATCGCGTCTGGTCGAAAACCCTCGATAAACCGGTGGGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGAt   <  4:101318/149‑1 (MQ=255)                                                                                                                                                   ccTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGa    <  7:151376/148‑1 (MQ=255)                                                                                                                                                    cTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATg  >  6:252022/1‑149 (MQ=255)                                                                                                                                                    cTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATg  >  3:219550/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ATACAAAAGGGGGAGAAAAGTATGGAACTCTACAGAGAATATCCTGCATGGCTTATCTTTTTACGCCGTACTTATGCGGTTGCAGCGGGCGTTCTGGCGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATG  >  NC_000913/729113‑729408

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