BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,388,343	T→C	T249A (ACC→GCC)	b4161 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,388,343	0	T	C	99.4%	1022.3 / ‑4.6	332	T249A (ACC→GCC)	b4161	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (153/177); total (155/177)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 360 total aligned reads displayed.

CCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTCAGCAGGCTGGATTTGC  >  NC_000913/4388253‑4388432

ccAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGgcgt                                                                                         >  1:1412663/1‑93 (MQ=255)

  aGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGc                                                                                       <  2:3829794/93‑1 (MQ=255)

   ggTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCt                                                                                      <  2:1962817/93‑1 (MQ=255)

   ggTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCt                                                                                      <  2:2233503/93‑1 (MQ=255)

    gTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGCCGTGGCGTGCTg                                                                                     >  2:880140/1‑93 (MQ=255)

     tGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGa                                                                                    >  1:7240608/1‑93 (MQ=255)

        cAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCg                                                                                 >  1:6857840/1‑93 (MQ=255)

        cAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCg                                                                                 >  2:5410797/1‑93 (MQ=255)

        cAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCg                                                                                 >  2:3439508/1‑93 (MQ=255)

        cAGAGGCCGAACTCACGCACTCCTGGGGAGTCAAGCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCg                                                                                 <  2:4235341/93‑1 (MQ=255)

           agGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGc                                                                              >  1:4638127/1‑93 (MQ=255)

           agGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGc                                                                              <  1:374549/93‑1 (MQ=255)

            ggCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGTGTGCTGACCGAGcc                                                                             >  1:3574225/1‑93 (MQ=255)

            ggCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGcc                                                                             <  1:2109112/93‑1 (MQ=255)

             gCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGccc                                                                            <  2:7711758/93‑1 (MQ=255)

              ccGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGcccc                                                                           <  2:5541869/93‑1 (MQ=255)

               cGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCg                                                                          >  1:5204512/1‑93 (MQ=255)

                gAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCgg                                                                         >  2:7106920/1‑92 (MQ=255)

                gAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGa                                                                         <  1:829770/93‑1 (MQ=255)

                gAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGa                                                                         <  1:5338848/93‑1 (MQ=255)

                 aaCTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCgcg                                                                        >  1:6014783/1‑91 (MQ=255)

                 aaCTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAg                                                                        >  2:1893221/1‑93 (MQ=255)

                  aCTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGt                                                                       >  1:3078285/1‑93 (MQ=255)

                  aCTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGt                                                                       <  2:2438800/93‑1 (MQ=255)

                   cTCCCGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGtt                                                                      <  1:1271760/93‑1 (MQ=255)

                   cTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGtt                                                                      >  2:2662734/1‑93 (MQ=255)

                    tCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTg                                                                     <  2:3608734/93‑1 (MQ=255)

                       cGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCa                                                                  <  2:6721278/93‑1 (MQ=255)

                       cGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCa                                                                  <  1:876336/93‑1 (MQ=255)

                        gCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCag                                                                 >  2:145454/1‑93 (MQ=255)

                          aCTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCagag                                                               >  1:1074776/1‑93 (MQ=255)

                          aCTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCagag                                                               >  1:5452448/1‑93 (MQ=255)

                           cTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAgaga                                                              <  1:4314986/93‑1 (MQ=255)

                           cTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAgaga                                                              >  1:7446241/1‑93 (MQ=255)

                            tCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGat                                                             >  2:3018883/1‑93 (MQ=255)

                            tCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGat                                                             >  1:6589369/1‑93 (MQ=255)

                            tCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGat                                                             <  1:3098410/93‑1 (MQ=255)

                             ccTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGata                                                            <  1:6253363/93‑1 (MQ=255)

                             ccTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGata                                                            <  2:474014/93‑1 (MQ=255)

                             ccTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGata                                                            >  1:5832634/1‑93 (MQ=255)

                             ccTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGata                                                            <  2:7418107/93‑1 (MQ=255)

                             ccTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGata                                                            <  2:6570829/93‑1 (MQ=255)

                              cTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGatat                                                           <  1:3342707/93‑1 (MQ=255)

                               tGGGGAGTCAATCACATCCCCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATc                                                          <  2:5964396/93‑1 (MQ=255)

                               tGGGGAGTCAATCACATCACCGCTGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATc                                                          <  2:5361312/93‑1 (MQ=255)

                               tGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATc                                                          <  2:159245/93‑1 (MQ=255)

                                ggggAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCg                                                         <  1:1742931/93‑1 (MQ=255)

                                  ggAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGtt                                                       <  2:7935206/93‑1 (MQ=255)

                                  ggAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGtt                                                       >  2:6041694/1‑93 (MQ=255)

                                  ggAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGtt                                                       <  1:7394676/93‑1 (MQ=255)

                                    aGTCAATCACACCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTgg                                                     <  1:4079806/93‑1 (MQ=255)

                                     gTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGGTATCGTTGGt                                                    >  2:6904109/1‑93 (MQ=255)

                                      tCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTc                                                   >  1:2441111/1‑93 (MQ=255)

                                      tCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGGTGTCAGAGATATCGTTGGTc                                                   >  1:7166586/1‑93 (MQ=255)

                                        aaTCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAg                                                 >  1:3217739/1‑93 (MQ=255)

                                        aaTCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAg                                                 >  2:7867983/1‑93 (MQ=255)

                                        aaTCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAg                                                 <  2:5190352/93‑1 (MQ=255)

                                        aaTCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAg                                                 <  2:3254549/93‑1 (MQ=255)

                                        aaTCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAg                                                 >  2:7282306/1‑93 (MQ=255)

                                         atcacatcacCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCCTGGCGTGCTGACGGAGCCCGGCGGTGTCAGAGATATCTTTGGGCAga                                                >  1:2695235/1‑92 (MQ=255)

                                             catcacCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGAtct                                            >  2:2332621/1‑93 (MQ=255)

                                             catcacCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGAtct                                            >  2:6102066/1‑93 (MQ=255)

                                              atcacCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGGTGGTCAGGAtctt                                           >  2:4598493/1‑92 (MQ=255)

                                                cacCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCtt                                         >  1:6506210/1‑93 (MQ=255)

                                                 acCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCttt                                        >  2:5258211/1‑93 (MQ=255)

                                                  ccgccgTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCtttt                                       <  1:6310237/93‑1 (MQ=255)

                                                   cgccgTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCttttt                                      <  2:3972020/93‑1 (MQ=255)

                                                    gccgTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTg                                     <  2:7736738/93‑1 (MQ=255)

                                                     ccgTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGc                                    <  2:570144/93‑1 (MQ=255)

                                                      cgTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCa                                   <  1:2227406/93‑1 (MQ=255)

                                                      cgTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCa                                   <  2:3108008/93‑1 (MQ=255)

                                                      cgTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCa                                   <  1:7360946/93‑1 (MQ=255)

                                                       gTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAg                                  <  2:2801603/93‑1 (MQ=255)

                                                        tGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGGAGc                                 >  1:94570/1‑93 (MQ=255)

                                                           gggAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGcccc                              <  2:6273836/93‑1 (MQ=255)

                                                            ggccgTGCTACGGCCGAGCGGCGGTCGTGGCGTGCTGCCCGAGCCCCGAGTTGTCGGAGGTCTCGTTGGTCAGGATCTCTTTTTGCCGCCCCa                             <  1:6998488/89‑1 (MQ=255)

                                                            ggAAGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCa                             <  1:4963307/93‑1 (MQ=255)

                                                              aaGTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGc                           >  2:4257866/1‑93 (MQ=255)

                                                                gTGATACAGCGGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCaa                         <  1:4343869/93‑1 (MQ=255)

                                                                gTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCaa                         <  2:1209443/93‑1 (MQ=255)

                                                                gTGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCaa                         <  2:5134681/93‑1 (MQ=255)

                                                                 tGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGtt                                                       >  1:7902855/1‑62 (MQ=255)

                                                                 tGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAAc                        >  2:2155564/1‑93 (MQ=255)

                                                                 tGATACAGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAAc                        <  2:2549798/93‑1 (MQ=255)

                                                                       aGCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTc                  <  2:3913001/93‑1 (MQ=255)

                                                                        gCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTca                 >  1:4197985/1‑93 (MQ=255)

                                                                        gCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGAAGCCCCCGCAACGCATTca                 >  1:6780979/1‑93 (MQ=255)

                                                                        gCCGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTAGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTca                 <  1:1927418/93‑1 (MQ=255)

                                                                          cGAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTcagc               >  1:2866010/1‑93 (MQ=255)

                                                                           gAGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTcagca              >  2:1091343/1‑93 (MQ=255)

                                                                           gAGCGGCGGTCGGGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTcagca              <  1:1626052/93‑1 (MQ=255)

                                                                            aGCGGCGGTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTcagcag             <  1:4827127/93‑1 (MQ=255)

                                                                               ggcggTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTCAGCAGGCt          >  2:1712717/1‑93 (MQ=255)

                                                                               ggcggTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTCAGCAGGCt          <  2:4937818/93‑1 (MQ=255)

                                                                                gcggTCTTGGCGTGCTGACCGAGCCACGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTCAGCAGGCTg         <  1:7791905/93‑1 (MQ=255)

                                                                                gcggTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTCAGCAGGCTg         >  1:3881463/1‑93 (MQ=255)

                                                                                  ggTCGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTCAGCAGGCTGGa       <  2:6391165/93‑1 (MQ=255)

                                                                                     cGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTCAGCAGGCTGGAttt    <  1:7677267/93‑1 (MQ=255)

                                                                                     cGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTCAGCAGGCTGGAttt    <  1:5675236/93‑1 (MQ=255)

                                                                                     cGTGGCGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTCAGCAGGCTGGAttt    <  1:217594/93‑1 (MQ=255)

CCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCAGAGATATCGTTGGTCAGGATCTCTTTTTGCAGCCCCAGCAACGCATTCAGCAGGCTGGATTTGC  >  NC_000913/4388253‑4388432

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