BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,678,204	G→A	P456L (CCG→CTG)	yhjJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,678,204	0	G	A	98.4%	673.4 / ‑2.4	190	P456L (CCG→CTG)	yhjJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/2); new base A (84/103); total (85/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 192 total aligned reads displayed.

TCATCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCGGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGCTA  >  NC_000913/3678114‑3678293

tCATCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGc                                                                                         >  2:3442348/1‑93 (MQ=255)

tCATCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGc                                                                                         >  1:32799/1‑93 (MQ=255)

 cATCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCt                                                                                        >  1:825138/1‑93 (MQ=255)

  aTCCGTGGCGACGGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTc                                                                                       >  2:2972727/1‑93 (MQ=255)

  aTCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTc                                                                                       <  1:241642/93‑1 (MQ=255)

   tCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATTATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTcc                                                                                      >  1:1789760/1‑92 (MQ=255)

   tCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCg                                                                                      >  1:1354957/1‑93 (MQ=255)

     gggggCGAGGGAGGTGGTCGCGGCGGCGGTAGAGGGGGCCATGATTTGACCCCGGCCCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGcc                                                                                    <  2:3117925/90‑1 (MQ=255)

     cGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGcc                                                                                    <  2:1438014/93‑1 (MQ=255)

       tGGCGACAGAGGTGGTCGCGGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCtt                                                                                  >  1:4599314/1‑93 (MQ=255)

         gCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTc                                                                                <  2:2618956/93‑1 (MQ=255)

         gCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTc                                                                                <  1:1985682/93‑1 (MQ=255)

         gCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTc                                                                                <  1:1985683/93‑1 (MQ=255)

         gCGACAAAGGTGGACACAACGGCGGTAGATGGGGCCATGATTGGATCCCAGACAGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTc                                                                                <  2:2457784/93‑1 (MQ=255)

            aCAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCGGCTCGCCTTTCGGc                                                                             <  2:4236776/93‑1 (MQ=255)

             cAGAGGTGGTCGCAGCGGCGGAAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCt                                                                            <  2:2226833/93‑1 (MQ=255)

              agagGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTg                                                                           >  2:3969504/1‑93 (MQ=255)

               gagGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgc                                                                          >  2:1310262/1‑93 (MQ=255)

               gagGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgc                                                                          >  1:4742837/1‑93 (MQ=255)

               gagGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgc                                                                          <  2:3925794/93‑1 (MQ=255)

                agGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgca                                                                         <  2:4327665/93‑1 (MQ=255)

                agGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgca                                                                         <  2:2832202/93‑1 (MQ=255)

                agGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgca                                                                         <  1:1943474/93‑1 (MQ=255)

                   tggtCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgcagca                                                                      <  2:4484211/93‑1 (MQ=255)

                   tggtCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTgcagca                                                                      <  1:692445/93‑1 (MQ=255)

                    ggtCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGcagcag                                                                     <  1:4365412/93‑1 (MQ=255)

                    ggtCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGcagcag                                                                     <  1:3832120/93‑1 (MQ=255)

                        gCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTAtt                                                                 <  1:2260552/93‑1 (MQ=255)

                          aGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTaa                                                               >  2:1713908/1‑93 (MQ=255)

                          aGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTaa                                                               >  2:2348720/1‑93 (MQ=255)

                              gcggTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCGGCTCGCCTTTCGGCTGCAGCAGTATTAACGcc                                                           >  2:1554108/1‑93 (MQ=255)

                              gcggTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGcc                                                           <  2:3455129/93‑1 (MQ=255)

                               cggTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCa                                                          <  1:3401910/93‑1 (MQ=255)

                                 gTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCata                                                        <  2:896590/93‑1 (MQ=255)

                                 gTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCata                                                        <  2:3003320/93‑1 (MQ=255)

                                 gTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCata                                                        >  1:2242704/1‑93 (MQ=255)

                                 gTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCata                                                        <  1:3082244/93‑1 (MQ=255)

                                   aGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATc                                                      <  1:1701445/93‑1 (MQ=255)

                                    gATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCa                                                     >  2:768010/1‑93 (MQ=255)

                                      tGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCGGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCAtt                                                   <  2:2147262/93‑1 (MQ=255)

                                      tGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCAtt                                                   <  2:1855333/93‑1 (MQ=255)

                                      tGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCAtt                                                   <  2:4613700/93‑1 (MQ=255)

                                       ggggCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTc                                                  <  1:2820480/93‑1 (MQ=255)

                                         ggCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGa                                                <  1:1331116/93‑1 (MQ=255)

                                             aTGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAAc                                            >  2:3519623/1‑93 (MQ=255)

                                             aTGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAAc                                            >  2:1510864/1‑93 (MQ=255)

                                                aTTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGt                                         <  1:2098/93‑1 (MQ=255)

                                                 tttGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGtt                                        <  1:41895/93‑1 (MQ=255)

                                                   tGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGa                                      <  1:3805006/93‑1 (MQ=255)

                                                    gATCCCAGCCCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGAc                                     <  2:1393821/93‑1 (MQ=255)

                                                    gATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGAc                                     <  1:4323393/93‑1 (MQ=255)

                                                    gATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGAc                                     <  2:4727035/93‑1 (MQ=255)

                                                      tCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGc                                   <  1:4409040/93‑1 (MQ=255)

                                                      tCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGc                                   >  1:3208873/1‑93 (MQ=255)

                                                        ccAGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAg                                 >  1:1244313/1‑93 (MQ=255)

                                                          aGACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGAt                               <  2:832195/93‑1 (MQ=255)

                                                           gACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATc                              >  1:1527288/1‑93 (MQ=255)

                                                           gACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATc                              >  1:3981816/1‑93 (MQ=255)

                                                           gACCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATc                              >  1:4100398/1‑93 (MQ=255)

                                                            aCCGCCTGCAACGCTTTCATGTTACATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATcc                             <  1:4209812/93‑1 (MQ=255)

                                                            aCCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATcc                             <  1:1062524/93‑1 (MQ=255)

                                                            aCCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATcc                             >  1:1900220/1‑93 (MQ=255)

                                                            aCCGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATcc                             >  1:2907461/1‑93 (MQ=255)

                                                             ccGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCt                            <  2:837196/93‑1 (MQ=255)

                                                             ccGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCt                            <  2:4529118/93‑1 (MQ=255)

                                                             ccGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCt                            <  1:4238275/93‑1 (MQ=255)

                                                             ccGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCt                            <  2:2563097/93‑1 (MQ=255)

                                                             ccGCCTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCt                            <  1:1017901/93‑1 (MQ=255)

                                                             ccGCCTGCAACGCTTGCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCt                            <  1:4344281/93‑1 (MQ=255)

                                                                ccTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGAt                         <  2:1178136/93‑1 (MQ=255)

                                                                ccTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGAt                         <  2:2995740/93‑1 (MQ=255)

                                                                 cTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGAtt                        >  1:4583765/1‑93 (MQ=255)

                                                                 cTGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGAtt                        <  2:3277270/93‑1 (MQ=255)

                                                                  tGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGAttt                       >  2:2504360/1‑93 (MQ=255)

                                                                  tGCAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGAttt                       >  2:1954496/1‑93 (MQ=255)

                                                                    cAACGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTaa                     <  2:29082/93‑1 (MQ=255)

                                                                     aaCGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAAc                    <  2:693760/93‑1 (MQ=255)

                                                                      aCGCTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACa                   >  2:3790722/1‑93 (MQ=255)

                                                                         cTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtct                <  2:211315/93‑1 (MQ=255)

                                                                         cTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtct                <  1:2716599/93‑1 (MQ=255)

                                                                         cTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtct                >  1:1313651/1‑93 (MQ=255)

                                                                         cTTTCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTCACGCAGATCCTGATTTAACAtct                >  2:3212956/1‑93 (MQ=255)

                                                                          tttCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtctc               <  1:748900/93‑1 (MQ=255)

                                                                          tttCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtctc               <  1:315195/93‑1 (MQ=255)

                                                                          tttCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtctc               >  1:4476411/1‑93 (MQ=255)

                                                                          tttCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGAAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACAtctc               <  2:4356056/93‑1 (MQ=255)

                                                                           ttCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATctct              <  2:1915215/93‑1 (MQ=255)

                                                                           ttCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATctct              >  1:3739616/1‑93 (MQ=255)

                                                                            tCATGTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTa             <  2:868474/93‑1 (MQ=255)

                                                                                tttAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGt         <  2:299296/92‑1 (MQ=255)

                                                                                gTTAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGAACACTGTTGACGCAGATCCTGATTTAACATCTCTACCGt         >  1:1891567/1‑93 (MQ=255)

                                                                                  tAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTTACGCAGATCCTGATTTAACATCTCTACCGTCa       >  2:811379/1‑93 (MQ=255)

                                                                                  tAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCa       <  1:3074782/93‑1 (MQ=255)

                                                                                  tAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCa       <  2:3414616/93‑1 (MQ=255)

                                                                                  tAAATTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCa       <  1:3477027/93‑1 (MQ=255)

                                                                                   aaaTTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAg      >  2:4108234/1‑93 (MQ=255)

                                                                                     aTTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGc    >  2:2764398/1‑93 (MQ=255)

                                                                                     aTTCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGc    <  2:583788/93‑1 (MQ=255)

                                                                                      ttCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGCt   <  2:3740053/93‑1 (MQ=255)

                                                                                       tCCAGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGCTa  >  1:948838/1‑93 (MQ=255)

TCATCCGTGGCGACAGAGGTGGTCGCAGCGGCGGTAGATGGGGCCATGATTTGATCCCAGACCGCCTGCAACGCTTTCATGTTAAATTCCGGCTCGCCTTTCGGCTGCAGCAGTATTAACGCCATATCATTCGACAACTGTTGACGCAGATCCTGATTTAACATCTCTACCGTCAGGCTA  >  NC_000913/3678114‑3678293

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