BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	512,212	T→G	5.4%	V191G (GTG→GGG)	glsA →	glutaminase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	512,212	0	T	G	5.4%	424.2 / 10.6	144	V191G (GTG→GGG)	glsA	glutaminase 1
Reads supporting (aligned to +/- strand): ref base T (67/69); new base G (3/5); total (70/74)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 151 total aligned reads displayed.

TATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGG  >  NC_000913/512063‑512361

tatCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgt                                                                                                                                                       >  1:817028/1‑150 (MQ=255)

tatCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgt                                                                                                                                                       >  1:1803434/1‑150 (MQ=255)

 atCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgtg                                                                                                                                                      <  2:2359417/150‑1 (MQ=255)

 atCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACgtg                                                                                                                                                      >  2:590038/1‑150 (MQ=255)

   ccAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGta                                                                                                                                                    <  2:705285/150‑1 (MQ=255)

    cAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGtat                                                                                                                                                   <  1:2165970/150‑1 (MQ=255)

    cAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGtat                                                                                                                                                   <  2:327854/150‑1 (MQ=255)

         gCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCg                                                                                                                                              <  1:657342/150‑1 (MQ=255)

         gCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCg                                                                                                                                              <  2:1201724/150‑1 (MQ=255)

         gCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCg                                                                                                                                              <  1:864091/150‑1 (MQ=255)

           aaCTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTc                                                                                                                                            >  1:184343/1‑150 (MQ=255)

           aaCTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTc                                                                                                                                            <  1:1118409/150‑1 (MQ=255)

              tggctggcGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGt                                                                                                                                         <  1:1729293/150‑1 (MQ=255)

              tggctggcGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGt                                                                                                                                         <  1:1402420/150‑1 (MQ=255)

              tggctggcGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGt                                                                                                                                         <  1:290698/150‑1 (MQ=255)

               ggctggcGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGGGTATACACGTCacag                                                                                                                                        <  2:6886/150‑4 (MQ=255)

                gctggcGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGc                                                                                                                                       <  1:1354351/150‑1 (MQ=255)

                 ctggcGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCt                                                                                                                                      >  1:1373631/1‑150 (MQ=255)

                 ctggcGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGGGTATACACGTCacaggc                                                                                                                                      <  2:1417600/150‑6 (MQ=255)

                     cGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCAc                                                                                                                                  >  1:712329/1‑150 (MQ=255)

                       aGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGc                                                                                                                                <  1:968821/150‑1 (MQ=255)

                          aGGTAGCGCTCTCTGACGAAGTCAAACAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATCCAATGGAAGCCTGTCACGTGTATACCCGTCAGTGCTCCACGctcc                                                                                                                             >  2:1605703/1‑150 (MQ=255)

                           ggTACCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAAACAACTTACATAACCGGGACATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctcct                                                                                                                            >  2:77213/1‑150 (MQ=255)

                            gTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGctcctc                                                                                                                           <  1:2355469/150‑1 (MQ=255)

                             tAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCa                                                                                                                          <  2:1840680/150‑1 (MQ=255)

                                   tctctGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTa                                                                                                                    >  2:915204/1‑150 (MQ=255)

                                    ctctGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTAt                                                                                                                   >  1:2254458/1‑150 (MQ=255)

                                      ctGACGATGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTg                                                                                                                 >  2:1822030/1‑150 (MQ=255)

                                       tGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGa                                                                                                                >  1:1097553/1‑150 (MQ=255)

                                         aCGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAAc                                                                                                              >  1:1910466/1‑150 (MQ=255)

                                           taaGTCGACCAGTCGTAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTg                                                                                                            <  1:2070802/149‑1 (MQ=255)

                                             aGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGc                                                                                                          >  1:649467/1‑150 (MQ=255)

                                              gTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCa                                                                                                         >  1:2077253/1‑150 (MQ=255)

                                                     aGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTg                                                                                                  >  2:220595/1‑150 (MQ=255)

                                                         ggAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcgc                                                                                              <  2:1280224/150‑1 (MQ=255)

                                                         ggAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcgc                                                                                              >  1:228374/1‑150 (MQ=255)

                                                         ggAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcgc                                                                                              >  1:856238/1‑150 (MQ=255)

                                                         ggAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGgcgc                                                                                              >  1:1101975/1‑150 (MQ=255)

                                                              aaacaaCCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGc                                                                                         >  2:215531/1‑150 (MQ=255)

                                                              aaacaaCCAACTTCCATAACAGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAAATGGCAACGCTTGGCGcgccgc                                                                                         >  2:830118/1‑150 (MQ=255)

                                                                   accaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcg                                                                                    <  1:1824559/150‑1 (MQ=255)

                                                                   accaacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcg                                                                                    >  1:1024204/1‑150 (MQ=255)

                                                                   accaacTTCCATAACCGGGCCATAGACTGGCTGCTGTACTCCGCCGGATATCTCTATTGTTATGCAATGTAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggct                                                                                    <  1:1523185/150‑2 (MQ=255)

                                                                     caacTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTggcggc                                                                                  >  2:952996/1‑150 (MQ=255)

                                                                         ttCCATAACCGGGACATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggt                                                                              >  2:524934/1‑150 (MQ=255)

                                                                           ccATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtgg                                                                            >  1:1731940/1‑150 (MQ=255)

                                                                            cATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAggtggt                                                                           >  1:2267766/1‑150 (MQ=255)

                                                                             aTAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGgtggtg                                                                          <  1:1708457/150‑1 (MQ=255)

                                                                                 ccGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGaa                                                                      >  2:2349119/1‑150 (MQ=255)

                                                                                       cATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGtt                                                                >  2:1207066/1‑150 (MQ=255)

                                                                                        aTAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTg                                                               <  1:2004970/150‑1 (MQ=255)

                                                                                            ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGc                                                           <  2:2148942/150‑1 (MQ=255)

                                                                                            ccTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTCGCAACGCTTGGCGCGACGCTGGCGGCAGTTGGTGTGAATCCGTTGACGc                                                           >  2:1292224/1‑150 (MQ=255)

                                                                                              tGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCAt                                                         >  1:335970/1‑150 (MQ=255)

                                                                                               ggCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCCGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTCTGAATCCGTTGACGCATa                                                        >  1:1426455/1‑150 (MQ=255)

                                                                                               ggCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATa                                                        <  2:1962830/150‑1 (MQ=255)

                                                                                               ggCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATa                                                        <  2:2267766/150‑1 (MQ=255)

                                                                                                     tgTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcgcg                                                  <  1:986190/150‑1 (MQ=255)

                                                                                                     tgTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAAcgcg                                                  <  2:1101975/150‑1 (MQ=255)

                                                                                                      ttACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTTACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTTGCAACTCTTTGCGCGACGCTGGCGGAAGGTTGTGTGAATCCGTTGTCGCATAAAGGCGt                                                 <  2:634660/149‑1 (MQ=255)

                                                                                                      gTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGt                                                 <  1:747259/150‑1 (MQ=255)

                                                                                                        aCTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACAGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGttc                                               <  1:480086/150‑1 (MQ=255)

                                                                                                         cTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGttct                                              >  1:1311939/1‑150 (MQ=255)

                                                                                                           ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGttcttc                                            >  2:1948586/1‑150 (MQ=255)

                                                                                                           ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGttcttc                                            <  2:564009/150‑1 (MQ=255)

                                                                                                           ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGttcttc                                            >  2:1845370/1‑150 (MQ=255)

                                                                                                           ccgccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGATTGGCGCGACGCTGGCGGCAGGTGTTGTGAATCCGTTGACGCATAAACGCCttcttc                                            >  1:2329167/1‑150 (MQ=255)

                                                                                                           acgcagGATATCTCTATTGAGATGAAATGGAATCCTGTGACGTGTATACAAGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGttcttc                                            >  2:1958317/2‑150 (MQ=255)

                                                                                                             gccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTACACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAg                                          >  1:516464/1‑150 (MQ=255)

                                                                                                             gccgGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGGGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCGGGCGGCAGGTGGGGTGAATCCGTTGACGCATAAACGCGTTCTTCAg                                          <  2:759002/150‑1 (MQ=255)

                                                                                                                ggATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGcc                                       <  1:2342067/150‑1 (MQ=255)

                                                                                                                 gATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCg                                      <  1:1416556/150‑1 (MQ=255)

                                                                                                                   tatCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGAc                                    >  1:1073696/1‑150 (MQ=255)

                                                                                                                     tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACaa                                  <  2:479480/150‑1 (MQ=255)

                                                                                                                     tctcTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACaa                                  >  1:2190694/1‑150 (MQ=255)

                                                                                                                      ctctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAAc                                 <  2:622480/150‑1 (MQ=255)

                                                                                                                      ctctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAAc                                 >  2:686085/1‑150 (MQ=255)

                                                                                                                      ctctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATAGAACTGCCAACGCTTGCCGCGACGCTGGCGGCAGCTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCCACAAc                                 >  2:1935075/1‑150 (MQ=255)

                                                                                                                       gctATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACg                                <  1:668371/149‑1 (MQ=255)

                                                                                                                           ttGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGcc                            >  1:1872217/1‑150 (MQ=255)

                                                                                                                              tgATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTa                         <  1:2107994/150‑1 (MQ=255)

                                                                                                                                aTGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACa                       >  1:1167867/1‑150 (MQ=255)

                                                                                                                                 tGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACAt                      >  1:90326/1‑150 (MQ=255)

                                                                                                                                 tGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACAt                      >  1:34155/1‑150 (MQ=255)

                                                                                                                                   cAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTc                    >  1:2391268/1‑150 (MQ=255)

                                                                                                                                    aaTGGAAGCCTGTGACGTGTATACCCGTCAGTGGTCCACGCTCCTCAATACTATTGATCTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                   <  1:1991854/150‑1 (MQ=255)

                                                                                                                                    aaTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                   <  2:1910466/150‑1 (MQ=255)

                                                                                                                                    aaTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCt                   >  1:1932000/1‑150 (MQ=255)

                                                                                                                                     aTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTg                  <  1:76516/150‑1 (MQ=255)

                                                                                                                                     aTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTg                  >  2:692300/1‑150 (MQ=255)

                                                                                                                                     aTGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTg                  <  2:1107638/150‑1 (MQ=255)

                                                                                                                                       ggAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCACCGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGc                >  2:1364093/1‑150 (MQ=255)

                                                                                                                                         aaGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCg              <  2:856238/150‑1 (MQ=255)

                                                                                                                                           gCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGaa            >  2:1028865/1‑150 (MQ=255)

                                                                                                                                            gcactgACGGGTATACACGTAAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGaaa           >  1:1092004/5‑150 (MQ=255)

                                                                                                                                            ccTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGaaa           >  1:163869/1‑150 (MQ=255)

                                                                                                                                            ccTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGaaa           >  2:913069/1‑150 (MQ=255)

                                                                                                                                                tgACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatgat       <  2:34155/150‑1 (MQ=255)

                                                                                                                                                 gACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatgatg      <  2:1373631/150‑1 (MQ=255)

                                                                                                                                                 gACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAatgatg      <  1:27595/150‑1 (MQ=255)

TATCCAACAGCAACTGGCTGGCGAGCAGGTAGCGCTCTCTGACGAAGTCAACCAGTCGGAACAAACAACCAACTTCCATAACCGGGCCATAGCCTGGCTGCTGTACTCCGCCGGATATCTCTATTGTGATGCAATGGAAGCCTGTGACGTGTATACCCGTCAGTGCTCCACGCTCCTCAATACTATTGAACTGGCAACGCTTGGCGCGACGCTGGCGGCAGGTGGTGTGAATCCGTTGACGCATAAACGCGTTCTTCAGGCCGACAACGTGCCGTACATTCTGGCCGAAATGATGATGG  >  NC_000913/512063‑512361

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