BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,778,879	T→C	7.1%	V494A (GTC→GCC)	lldP →	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,778,879	0	T	C	7.1%	445.1 / 17.2	141	V494A (GTC→GCC)	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand): ref base T (58/72); new base C (5/5); total (63/78)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 143 total aligned reads displayed.

TCGCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATCACCACGCTTCAGGCT  >  NC_000913/3778734‑3779027

tcGCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCAc                                                                                                                                                   <  3:912226/149‑1 (MQ=255)

tcGCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCAc                                                                                                                                                   >  4:281245/1‑149 (MQ=255)

  gCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCg                                                                                                                                                 <  3:463945/149‑1 (MQ=255)

  gCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCg                                                                                                                                                 >  3:397154/1‑149 (MQ=255)

  gCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCg                                                                                                                                                 <  3:93101/149‑1 (MQ=255)

     gTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTa                                                                                                                                              >  1:855528/1‑149 (MQ=255)

     gTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGACGCCACCGGtg                                                                                                                                              <  1:56201/149‑2 (MQ=255)

       tCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCACCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGACGCCACCGGtggt                                                                                                                                            <  4:62376/149‑4 (MQ=255)

        ccTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAga                                                                                                                                           <  2:135543/149‑1 (MQ=255)

         cTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAgat                                                                                                                                          >  4:845371/1‑149 (MQ=255)

            ggctggctGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAgatgat                                                                                                                                       <  3:14509/149‑1 (MQ=255)

                ggctggGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTcc                                                                                                                                   <  4:363261/149‑1 (MQ=255)

                  ctggGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTcccc                                                                                                                                 <  3:180427/149‑1 (MQ=255)

                    gggggTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGc                                                                                                                               <  1:223735/149‑1 (MQ=255)

                           ttCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATctat                                                                                                                        <  1:47193/149‑1 (MQ=255)

                            tCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATctatc                                                                                                                       <  2:650618/149‑1 (MQ=255)

                            tCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCACGCAATctatc                                                                                                                       <  1:694774/149‑1 (MQ=255)

                               tGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATctatcgct                                                                                                                    >  4:421718/1‑149 (MQ=255)

                                   cGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATctatcgctatcg                                                                                                                >  3:193474/1‑149 (MQ=255)

                                    gggTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCtatcgctatcgc                                                                                                               <  1:198384/149‑1 (MQ=255)

                                      gTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCt                                                                                                             <  2:52698/149‑1 (MQ=255)

                                      gTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCACAATCTATCGCTATCGCCt                                                                                                             <  4:110950/149‑1 (MQ=255)

                                          gATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTgc                                                                                                         <  4:879107/149‑1 (MQ=255)

                                             aCCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTgcggc                                                                                                      <  1:206824/149‑1 (MQ=255)

                                              ccTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTACCCGCAATCTATCGCTATCGCCTGTgcggcg                                                                                                     >  4:553260/1‑149 (MQ=255)

                                               cTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGcggcgg                                                                                                    >  1:638330/1‑149 (MQ=255)

                                               cTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGcggcgg                                                                                                    <  4:114609/149‑1 (MQ=255)

                                               cTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGcggcgg                                                                                                    >  4:16839/1‑149 (MQ=255)

                                                tCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGt                                                                                                   >  4:536496/1‑149 (MQ=255)

                                                tCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGt                                                                                                   >  1:212842/1‑149 (MQ=255)

                                                tCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGt                                                                                                   >  2:787858/1‑149 (MQ=255)

                                                tCATCTAACGCCCTGTTAGCCGCGCTGCAAGCCACCGCAGCGCCACAAACTGGCGTCTGTGATCTGTTGATGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCACCCTGCAATCTATCGCTATCGCGTGTGCGGCAGt                                                                                                   >  4:570517/1‑149 (MQ=255)

                                                 cATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTa                                                                                                  <  2:148639/149‑1 (MQ=255)

                                                 cATCTAACGCCCTGTTCGCCGCGCTGCAAGACACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTa                                                                                                  >  4:423871/1‑149 (MQ=255)

                                                  aTCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAg                                                                                                 <  1:105296/149‑1 (MQ=255)

                                                      aaCGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCt                                                                                             >  1:145347/1‑149 (MQ=255)

                                                         gCCCTGTTCGCTGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggt                                                                                          >  1:694178/1‑149 (MQ=255)

                                                         gCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggt                                                                                          <  3:43495/149‑1 (MQ=255)

                                                         gCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggt                                                                                          <  2:785184/149‑1 (MQ=255)

                                                          cccTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggtg                                                                                         >  4:502871/1‑149 (MQ=255)

                                                           cttGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggtgg                                                                                        >  1:215884/3‑149 (MQ=255)

                                                           ccTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggtgg                                                                                        <  2:145347/149‑1 (MQ=255)

                                                           ccTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggtgg                                                                                        <  3:154233/149‑1 (MQ=255)

                                                           ccTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggtgg                                                                                        >  1:385739/1‑149 (MQ=255)

                                                           ccTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTATGCGGCGGTAGGCCtggtgg                                                                                        <  3:436588/149‑1 (MQ=255)

                                                            cTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTggg                                                                                       <  3:823749/149‑1 (MQ=255)

                                                                         tGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGAt                                                                          >  2:17249/1‑149 (MQ=255)

                                                                              gCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGtt                                                                     <  3:405988/149‑1 (MQ=255)

                                                                               ccACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTc                                                                    >  2:858386/1‑149 (MQ=255)

                                                                                  ccGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGc                                                                 <  1:476064/149‑1 (MQ=255)

                                                                                    gcagcaCAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCtt                                                               >  2:231684/1‑149 (MQ=255)

                                                                                       gcaCAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTAc                                                            <  2:5262/149‑1 (MQ=255)

                                                                                          caacaaATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGt                                                         <  2:636469/149‑1 (MQ=255)

                                                                                           aacaaATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTc                                                        >  2:760400/1‑149 (MQ=255)

                                                                                            acaaATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCa                                                       >  3:564828/1‑149 (MQ=255)

                                                                                            acaaATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCa                                                       <  3:261634/149‑1 (MQ=255)

                                                                                              aaaTTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCaaa                                                     <  4:757825/149‑1 (MQ=255)

                                                                                              aaaTTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCaaa                                                     <  4:538239/149‑1 (MQ=255)

                                                                                              aaaTTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCaaa                                                     <  1:528276/149‑1 (MQ=255)

                                                                                                 ttGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAacac                                                  >  2:474890/1‑149 (MQ=255)

                                                                                                 ttGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAacac                                                  <  3:281245/149‑1 (MQ=255)

                                                                                                   ggCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAg                                                <  1:17249/149‑1 (MQ=255)

                                                                                                    gCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGa                                                                           >  4:514116/1‑121 (MQ=255)

                                                                                                    gCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGa                                                                           <  3:514116/121‑1 (MQ=255)

                                                                                                      gTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATctatc                                                                                                                       >  2:869152/1‑75 (MQ=255)

                                                                                                        ctctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGa                                           <  2:67418/149‑1 (MQ=255)

                                                                                                         tctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGAt                                          >  4:557171/1‑149 (MQ=255)

                                                                                                         tctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGAt                                          >  1:122948/1‑149 (MQ=255)

                                                                                                         tctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGAt                                          >  1:345792/1‑149 (MQ=255)

                                                                                                         tctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGAt                                          <  2:215884/149‑1 (MQ=255)

                                                                                                          ctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACTGCCTGATc                                         <  4:149955/149‑1 (MQ=255)

                                                                                                          ctGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTACGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATc                                         <  2:919742/149‑1 (MQ=255)

                                                                                                           tGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCt                                        >  3:108272/1‑149 (MQ=255)

                                                                                                           tGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCt                                        >  1:778681/1‑149 (MQ=255)

                                                                                                           tGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCt                                        <  2:385739/149‑1 (MQ=255)

                                                                                                             aTCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTc                                      >  1:501395/1‑149 (MQ=255)

                                                                                                               cTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCaaa                                                                                   <  4:265633/102‑1 (MQ=255)

                                                                                                               cTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCAc                                    <  1:739947/149‑1 (MQ=255)

                                                                                                                 gTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCt                                  >  2:326060/1‑149 (MQ=255)

                                                                                                                 gTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCt                                  >  2:356711/1‑149 (MQ=255)

                                                                                                                        gTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTc                                                        <  1:609828/120‑1 (MQ=255)

                                                                                                                        gTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTc                                                        >  2:609828/1‑120 (MQ=255)

                                                                                                                        gTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGt                           >  3:377965/1‑149 (MQ=255)

                                                                                                                         ttGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTg                          <  1:173047/149‑1 (MQ=255)

                                                                                                                           gccgccAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTggg                        <  2:501395/149‑1 (MQ=255)

                                                                                                                              gccAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGt                     <  1:858386/149‑1 (MQ=255)

                                                                                                                                  aTACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATc                 <  2:168744/149‑1 (MQ=255)

                                                                                                                                    cttaccGGTGACGCCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATcac               >  3:270966/4‑149 (MQ=255)

                                                                                                                                    accaccGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATcac               <  3:536496/149‑1 (MQ=255)

                                                                                                                                     ttaccGGTGACGCCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTg                                            <  2:46502/117‑1 (MQ=255)

                                                                                                                                     ttaccGGTGACGCCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTg                                            >  1:46502/3‑119 (MQ=255)

                                                                                                                                      taccGGTGACGCCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATcacca             >  4:814969/2‑149 (MQ=255)

                                                                                                                                      caccGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCtggtgg                                                                                        >  1:274535/1‑74 (MQ=255)

                                                                                                                                       accGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTTAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATcaccac            >  3:587796/1‑149 (MQ=255)

                                                                                                                                       accGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATcaccac            <  2:745309/149‑1 (MQ=255)

                                                                                                                                       accGGTGACGCCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATcaccac            >  2:403424/1‑149 (MQ=255)

                                                                                                                                            tGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATCACCACGCTTc       <  1:309170/149‑1 (MQ=255)

                                                                                                                                             ggCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATCACCACGCTTCa      <  4:610632/149‑1 (MQ=255)

                                                                                                                                               cGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTTACACACAGCCTGATCTACACCTGTATAGTGGGCGTGATCACCACGCTTCAgg    <  2:287783/149‑1 (MQ=255)

                                                                                                                                                 tCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGCTTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCCGTATAGTGGGCGTGATCACCACGCTTCAGGCt  >  2:736534/1‑149 (MQ=255)

TCGCCGTTCCTCGGCTGGCTGGGGGTATTCCTGACCGGGTCGGATACCTCATCTAACGCCCTGTTCGCCGCGCTGCAAGCCACCGCAGCACAACAAATTGGCGTCTCTGATCTGTTGCTGGTTGCCGCCAATACCACCGGTGGCGTCACCGGTAAGATGATCTCCCCGCAATCTATCGCTATCGCCTGTGCGGCGGTAGGCCTGGTGGGCAAAGAGTCTGATTTGTTCCGCTTTACTGTCAAACACAGCCTGATCTTCACCTGTATAGTGGGCGTGATCACCACGCTTCAGGCT  >  NC_000913/3778734‑3779027

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: `‑`