BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,969,464	C→T	E5K (GAA→AAA)	ygdT ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,969,464	0	C	T	97.4%	1125.0 / 10.4	268	E5K (GAA→AAA)	ygdT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (2/5); new base T (128/133); total (130/138)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 276 total aligned reads displayed.

TAAGTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAA  >  NC_000913/2969366‑2969561

tAAGTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTtg                                                                                                 >  2:4509841/1‑101 (MQ=255)

  aGTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTtgtg                                                                                               >  2:4509844/1‑101 (MQ=255)

   gTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTgg                                                                                              >  1:4509845/1‑101 (MQ=255)

    ttAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGa                                                                                             >  1:4509850/1‑101 (MQ=255)

    ttAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGa                                                                                             >  2:4509848/1‑101 (MQ=255)

      aaGCACTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATa                                                                                           >  1:4509854/1‑101 (MQ=255)

      aaGCACTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATa                                                                                           <  2:4509523/101‑1 (MQ=255)

         cACTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACa                                                                                        >  2:4509857/1‑101 (MQ=255)

           cTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATg                                                                                      >  1:4509860/1‑101 (MQ=255)

           cTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATg                                                                                      <  1:4509281/101‑1 (MQ=255)

            tGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATgg                                                                                     >  2:4509864/1‑101 (MQ=255)

            tGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATgg                                                                                     <  1:4509408/101‑1 (MQ=255)

              aaCCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGgtg                                                                                   <  1:4509314/101‑1 (MQ=255)

               aCCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATgg                                                                                     <  2:4509148/98‑1 (MQ=255)

               aCCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGgtgt                                                                                  <  1:4509417/101‑1 (MQ=255)

                ccACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTa                                                                                 <  1:4509303/101‑1 (MQ=255)

                ccACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTa                                                                                 >  2:4509868/1‑101 (MQ=255)

                ccACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTa                                                                                 <  1:4509254/101‑1 (MQ=255)

                 cACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTaa                                                                                >  1:4509869/1‑101 (MQ=255)

                 cACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTaa                                                                                <  1:4509463/101‑1 (MQ=255)

                 cACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTaa                                                                                <  2:4509453/101‑1 (MQ=255)

                  aCTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAg                                                                               >  1:4509873/1‑101 (MQ=255)

                  aCTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAg                                                                               <  1:4509189/101‑1 (MQ=255)

                     aaaaaCTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACATGGTGTAAGATc                                                                            <  2:4509285/101‑1 (MQ=255)

                      aaaaCTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATcc                                                                           <  1:4509342/101‑1 (MQ=255)

                      aaaaCTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATcc                                                                           >  1:4509876/1‑101 (MQ=255)

                       aaaCTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCt                                                                          >  2:4509878/1‑101 (MQ=255)

                        aaCTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTg                                                                         >  2:4509880/1‑101 (MQ=255)

                         aCTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGt                                                                        <  1:4509339/101‑1 (MQ=255)

                            ggAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGa                                                                              <  2:4509289/92‑1 (MQ=255)

                            ggAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTa                                                                     <  1:4509329/101‑1 (MQ=255)

                             gAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTAt                                                                    <  2:4509228/101‑1 (MQ=255)

                             gAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTAt                                                                    <  2:4509328/101‑1 (MQ=255)

                              aGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTAtt                                                                   >  2:4509884/1‑101 (MQ=255)

                              aGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTAtt                                                                   >  1:4509886/1‑101 (MQ=255)

                                   cgtcgCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAg                                                              >  1:4509891/1‑101 (MQ=255)

                                   cgtcgCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAg                                                              >  2:4509890/1‑101 (MQ=255)

                                    gtcgCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGt                                                             >  1:4509896/1‑101 (MQ=255)

                                    gtcgCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGt                                                             <  2:4509402/101‑1 (MQ=255)

                                    gtcgCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGt                                                             >  2:4509895/1‑101 (MQ=255)

                                      cgCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGa                                                           >  2:4509898/1‑101 (MQ=255)

                                      cgCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGa                                                           <  2:4509433/101‑1 (MQ=255)

                                       gCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTAAGTGAc                                                          <  1:4509507/101‑1 (MQ=255)

                                        cACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGt                                                                        <  2:4509385/86‑1 (MQ=255)

                                        cACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGAcc                                                         >  2:4509901/1‑101 (MQ=255)

                                         aCGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCa                                                        >  2:4509904/1‑101 (MQ=255)

                                            tCAAGGCTGTAAATGGAAACAGTAGTTGAGGGTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGAc                                                          >  1:4509909/1‑96 (MQ=255)

                                            tCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGAt                                                     <  2:4509753/101‑1 (MQ=255)

                                             cAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGAtt                                                    >  2:4509910/1‑101 (MQ=255)

                                               aGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTg                                                  >  2:4509912/1‑101 (MQ=255)

                                                ggCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTgg                                                 <  2:4509312/101‑1 (MQ=255)

                                                  cTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGaa                                               >  2:4509916/1‑101 (MQ=255)

                                                   tGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGaaa                                              <  2:4509379/101‑1 (MQ=255)

                                                    gTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGaaaa                                             >  1:4509919/1‑101 (MQ=255)

                                                     tAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAAc                                            <  1:4508953/101‑1 (MQ=255)

                                                     tAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAAc                                            >  2:4509921/1‑101 (MQ=255)

                                                      aaaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAAcc                                           >  2:4509927/1‑101 (MQ=255)

                                                      aaaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAAcc                                           <  1:4509368/101‑1 (MQ=255)

                                                      aaaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAAcc                                           >  1:4509925/1‑101 (MQ=255)

                                                       aaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAAccc                                          <  2:4509372/101‑1 (MQ=255)

                                                       aaTGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAAccc                                          >  2:4509929/1‑101 (MQ=255)

                                                          ggAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGtt                                       <  2:4509337/101‑1 (MQ=255)

                                                          ggAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGtt                                       >  2:4509933/1‑101 (MQ=255)

                                                           gAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTg                                      >  1:4509934/1‑101 (MQ=255)

                                                               cAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAgt                                  >  2:4509938/1‑101 (MQ=255)

                                                                agtagtGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGaaa                                                                                           >  1:4509944/1‑41 (MQ=38)

                                                                agtagtGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAgtg                                 >  2:4509941/1‑101 (MQ=255)

                                                                agtagtGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAgtg                                 <  1:4509591/101‑1 (MQ=255)

                                                                agtagtGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAgtg                                 <  2:4509515/101‑1 (MQ=255)

                                                                  tagtGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGtt                               >  1:4509947/1‑101 (MQ=255)

                                                                   agtGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGtt                               <  2:4509251/100‑1 (MQ=255)

                                                                   agtGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTg                              <  1:4509532/101‑1 (MQ=255)

                                                                    gtGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTgc                             <  1:4509438/101‑1 (MQ=255)

                                                                     tGGAGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTgcg                            <  1:4509395/101‑1 (MQ=255)

                                                                        aGGTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCaa                         <  2:4509357/101‑1 (MQ=255)

                                                                         ggTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAAc                        >  1:4509951/1‑101 (MQ=255)

                                                                         ggTTTTTCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAAc                        >  2:4509950/1‑101 (MQ=255)

                                                                           tttttCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTc                      <  2:4509366/101‑1 (MQ=255)

                                                                           tttttCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTc                      <  2:4509470/101‑1 (MQ=255)

                                                                           tttttCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTc                      >  1:4509954/1‑101 (MQ=255)

                                                                            ttttCACAGTTATCCCAGCTTTTTGTGGATAACATGgtgca                                                                                 >  2:4509956/1‑39 (MQ=39)

                                                                             tttCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGa                    <  1:4509153/101‑1 (MQ=255)

                                                                             tttCACAGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGa                    <  1:4509521/101‑1 (MQ=255)

                                                                              ttCACAGTTATCCCAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGAt                   <  1:4509570/101‑1 (MQ=255)

                                                                                cacaGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTa                 >  2:4509960/1‑101 (MQ=255)

                                                                                 acaGTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTAc                >  1:4509962/1‑101 (MQ=255)

                                                                                    gTTATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCgg             >  1:4509967/1‑101 (MQ=255)

                                                                                     ttATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGc            <  1:4509403/101‑1 (MQ=255)

                                                                                     ttATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGc            <  2:4509108/101‑1 (MQ=255)

                                                                                      tATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCa           <  2:4509653/101‑1 (MQ=255)

                                                                                      tATCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCa           <  2:4509423/101‑1 (MQ=255)

                                                                                        tCCCAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAAc         >  2:4509968/1‑101 (MQ=255)

                                                                                           cAGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTa      <  2:4509499/101‑1 (MQ=255)

                                                                                            aGCTTTTTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTaa     <  1:4509387/101‑1 (MQ=255)

                                                                                            aGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTaa     <  1:4509655/101‑1 (MQ=255)

                                                                                               ttttTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTaaaaa  >  1:4509972/1‑101 (MQ=255)

                                                                                               ttttTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTaaaaa  <  1:4509482/101‑1 (MQ=255)

                                                                                               ttttTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTaaaaa  <  2:4509572/101‑1 (MQ=255)

                                                                                               ttttTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTaaaaa  <  2:4509600/101‑1 (MQ=255)

                                                                                               ttttTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTaaaaa  <  2:4509630/101‑1 (MQ=255)

TAAGTTAAGCACTGAACCACTAAAAACTGGAGTTTCGTCGCACGTCAAGGCTGTAAATGGAAACAGTAGTGGAGGTTTTTCACAGTTATCCCAGCTTTCTGTGGATAACATGGTGTAAGATCCTGTTTATTTTCAGTGACCAGATTTGGAAAACCCGTTGCAGTGTTGCGCAACTCGATTACCGGCAACCTAAAAA  >  NC_000913/2969366‑2969561

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