BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,796,375	T→A	6.7%	I138F (ATT→TTT)	ygaU ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,796,375	0	T	A	6.7%	420.3 / 16.2	120	I138F (ATT→TTT)	ygaU	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base T (57/55); new base A (2/6); total (59/61)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 126 total aligned reads displayed.

CCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTT  >  NC_000913/2796228‑2796523

ccTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCTCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAt                                                                                                                                                      <  3:919474/148‑1 (MQ=255)

  tGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAtttt                                                                                                                                                   >  2:447651/1‑149 (MQ=255)

  tGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAttt                                                                                                                                                    >  4:104795/1‑148 (MQ=255)

         ggcggTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCgg                                                                                                                                             <  4:112247/148‑1 (MQ=255)

         ggcggTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAAATTTATCCgg                                                                                                                                             <  3:963891/148‑1 (MQ=255)

          gcggTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCggg                                                                                                                                            >  4:52055/1‑148 (MQ=255)

            ggTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCt                                                                                                                                          <  3:52055/148‑1 (MQ=255)

             gTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCttt                                                                                                                                        >  2:573870/1‑149 (MQ=255)

               aCTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCttttt                                                                                                                                      <  2:998284/149‑1 (MQ=255)

                   gcgcGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCa                                                                                                                                  >  1:388455/1‑149 (MQ=255)

                     gcgAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCAt                                                                                                                                 <  3:288679/148‑1 (MQ=255)

                        agaTACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCgg                                                                                                                              <  4:163320/148‑1 (MQ=255)

                            aCCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTAt                                                                                                                         <  2:584247/149‑1 (MQ=255)

                              cgcgcgcgCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTc                                                                                                                       <  2:672790/149‑1 (MQ=255)

                              cgcgcgcgCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTc                                                                                                                       >  1:527366/1‑149 (MQ=255)

                                   gcgcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGATGGCAATTAACTATTCGTGCGTAATGCGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTc                                                                                                                  >  2:769073/1‑149 (MQ=255)

                                     gcCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGa                                                                                                                >  2:1055179/1‑149 (MQ=255)

                                      ccAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCACTTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGa                                                                                                                >  4:312675/1‑148 (MQ=255)

                                       cAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGaa                                                                                                               >  4:376461/1‑148 (MQ=255)

                                          tGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGAtt                                                                                                           <  1:1055241/149‑1 (MQ=255)

                                          tGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGAt                                                                                                            >  4:109620/1‑148 (MQ=255)

                                             gttcaggcgcCAGCATTTACTGACGCACATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGAtttt                                                                                                         <  4:683856/142‑1 (MQ=255)

                                                aaTGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTg                                                                                                     <  1:78528/149‑1 (MQ=255)

                                                 aTGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTg                                                                                                     >  3:173875/1‑148 (MQ=255)

                                                  tGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTa                                                                                                   <  2:782679/149‑1 (MQ=255)

                                                    cgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTAc                                                                                                  <  3:517622/148‑1 (MQ=255)

                                                    cgcCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTAc                                                                                                  <  3:376461/148‑1 (MQ=255)

                                                       cAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGAt                                                                                              >  1:949069/1‑149 (MQ=255)

                                                       cAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGAt                                                                                              <  1:96522/149‑1 (MQ=255)

                                                          cATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAg                                                                                           <  2:869826/149‑1 (MQ=255)

                                                           aTTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGTATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAg                                                                                           >  4:723646/1‑148 (MQ=255)

                                                            tttACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCg                                                                                         >  2:667058/1‑149 (MQ=255)

                                                                 tGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTCTCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTAcc                                                                                    <  2:185670/149‑1 (MQ=255)

                                                                   aCGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACCGATTAGCGTTACCg                                                                                   <  1:807103/148‑1 (MQ=255)

                                                                   aCGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGt                                                                                  >  2:68899/1‑149 (MQ=255)

                                                                    cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTa                                                                                 <  1:573870/149‑1 (MQ=255)

                                                                    cGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTa                                                                                 <  2:388455/149‑1 (MQ=255)

                                                                     gCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAg                                                                                >  2:566808/1‑149 (MQ=255)

                                                                      ccAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAg                                                                                >  4:776926/1‑148 (MQ=255)

                                                                      ccAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAg                                                                                <  3:35745/148‑1 (MQ=255)

                                                                            cccATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGtt                                                                         <  1:315567/149‑1 (MQ=255)

                                                                            cccATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGt                                                                          <  4:313947/148‑1 (MQ=255)

                                                                             ccATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGtt                                                                         <  3:312675/148‑1 (MQ=255)

                                                                              cATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTg                                                                       >  1:273960/1‑149 (MQ=255)

                                                                              cATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTg                                                                       >  1:414545/1‑149 (MQ=255)

                                                                               aTTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTg                                                                       <  3:104766/148‑1 (MQ=255)

                                                                               aTTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCCTTACCGTAGACCTGTTTg                                                                       <  3:645726/148‑1 (MQ=255)

                                                                                ttATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGa                                                                     >  1:435873/1‑149 (MQ=255)

                                                                                 tATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGaa                                                                    >  2:369756/1‑149 (MQ=255)

                                                                                   tCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGaaa                                                                   <  4:173875/148‑1 (MQ=255)

                                                                                   tCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGaaa                                                                   <  4:221727/148‑1 (MQ=255)

                                                                                     cagcagGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATgg                                                                <  1:667058/149‑1 (MQ=255)

                                                                                     cagcagGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATg                                                                 <  4:434091/148‑1 (MQ=255)

                                                                                      agcagGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATgg                                                                >  3:383891/1‑148 (MQ=255)

                                                                                              aTGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCaga                                                        >  4:861409/1‑148 (MQ=255)

                                                                                                  ctggcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTc                                                   <  1:340683/149‑1 (MQ=255)

                                                                                                     gcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGcc                                                <  1:20454/149‑1 (MQ=255)

                                                                                                     gcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGcc                                                <  1:424622/149‑1 (MQ=255)

                                                                                                     gcAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGc                                                 <  3:743500/148‑1 (MQ=255)

                                                                                                        aTTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAg                                              >  3:625121/1‑148 (MQ=255)

                                                                                                         ttAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGa                                             >  3:635038/1‑148 (MQ=255)

                                                                                                         ttAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGAc                                            <  1:617016/149‑1 (MQ=255)

                                                                                                              tactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTaa                                        >  3:862255/1‑148 (MQ=255)

                                                                                                              tactCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAAc                                       >  2:594189/1‑149 (MQ=255)

                                                                                                              tactCTCCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTaa                                        >  4:181636/1‑148 (MQ=255)

                                                                                                               actCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACg                                      >  1:777433/1‑149 (MQ=255)

                                                                                                                  ctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGCCCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGta                                   >  2:1009947/1‑149 (MQ=255)

                                                                                                                  ctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGa                                             >  3:540678/1‑139 (MQ=255)

                                                                                                                  ctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGa                                             <  4:540678/139‑1 (MQ=255)

                                                                                                                  ctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGt                                    <  3:546970/148‑1 (MQ=255)

                                                                                                                  ctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGt                                    <  3:861409/148‑1 (MQ=255)

                                                                                                                  ctTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGt                                    >  4:785339/1‑148 (MQ=255)

                                                                                                                   ttCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGtat                                  >  1:351967/1‑149 (MQ=255)

                                                                                                                   ttCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGtat                                  >  1:722529/1‑149 (MQ=255)

                                                                                                                          aTACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActg                           <  1:369756/149‑1 (MQ=255)

                                                                                                                          aTACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAAct                            >  3:690713/1‑148 (MQ=255)

                                                                                                                            aCGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActggc                         >  1:529539/1‑149 (MQ=255)

                                                                                                                            aCGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActggc                         <  2:749013/149‑1 (MQ=255)

                                                                                                                             cGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActggct                        >  1:1055638/1‑149 (MQ=255)

                                                                                                                              gCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActggctg                       >  1:964442/1‑149 (MQ=255)

                                                                                                                              gCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActggctg                       >  2:688684/1‑149 (MQ=255)

                                                                                                                              gCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAActggctg                       >  1:684209/1‑149 (MQ=255)

                                                                                                                               cAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCg                                                                                                                      <  3:1002778/53‑1 (MQ=255)

                                                                                                                               cAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCg                                                                                                                      >  4:1002778/1‑53 (MQ=255)

                                                                                                                                    cTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATaaa                               <  4:314107/135‑1 (MQ=255)

                                                                                                                                        ccccGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGGACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACACAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCt              <  4:257138/148‑1 (MQ=255)

                                                                                                                                          ccGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTgg            <  3:199193/148‑1 (MQ=255)

                                                                                                                                          ccGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTgg            <  4:690713/148‑1 (MQ=255)

                                                                                                                                          ccGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTgct           <  2:435873/149‑3 (MQ=255)

                                                                                                                                          ccGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGgt           <  1:447985/149‑1 (MQ=255)

                                                                                                                                          ccGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGgt           >  2:203517/1‑149 (MQ=255)

                                                                                                                                          ccGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGgt           >  2:320713/1‑149 (MQ=255)

                                                                                                                                          ccGGATAAAATTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTAtt                                                                                                      <  1:636279/58‑1 (MQ=255)

                                                                                                                                           cGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTgg            >  3:756831/1‑147 (MQ=255)

                                                                                                                                           cGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTgg            <  4:756831/147‑1 (MQ=255)

                                                                                                                                           cGGATAAAATTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGgt           >  4:1026657/1‑148 (MQ=255)

                                                                                                                                            ggATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGgtg          >  3:319822/1‑148 (MQ=255)

                                                                                                                                              aTAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCg       <  1:688684/149‑1 (MQ=255)

                                                                                                                                               tAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGc      <  1:713697/149‑1 (MQ=255)

                                                                                                                                                   ttttATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGtt  >  1:409357/1‑149 (MQ=255)

CCTGACTGCGGCGGTACTGGCGCGAGATACCGCGCGCGCCAGTGAGTTAATGCGCCAGCATTTACTGACGCCAATTCCCATTATCCAGCAGGCGATGGCTGGCAATTAACTACTCTTCCGGAATACGCAACACTTGCCCCGGATAAATTTTATCCGGGCTTTTTAGCATCGGTTTATTCGCTTCGAAGATTTTATTGTACAGATTAGCGTTACCGTAGACCTGTTTGGAAATGGCACTCAGAGTGTCGCCAGACTTAACGGTATAAAACTGGCTGGCAGTGGCTGGTGTCGCCGTT  >  NC_000913/2796228‑2796523

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